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ABSTRACT 


The  current  concepts  of  the  tarnation  of  fall¬ 
out  from  land  surface  nuclear  detections  have  been 
reviewed.  Thermodynamic  equations  have  been  devel¬ 
oped  for  a  part  of  the  overall  cond-  r.satlo.n  process 
to  account  for  fractionation  of  the  radioactive 
species.  Empirical  functions  for  «v:«  of  the  fire¬ 
ball  parameters  have  been  developed  rrom  available 
data  and  assumptions  about  the  utix.zaticn  of  the 
energy  released  in  a  nuclear  explosj.-.r>  •  These 
functions  and  available  data  on  the  -apor  pressure 
of  fission  product  elements  and  compounds  (mainly 
oxides)  and  on  the  fission  yields  i>om  fission  of 
U235,  u23«,  and  Pu239  were  utilized  to  compute 
decay  curves  for  the  unfracti.-nated  mixtures  and 
for  an  idealized  fallout  condition  from  a  surface 
nuclear  detonation. 


SUMMARY 


The  Problem 


The  radiological  hazard  from  nuclear  detonations  arises  from 
the  fact  that  radioactive  materials  produced  in  the  explosion  become 
associated  with  rather  large  particles  of  e pvt rfr ruriental  materials 
that  fall  back  to  earth  shortly  afterward.  The  accumulation  of 
these  particles  on  the  earth's  surface  creates  a  gamma  radiation 
hazard  because  of  radioactive  content  of  the  particles. 

The  nature,  degree,  and  type  of  hazard  over  a  prolonged  time 
period  depends  upon  the  radioactive  composition  of  the  particles . 

Tt  is  known  that  the  decay  rates  of  the  fission  products  produced 
in  weapons  tests  differs  from  one  shot  to  another  and  also  differs 
from  one  location  to  another  for  a  given  shot.  It  is  also  known 
that,  in  a  nuclear  war  for  which  countermeasures  are  to  be  developed 
and  planned,  the  materials  at  likely  targets  hre  different  from 
those  at  weapons  test  sites .  Therefore  to  explain  the  causes  of 
the  different  observed  radiation  decay  of  fallout  and  better  define 
the  radiological  hazards  and  countermeasure  requirements,  a  better 
understanding  of  the  mechani3m(s)  of  fallout  formation  is  required. 

Findings 


In  this  report  an  attempt  is  made  to  describe  mathematically  some 
of  the  processes  of  fallout  formation.  A  general  process  of  fallout 
formation  is  outlined  based  on  the  observed  structures  and  types  of 
fallout  particles  produced  by  nuclear  explosions.  Thermodynamic 
equations  are  developed  to  describe  some  of  the  possible  types  of 
condensation  that  could  occur.  The  conditions  of  temperature,  pres¬ 
sure,  volume,  time,  and  energy  utilization  for  establishing  the 
boundary  conditions  of  the  condensation  processes  are  obtained  from 
some  of  the  data  in  The  Effects  of  Nuclear  Weapons!^-  together  with 
numerous  assumptions  in  order  to  make  estimates  for  land  surface 
detonations .  The  boundary  conditions  aid.  thermal  data  are  used  to 
estimate  the  amount  of  melted  soil  present,  i.n  the  fireball  when  the 
temperature  1ms  decreased  to  the  soil  melting  point. 

Data  on  the  vapor  pressures  of  the  fission  product  elements  and 
compounds  (mainly  oxides)  is  summarized.  The  fission  yields  for  slow 
and  fast  fission  of  u235_,  \j23 8,  and  Pu-39  are  discussed  and  summarized. 
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^or  some  of  the  fission  processes,  estimates  of  the  fission  yield  were 
made  to  complete  the  yield  curve  for  use  in  computations. 

Using  an  ideal  non-reactive  soil  with  a  defined  melting  point 
of  3.Jj00°C  and  which  formed  ideal  dilute  solutions  with  the  fission 
product  (usually  as  the  oxide)  elements,  computations  are  made  for 
a  surface  land  detonation  with  an  estimated  yield  of  2*3  MT  and  in 
which  the  temperature  of  lbOO°C  occurs  at  60  sec  after  detonation. 

For  this  case,  it  is  estimated  that  6.3  $  of  the  total  energy  was 
contained  in  the  liquid  particles  Just  prior  to  their  solidification 
at  the  melting  point • 

Calculations  are  made  of  the  disintegration  rates,  photon  emis¬ 
sion  rates,  photon-energy  emission  rates,  and  air-ionization  rates 
from  the  unfractionated  or  normal  mixture  of  figsion  products  from 
the  slow-  and  fast-neutron  fission  of  U235,  U23^  and  Pu239 .  The 
six  ionization  rate  decay  curves  gave  ratios  of  the  r/hr  at  1  hr 
(at  3  ft  above  an  infinite  smooth  plane)  to  KT  equivalents  per  sq 
mi  that  range  from  3^0  to  3950. 

For  the  idealized  detonation  conditions  in  which  the  computations 
assume  only  one  particle  size  (i.e.,  that  none  fall  out  of  the  fire¬ 
ball  up  to  60  sec,  and  then  all  the  particles  leave  in  solid  form), 
computations  of  the  disintegration  rates  and  air-ionization  rates  for 
thermal  and  fission-neutron  fission  of  U235,  8  Mev-neutron  (broad  hand) 
fission  of  U23”,  and  fission-neutron  fission  of  Pu239.  in  these  cases, 
the  ratios  of  the  r/hr  at  1  hr  (at  3  ft  above  an  infinite  smooth  plane) 
to  KT  equivalents  per  sq  mi  are  found  to  range  from  1^80  to  1560 ■  If 
induced  activities,  instrument  response,  and  terrain  roughness  factors 
are  considered,  this  ratio  could  be  as  low  as  1000.  However,  for  a 
fractionated  mixture,  the  ratio  has  no  real  meaning.  In  addition,  the 
ratio  is  net  constant,  except  for  a  single  computation;  in  a  read,  fall¬ 
out  area  it  will  change  with  particle  size  and  distance  from  shot  point. 
The  report  gives  some  suggestions  for  further  improvements  in  the  theory 
and  suggests  possible  methods  for  making  more  realistic  and  comprehen¬ 
sive  computations. 
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SECTION  1 


INTRODUCTION 


1.1  OBJECTIVES 


The  major  objectives  of  this  report  are  to  formalize  some  of 
the  current  concepts  of  formation  of  fallout  from  land-surface 
detonations  and  to  develop  approximating  formulae,  by  use  of 
thermodynamics  logic  and  empirical  functions,  for  the  purpose  of 
computing  -  or  estimating  -  the  radioactive  decay  of  such  fallout. 
Test  detonations  from  which  this  type  of  Information-  could  be  ob¬ 
tained  have  net  yet  been  conducted.  Another  objective  was  to 
develop  a  procedure  for  computing  decay  curves  and  to  demonstrate 
its  use  by  making  sample  calculations . 


1.2  SCOPE 


The  report  is  divided  into  seven  sections.  In  this  section 
( INTRODUCTION ) ,  the  general  background  material  is  summarized. 
Inferences  are  made  from  analyses  of  fallout  particles  with  re¬ 
gard  to  how  they  must  have  been  formed>-  and  some  definitions  of 
fractionation  are  given. 

In  Section  2,  the  condensation  process  in  the  formation  of 
fallout  is  described  and  equations  are  derived  for  several  alter¬ 
nate  processes. 

In  Section  3;  data  on  the  vapor  pressures  of  the  fission 
product  oxides,  elements,  and  some  compounds  together  with  the 
oxides  of  a  few  other  elements  are  summarized.  Some  equations 
are  given  to  illustrate  the  use  of  the  data  in  the  computations. 
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In  Section  4,  some  empirical  functions  are  derived  for  esti¬ 
mating  the  amount  of  energy  from  a  land  surface  nuclear  explosion 
that  is  available  for  vaporizing  and/or  melting  soil  particles. 

In  the  development  of  the  material  for  this  section,  correlation 
of  the  information  given  in  The  Effects  of  Nuclear  Weapons , which 
is  the  only  unclassified  source  available,  led  to  the  finding  of  a 
number  of  inconsistencies  in  data  and  data  treatment  in  that  docu¬ 
ment.  However,  it  is  emphasized  that  the  scope  of  this  report  does 
not  include  criticizing  that  document  but  rather  utilizing  its  con¬ 
tents  to  establish  the  approximating  formulae  necessary  for  this 
study . 


In  Section  5j  available  data  on  the  fission  yields  for  several 
fissile  materials  are  summarized.  Correlation  techniques  were  used 
to  estimate  the  yields  of  mass  chains  for  which  data  were  not  available. 

In  Sections  6  and  7,  the  computational  methods  and  sample  calcu¬ 
lations  of  the  decay  curves  are  described  for  both  unfractionated  and 
fractionated  fission  products .  The  calculations  are  made  for  diff¬ 
erent  fissile  materials.  The  comparisons  given  for  the  fractionated 
mixture  are  based  on  calculations  for  a  non-reactive  soil  that  formed 
an  ideal  solution  with  all  fission  product  oxides  in  the  liquid  state. 
Further  study  and  experimental  results  will  be  required  before  the 
calculations  can  be  made  for  a  process  that  more  closely  corresponds 
to  the  real  one. 


1.3  LIMITATIONS 


The  many  items  considered  in  the  general  treatment  all  have  a 
common  limitation:  the  lack  of  sufficient  available  unclassified 
data  for  testing  the  validity  of  individual  formula  or  the  completed 
computation. 

In  the  text  itself,  the  assumptions  and  postulates  are  stated 
usually  only  once  in  order  to  minimize  undue  interruption  of  the 
main  line  of  the  argument  and  development  of  the  material.  Conclu¬ 
sions  and  comparisons  are  often  stated  without  definite  repeated 
reference  to  the  original  assumptions.  This  limitation  is  given 
here  as  a  caution  to  the  reader  in  the  use  of  the  computational 
results  given  in  this  report. 

The  computational  results  were  made  for  illustrative  use  only 
since  an  idealized  chemical  system  was  used.  A  great  deal  of  addi¬ 
tional  work  can  be  done  to  improve  the  quantitative  aspects  of  the 
theory.  This  report  may  be  most  helpful  in  identifying  the  type  of 
information  needed  for  improving  both  the  conceptual  basis  for  deriv¬ 
ing  the  mathematics  and  the  quantitative  use  of  the  theory. 
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1.4  BACKGROUND 


In  a  nuclear  detonation  employing  fissionable  materials  about 
90  fission-product  mass  chains  (to  elements)  are  produced.  The 
high  temperatures  immediately  following  the  detonation  virtually 
assure  that  all  of  these  elements  are  immediately  vaporized  so 
that  the  fallout  formation  process  is  one  in  which  these  elements 
condense  from  a  vapor  phase  during  the  cooling  of  the  fireball. 
Among  the  90-odd  mass-chains  produced,  with  specific  fission  yields 
depending  on  the  fissile  material  and  the  neutron  energy,  there  are 
many  atomic  species  and  compounds  whose  volatilities  at  high  tem¬ 
peratures  are  very  different.  Hence,  at  a  given  temperature,  the 
range  of  values  of  the  equilibrium  partial  vapor  pressure  of  the 
various  species  and  compounds  will  be  large.  A  condensation  pro¬ 
cess  that  occurs  in  a  system  of  a  rapidly  decreasing  high  temper¬ 
ature  can  be  expected  therefore  to  result  in  the  preferential 
condensation  of  the  less  volatile  elements.  This  should,  in  turn, 
result  in  an  alteration  of  the  relative  abundance  of  the  fission 
product  mass  chains  or  radionuclides  as  found  in  fallout  compared 
to  the  amounts  initially  produced.  Any  such  alteration,  when 
observed,  is  usually  called  fractionation. 


A  general  description  of  3ome  of  the  fallout  particles  found 
in  local  fallout  from  surface  and  tower  shots  and  of  how  they  could 
have  originated  has  been  given  by  Adams  and  coworkers. 3*4,5, o,7 


The  essential  details  of  this  work  that  have  been  used  in  this  report 
are:  (l)  condensation  of  the  fission  product  elements  can  begin  at 
the  highest  temperature  at  which  a  macroscopic  liquid  phase  can 
exist  in  the  fireball  and  this  liquid  phase  will  essentially  consist 
of  substances  such  as  iron  oxide  (tower  shots)  or  aluminum  and  sili¬ 
con  oxides  (surface  ground  shot),  and  (2)  the  fission  products  that 
co -condense  with  or  into  liquid  particles  are  dissolved -into  the 
melts  and  remain  there  as  compounds  or  in  solid  solution  when  the 
particle's  solidify. 


In  low  tower  shots,  it  appears  that  the  very  small  drops  of 
vapor -condensed  iron  oxide  are  later  dissolved  by  melted  soil  par¬ 
ticles  to  form  a  glassy  iron  oxide -silicon  oxide  solution.  ^Data 
on  the  solubility  of  fission  product  elements  into  dilute  acids  and 
complexing  agents,  obtained  by  Fuller^  from  the  larger  particles 
from  a  low  tower  shot  that  fell  close  to  ground  zero,  show  that 
only  a  very  small  fraction  of  the  fission  products  are  soluble. 
However,  data  on  the  solubility,  obtained  by  Larson^  from  smaller 
particles  collected  at  greater  distances  from  shot  point,  give 
larger  fractions  of  soluble  fission  products;  his  data  show  also 
that  radioelements  other  than  just  the  daughter  products  of  rare 
gas  nuclides  are  among  the  soluble  group  of  elements. 
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The  essential  facts  from  these  data  are  that  some  of  each  of 
the  fission  products  condenses  into  liquid  particles  and  that  some 
of  each  condensed  onto  the  surfaces  of  solidified  particles.  Also, 
the  particles  that  fall  in  the  local  areas  downwind  from  shot  point 
(and  where  the  amount  that  falls  out  is  large  enough  to  produce  a 
significant  radiation  hazard)  are  too  large  to  have  teen  produced 
"by  a  vapor  condensation  process.  These  particles  therefore  are 
formed  either  from  the  break-up  of  a  bulk  liquid  melt  or  from  the 
melting  of  single  grains  of  soil  that  enter  the  fireball  after  it 
has  cooled  to  some  given  temperature . 

The  general  condensation  process  can  therefore  be  divided  into 
two  time  periods.  The  first  period  of  the  process  is  characterized 
by  the  presence  of  gas  and  liquid  phases  and  the  second  period  by 
the  existence  of  gas  and  solid  phases.  The  first  period  of  conden¬ 
sation  ends  when  the  bulk  carrier  or  substrate  material  of  the  par¬ 
ticles  solidifies.  The  degree  to  which  the  change  of  phase  of  the 
carrier  affects  the  continuity  of  the  process  and  the  distribution 
of  each  fission  product  element  in  or  on  a  particle  will  depend  on 
the  thermo -chemical  properties  of  the  interacting  materials.  Some 
fission  product  elements  may  condense  by  sublimation  on  the  surface 
of  the  particle  and  be  readily  available  for  solution  upon  contact 
with  water.  Others  may  react  with  the  substrate  material  and/or 
diffuse  through  the  surface  of  the  particle. 

One  particular  aspect  of  first  period  of  condensation  is  im¬ 
portant.  This  is  that  the  fission  product  vapors  condense  into 
the  liquid  phase  of .  the  carrier  material  to  form  a  very  dilute 
solution. 

The  fraction  of  each  fission  product  element'  that  condenses 
into  the  liquid  carrier  particles  depends  on  the  melting  point  of 
the  carrier  and  the  time  after  fission  at  which  that  temperature 
occurs.  If  the  melting  point  of  the  carrier  is  high,  the  fractions 
condensed  will  be  low.  If  the  melting  point  of  the  carrier  is  low, 
the  fractions  condensed  will  be  high.  The  fraction  which  had  not 
condensed  into  the  liquid  phase  of  the  carrier  can  condense  on  or 
react  with  the  surface  of  the  solid  particles.  These  could  consist 
of  the  smaller  of  the  melted  particles  or  of  unmelted  particles  that 
enter  the  gas  volume  at  later  times. 


1.5  FACTORS  CONTROLLING  FRACTIONATION  IN  FALLOUT 


In  an  overall  sense,  the  amount  of  each  fission  product  found 
in  fallout  relative  to  some  standard  of  comparison  depends  on  five 
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main  factors:  (l)  the  original  fission  yields  or  relative  abundance, 
of  the  fission  products,  (2)  neutron  capture  by  the  fission  products 
themselves,  (3)  the  degree  to  which  each  fission  product  condenses 
into  or  onto  the  carrier  particles,  (b)  neutron  emissions  in  the 
decay  chain,  and  (5)  the  radiochemical  standards  used  to  measure  the 
fractionation . 

The  relative  abundance  of  each  fission  product  element  origin¬ 
ally  produced  depends  on  the  fissile  material  used  —  i.e.  whether 
the  material  is  11^35,  u230j  Pu239,  or  some  other  element.  The  fission 
chain  yields  also  depend  on  the  energy  spectrum  of  the  incident  neu¬ 
trons.  In  comparison  with  the  fission  yields  from  thermal  neutrons 
on  U235,  usually  taken  as  the  reference  standard,  the  yields  of  the 
elements  in  the  lighter-mass  peak  for  other  types  of  fission  in  gen¬ 
eral,  shift  more  than  do  those  in  the  heavier  mass  peak.  With  the 
heavier  fissile  elements,  the  center  of  the  lighter-mass  group  moves 
toward  the  higher  mass  numbers.  As  the  incident  neutron  energy  in¬ 
creases,  the  yields  of  the  valley  elements  and  those  of  highest  and 
lowest  mass  numbers  rise  and  the  neutron  yield  per  fission  Increases. 
The  increase  in  neutron  yield  per  fission  tends  to  spread  the  two 
peaks  farther  apart  and,  again,  the  lighter-mass  group  is  shifted 
more  than  the  heavy  group. 

Neutron  capture  by  the  fission  product  elements  would  result  in 
a  general  shift  of  the  whole  yield  curve  to  higher  mass  numbers.  The 
result  would  be  a  decrease  in  the  yields  of  the  elements  with  the 
smaller  mass  numbers  (left  side  of  the  peaks)  of  both  groups,  and  an 
increase  in  the  yields  of  the  elements  with  the  larger  mass  numbers 
of  both  groups.  Relatively  little  change  would  result  in  the  yields 
of  elements  in  the  peaks  excepting  for  those  that  may  have  extremely 
high  capture  cross-sections.  The  subject  is  not  discussed  further  in 
this  report  because  of  insufficient  data. 

Section  2  of  this  report  discusses  the  role  of  the  condensation 
process  in  fallout  formation  and  in  fractionation.  This  process  is 
often  assumed  (and  erroneously)  to  be  the  only  cause  of  observed  frac¬ 
tionation  in  fallout. 

Neutron  emission  during  the  decay  process  results  in  a  product 
nuclide  with  a  mass  number  one  unit  less  than  the  parent.  This  chain 
"shift"  can  be  accounted  for  if  the  decay  scheme  is  known.  However 
for  many  of  the  short-lived  radionuclides  there  is  insufficient  data 
for  its  further  consideration  in  this  report. 

The  experimental  measure  of  fractionation  is  most  often  given 
as  an  "R"  factor  or  value  relative  to  thermal  fission  of  U235  and 
a  selected  radionuclide.  The  most  commonly  selected  radionuclide 
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for  comparison  is  Mo99.  Thus,,  relative  to  these  standards,  a  radio¬ 
chemical  assay  of  a  fallout  sample  that  gives  an  "R"  value  different 
from  one  does  not  necessarily  mean  that  the  nuclide  in  question  has, 
in  fact,  teen  fractionated.  Knowledge  of  the  true  initial  fission 
yields  are  required  to  correct  the  observed  assay  data  to  determine 
whether  fractionation  has,  in  fact,  occurred. 
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SECTION  2 


THE  CONDENSATION  PROCESS 


2.1  THE  FIRST  IERIOD  OF  CONDENSATION 


2.1.1  Description  of  the  Formation  Process 


As  is  implied  by  the  general  description  of  the  formation 
process  given  in  Section  1,  no  single  well-defined  condensation 
time  for  the  fission  product  mixture  can  exist.  Rather*  it  is  a 
continuing  process  that  can  be  separated  into  the  two  more  or  less 
well-defined  time  periods  described  above.  In  the  first  time  period 
(during  the  cooling  of  the  fireball),  the  major  feature  of  the  con¬ 
densation  process  is  the  existence  of  vapor-liquid  phase  equilibria. 
This  period  of  condensation  ends  when  the  carrier  material  solidi¬ 
fies,  with  the  fission  products  being  either  fixed  in  a  solid  solu¬ 
tion  matrix  or  compounded  with  the  carrier  material. 

The  major  feature  of  the  second  period  of  condensation  is  the 
existence  of  vapor-solid  phase  equilibria  in  which  the  remaining 
fission  product  elements  condense  at  lower  temperatures  on  the  sur¬ 
face  of  solid  particles.  The  second  period  of  condensation  never 
ceases  in  an  absolute  sense,  except  for  those  particles  which'  fall 
out  or  otherwise  leave  the  space  containing  the  residual  gases. 
Actually,  the  process  can  reverse  for  a  fission  product  element 
that  later  decays  to  a  more  volatile  element;  for  example,  elements 
like  iodine  and  the  rare  gases  could  sublime  as  fast  as  they  form 
from  non-volatile  precursors  (i.e.,  at  ordinary  temperatures)  con¬ 
densed  on  the  surface  of  fallout  particles.  This  process  is  unlikely 
when  the  fission  products  are  trapped  within  a  glassy  matrix,  since 
the  vapor  pressures  due  to  the  lev  concentration  of  dissolved  fission 
products  (ca.  10 "10  moles/mole)  would  be  extremely  low,  and  diffusion 
through  the  glass  solid  would  be  very  slow. 


7 


The  essential  problem  in  the  theory  for  the  process  during 
the  first  period  of  condensation  is  to  establish  the  vapor-liquid 
phase  equilibria  of  each  fission  product  element  at  the  time  that 
the  carrier  material  solidifies;  that  is,  to  determine  the  fraction 
of  each  element  present,  condensed,  Inside  the  carrier  melt  when  it 
solidifies. 


When  one  of  two  phases  in  contact  is  a  gas,  simple  kinetic 
theory  can  be  used  to  show  that  condensation-vaporization  equili¬ 
brium  can  be  established  within  a  small  fraction  of  a  second  at 
temperatures  above  1000°K.  Thus  those  gaseous  species  of  each 
fission  product  element„that  do  not  react  with  the  liquid  carrier 
but  dissolve  into  it  should  obey  Henry's  law  of  dilute  solutions. 

In  fact,  the  solutions  should  be  sufficiently  dilute  as  to  result 
in  no  change  in  the  free  energy  of  the  liquid  carrier  so  that  the 
free  energy  of  each  element  in  the  solution  should  be  independent 
of  any  other.  Therefore,  it  is  possible  to  consider  the  solubil¬ 
ity  of  each  element  as  making  a  two-component  system  with  the 
carrier.  Moreover,  there  should  be  no  appreciable  surface  loading 
(large  excess  surface  concentrations)  during  the  condensation  pro¬ 
cess  if  the  temperature  range  over  which  the  liquid  carrier  exists 
exceeds  200  or  300°C -  A  concentration  gradient  in  a  particle,  how¬ 
ever  should  exist  especially  for  the  larger  particles  of  which  some 
may  not  be  melted  in  the  center  when  the  air  or  gas  temperature 
about  the  particle  falls  below  the  melting  point  of  the  bulk  carrier. 


2.1.2  Condensation  Thermodynamics 


For  the  condensation  of  dnj  moles  of  element  j  from  a  gaseous 
mixture  tg  a  liquid  solution,  the  change  in  the  partial  molar  free 
energy,  Fj,  of  the  element  in  the  gas  is  given  by 
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and  the  corresponding  change  in  its  partial  molar  free  energy,  Fj,  in 
the  solution  is 
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where  P  is  the  total  pressure,  T  is  the  temperature,  N.  is  the  mole 
fraction  of  element  j  in  the  gas  mixture,  Nj  is  its  mole  fraction 
in  the  liquid,  f?  is  its  fugacity  (idealized  pressure)  in  the  gas 
phase  and  fj  is  its  fugacity  in  the  liquid  phase.  At  moderate  and 
low  total  pressures,  the  fugacity  of  the  element  in  the  gas  phase 
is  given  by 
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in  which  kj  is  the  Henry's  law  constant  at  a  given  temperature  and 
total  pressure;  hence 


In  an  equilibrium  process,  the  two  changes  in  the  partial  molar 
free  energies  for  the  transfer  of  dn<  moles  from  the  gas  mixture  to 
the  liquid  solution  are  equal.  Equating  Eqs.  1  and  2  and  substituting 
the  appropriate  thermodynamic  equivalents  for  the  indicated  partial 
differentials  gives 

-O  -O  Vr  — 

vi  Li  .  n  vi 

dP sL,  dT  +  d£nN^  =  -1  dP  -  -X  dT  +  d/nl^  (7) 

RT  RT2  0  RT  rt2  J  k  ‘ 

in  which  Vj  is  the  partial  molar  volume  of  element  j  in  the  gas  mix¬ 
ture,  Lj  is  its  relative  partial  molar  heat  content  in  the  gas  mix¬ 
ture,  Vj  is  its  partial  molar  volume  in  the  liquid  solution,  Lj  is 
its  relative  partial  molar  heat  content  in  the  liquid  solution,  R  is 
the  molar  gas  constant,  P  is  total  pressure,  and  T  is  the  temperature 
in  °K.  For  dilute  solutions. 
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vhere^Hy  is  the  heat  of  vaporization  of  the  condensing  specie  of 
element  j .  Actually,  AHV  is  the  heat  of  the  reaction  for  the  vapor¬ 
ization  of  the  gas  from  its  form  in  solution.  If  it  exists  as  a  dif¬ 
ferent  compound  in  solution,  the  heat  of  formation  of  this  compound  is 
included  in  the  value  of  AHV.  For  an  ideal  solution,  or  one  in  which 
there  is  no  heat  of  dilution  or  compound  formation,  4HV  is  the  actual 
heat  of  vaporization.  Since  Vj>=*Vj,  (Vj/RT)cLP  (or,dinf j/6p)  can  be 
neglected.  For  an  ideal  gas,  r°  can  be  replaced  by  N°P.  Then 

l  6  p 

Substituting  Eqs.  8  and  9  into  Eq.  7  and  integrating  gives 
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in  which  k.  is  an  integration  constant  and  where  the  term  k.e 
is  identifiable  via  Eqs.  3  and  5  (f9  »  =  p<,  the  partial^pressure 

of  element  j;  and  fj  =  P)  as  the  Henry 'sdlaw  constant.  It  may  be 
noted  from  Eq.  10  that  an  increase  in  P  results  in  a  decrease  in  the 
ratio  or  an  increase  in  the  mole  fraction  of  the  minor  constitu¬ 

ent  j,  in  the  liquid  phase  relative  to  its  mole  fraction  in  the  gas 
phase.  For  early  time  condensations  (high  M.P.  of  carrier),  when  the 
temperature  is  high  and  the  gas  volume  not  fully  expanded,  the  total 
pressure  should  be  high;  thus  the  effect  of'  the  pressure  alone  would 
tend  to  produce  a  more'  complete  condensation  at  shorter  times  after 
detonation.  However,  the  Henry's  law  constant  also  is  larger  at 
higher  temperatures  so  that  the  two  terms  in  Eq.  10  would  decrease 
(or  increase)  simultaneously  and  the  change  in  the  ratio  of  the  two 
mole  fractions  with  time  (or  temperature)  in  the  expanding  gas  volume 
will  be  less  than  for  just  the  decrease  in  temperature. 


2.1.3  Effect  of  Carrier  Particle  Size  on  Mole  Fraction  Ratio 


A  dependence  of  the  ratio  of  the  mole  fractions  on  the  size 
of  the  liquid  drop,  itself  may  exist  because  the  dependence  of  the 
total  pressure  about  each  drop  on  drop-size. 
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This  dependence  is  given  by 

RTinp/p0  =  hybl/fid. 

in  which  p0  is  the  vapor  pressure  of  the  carrier  material  over  the 
liquid  with  a  flat  surface.,  p  is  the  pressure  over  the  drop  of  dia¬ 
meter,  d,  y  is  the  surface  tension  of  the  drop  (assumed  independent 
of  Tj,  M  is  the  molecular  weight  of  the  carrier,  and Pis  the  density 
of  the  liquid.  For  AI2O3  at  2050°C  and  SiOg  at  1800°C,  the  value  of 
7  is  690  dynes/cm  and  30 1  dynes/cm,  respectively. The  value  of 
the  ratio,  p/p0,  for  these  values  of  7  is  not  very  different  from 
unity  for  diameters  larger  than  a  few  tenths  of  a  micron.  Hence, 
unless  the  carrier  material  has  a  surface  tension  larger  than  do 
these  two  oxides  by  more  than  two  orders  of  magnitude,  the  increased 
vapor  pressure  of  the  carrier  material  over  the  larger  drops  should 
not  be  enough  to  result  in  larger  values  of  Nj.  If  the  surface  ten¬ 
sion  were  extremely  large,  the  fission  product  elements  most  likely 
to  be  preferentially  condensed  on  the  smaller  particles  are  those 
whose  volatilities  are  the  same  or  lower  than  that  of  the  carrier 
itself.  These  would  co-condense  with  the  volatilized  carrier 
molecules  as  soon  as  the  temperature  dropped  to  the  carrier  boil¬ 
ing  point,  since  at  this  time  the  vapor  pressure  of  the  carrier 
material  would  be  a  large  fraction  of  the  total  pressure.  At  the 
melting  point  of  the  carrier,  its  own  vapor  pressure  would  be  so 
small  a  fraction  of  the  total  pressure  that  it  could  hot  influence 
the  mole  fraction  ratio  even  if  the  surface  tension  were  extremely 
large . 

For  particles  with  a  fairly  large  range  in  sizes,  the  mole 
fraction,  Nj,  needs  to  be  precisely  defined.  Carrier  materials  such 
as  silicate  soils  containing  metal  oxides  are  refractory  materials 
with  low  heat  conductance.  As  mentioned  above,  a  particle  need  not 
be  completely  melted  throughout  its  volume  in  order  to  condense 
gaseous  molecules.  Only  a  liquid  layer  of  film  on  its  exterior  is 
required  for  the  process.  Also,  since  the  time  spent  in  the  liquid 
state  is  short,  the  condensates  may  not  diffuse  more  than  a  short 
distance  from  the  surface  of  the  particle  before  it  solidifies. 

Thus  for  particles  of  a  given  temperature  history,  a  maximum  size 
particle  should  exist  that  is  completely  melted  before  it  resolidi¬ 
fies.  Elements  or  nuclides  that  condense  first  (less  volatile  than 
the  carrier)  should  penetrate  somewhat  farther  into  the  completely 
melted  particles  than  those  that  are  more  volatile  and  condense 
more  rapidly  just  before  the  particle  solidifies  if  the  penetration 
depended  only  on  diffusion. 

The  penetration  of  the  condensates  into  the  liquid  drop  and 
its  rate  of  distribution  throughout  the  volume  would  be  more  rapid 
if  caused  by  turbulence  and  convection  -  especially  in  the  peripheral 
regions  of  the  liquid  particles.  The  general  uniformity  of  the 
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radioactive  concentration  in  the  silicate  fallout  particles  indicate 
a  process  of  distribution  other  than  by  diffusion.  The  resultant 
process  appears  to  be  one  in  which  the  condensates  are  deposited  to 
a  given  depth  from  the  surface  of  the  very  largest  particles  and  in 
a  more  or  less  uniform  concentration  throughout  the  medium-size  and 
smaller  particles.  The  more  volatile  elements  may  form  solutions 
with  the  larger  surface  concentration  excesses  if  they  are  conden¬ 
sing  in  large  amounts  just  as  the  particles  are  solidifying. 

If  the  average  depth  of  the  surface  layer  is  designated  as 
h  and  is  assumed  to  be  the  same  for  all  particles,  then  the  number 
of  moles  of  carrier  per  particle  (spherical)  that  was  involved  in 
forraing  the  dilute  solution  (i.e.  was  melted  during  the  formation 
process)  is  given  by 


n(£,p)  =  UP*  [d(d-2h)  +  (4/3)h2]  ,  d^2h  (12) 

in  which  n(£ , p)  is  the  number  of  moles  of  she  carrier  (that  was  melted) 
in  the  surface  layer  of  the  particle. 

For  dilute  -solutions,  the  mole  fraction  is 

Nj  =  n  (p)/n(£,p)  (13) 

in  which  n^(p)  is  the  number  of  moles  of  element  j  in  the  particle; 
the  latter0 is  then  given  by 

nj(p)  =  Nj  ^  [d(d-2h)  +  (4/3)h2]  ,  d^>2h  (l4) 

If  n9(p)  is  taken  as  the  total  number  of  moles  of  the  carrier  in  the 
particle,  then 

n^(p)/nQ(p)  =  6.00  Njh  fe..d.~2h)^..  ( V3)hg]  ^  d>.  2h  (15) 

nj(p)/n0(p)  =  N^,  d^2h  (16) 

The  mole  fraction  in  Eq.  16  is  for  the  particles  that  were  completely 
melted.  The  ratio,  n,(p)/n (p)N,,  for  Eqs.  15  and  16  is  shown  in 
Fig.  1  for  several  values  of  h.  if  the  value  of  Nj  for  each  fission 
product  is  the  same  for  all  sizes  via  Eq.  10  and  discussion  of  Eq.  11, 
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then  the  gross  average  concentration  of  fission  products  per  unit 
mass  of  particles  should  he  a  multiple  of  a  curve  similar  to  those 
in  Fig.  1  for  some  value  of  h.  The  value  of  h  could  be  determined 
from  an  experimentally  determined  curve  by  extending  the  straight 
line  portions  of  the  curve  with  logarithmic  slopes  of  0  and  -1  to 
their  intercept  at  a  size,  which  is  6  times  the  value  of  h. 

The  case  for  a  non-uniform  concentration  in  the  surface  layer 
or  throughout  the  particle  was  not  treated  for  presentation  in  this 
report.  Its  consideration  requires  estimates  of  the  gaseous  diffu¬ 
sion  rates  for  condensation  at  the  particle  surface  as  well  as  esti¬ 
mates  of  subsequent  mixing  rates  of  the  condensates  in  the  liquid 
layer  of  the  particle.  Such  estimates  may  he  required  for  a  more 
complete  treatment  of  the  condensation  process. 

There  are  two  additional  factors  that  could  have  some  bearing 
on  the  relative  amount  of  an  element  condensing  into  the  liquid  par¬ 
ticles.  The  larger  particles  may  either  fall  out  of  the  gaseous 
volume  before  they  solidify  or  may  be  melted  only  on  one  side  (not 
sufficiently  melted  to  form  a  spherical  particle) .  Both  of  these 
two  groups  of  particles  would  condense  smaller  amounts  of  all  ele¬ 
ments  per  particle  than  would  be  estimated  from  Eq.  15  and  also 
would  be  most  deficient  in  the  more  volatile  elements. 

The  origins  of  the  melted  particles  in  the  fireball  may  be 
several  in  number.  The  small  vapor-condensed  particles  originating 
from  vaporized  soil  have  been  mentioned.  Some  of  the  fallout  par¬ 
ticles  are  particle's  that  were  originally  lying  on  the  ground  out 
to  some  distance  from  shot  point,.  Some  of  these  probably  are  melted 
by  the  heat  absorbed  from  the  radiant  energy  emitted  at  detonation. 
These  particles  are  then  drawn  into  the  fireball  as  it  rises.  This 
mechanism  probably  is  the  dominant  one  for  low  air  bursts  where  a 
small  crater  is  formed.  In  this  case  the  size  distribution  of  the 
fallout  particles  that  are  produced  should  be  the  same  as  the  original 
size  distribution  of  the  surface  soil.  In  the  detonation;,  the  blast 
wave  would  powder  the  soil  to  some  depth  and  the  resulting  dust  par¬ 
ticles,  like  the  surface -melted  particles,  would  retain  their  sizes 
in  the  formation  process.  However,  these  particles  would  not  be 
melted  until  after  they  enter  the  fireball. 

For  surface  detonations ,  a  layer  of  soil  may  be  melted  in  the 
process  of  crater  formation.  It  would  seem  reasonable  that  a  liquid 
layer  of  soil  would  intervene  between  the  vapor  in  the  fireball  and 
a  layer  of  shock-powdered  earth  as  long  as  the  fireball  remained  in 
contact  with  the  earth's  surface.  If  enough  "fluxing"  material  such 
as  the  carbonates  are  present,  a  rather  large  amount  of  fluid  material 
could  be  formed.  As  the  fireball  lifts,  the  liquid  would  break  up 
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into  drops  and  enter  the  fireball,  followed  by  the  powdered  soil. 
The  mechanical  break-up  of  the  fluid  mass  should  produce  about 
the  same  size  distribution  of  particles  from  liquids  that  have 
about  the  same  surface  tension.  If  most  of  the  fallout  particles 
originate  from. the  bulk  liquid  layer,  then  the  particle  size  dis¬ 
tributions  in  the  fireball  from  detonations  on  different  soils 
would  be  more  nearly  alike  than  if  most  of  them  originate  from 
the  surface  soils  and  the  shock -powdered  dust  from  the  crater. 

There  is  no  unique  method,  then,  of  selecting  a  size  distri¬ 
bution  for  the  liquid  particles  present  in  the  fireball  up  to  the 
time  they  solidify.  But  since  Eqs.  12  through  16  will  not  lead  to 
a  variation  in  the  relative  amount  condensed  in  the  particles  for 
different  elements  because  of  size,  it  is  convenient  to  consider 
all  the  melted  particles  as  one  bulk  liquid  phase  in-  a  volume  of 
space  and  designate  the  distribution  in  general  terms  only  as 


i=m 

n(jt)  =  ^~>J1n1(l,p) 

i=l 


(IT) 


where  n(£)  is  the  total  moles  of  liquid  carrier  in  the  volume  V  at 
any  time,  ^  is  the  number  of  particles  of  size  i,  and  n^(/,p)  is 
the  same  as  n(j0,p)  for  the  size  i.  With  Eq.  IT,  the  mole  fraction 
of  element  j  in  the  liquid  phase  is 

Nj  =  nj/n(i)  (18) 

where  n^  is -the  number  of  moles  of  element  j  dissolved  in  the  n(/) 
moles  of  liquid  carrier.  The  mole  fraction  of  element  j  in  the 
vapor  phase  is  given  by 


=  n9/n 


(19) 


where  is  the  number  of  moles  of  element  j  mixed  with  n  moles  of 
vapor.  If  the  perfect  gas  law  is  used  for  n,  Eq.  19  becomes 


o 

N,  =  _ 

J  PV 


(20) 


Substitution  of  Eqs.  18  and  20  into  Eq.  10  gives 


M°  _  njkj 
j  ”  ( n(/)/V)RT 


(21) 
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where  k<  is  the  Henry’s  law  constant.  It  may  he  noted  that  the  ratio, 
n9/nj, depends  on  kj,  T,  and  the  number  of  moles  of  liquid  carrier  per 
unit  volume  in  the  fireball.  Since  the  mole  fractions  are  small,  the 
nj/nj  ratio  is  independent  of  the  total  amount  of  the  element  present. 
Thus J the  same  fraction  of  the  element  is  condensed  for  a  100  $  fission 
weapon  as  for  a  "clean"  thermonuclear  weapon  of  the  same  total  yield. 
The  value  of  n(£)/v  depends  only  on  the  total  yield  (see  Section  4) . 

2.1.4  The  Material  Balance  and  "R"  Factor  Equations 


Since  the  elements  considered  are  radioactive,  the  number  of 
moles  of  each  is  constantly  changing;  therefore  the  material  balance 
of  a  given  element  also  changes  with  time.  On  the  other  hand,  the 
fission  yield  of  a  given  mass  chain  (except  for  neutron  emitters)  is 
constant.  If  the  amount  of  a  radionuclide  of  element  j  and  mass 
number  A  present  at  the  time,  t,  after  fission  is  given  by  y^(t), 
then  the  total  amount  of  element  j  present  at  time  t  is  J  * 


Y,(t)  y,*(t)  atoms  or  moles 

J  A 


(22) 


The  corresponding  sum  for  the  chain  yield  of  mass  number  A  is 

Ya  =1,  yjA(t)  atoms  or  moles  (23) 

J 

in  which  YA  is  constant  except  for  the  mass  chains  containing  neutron 
emitters.  The  material  balance  for  element  j  between  gas  and  liquid 
phases  is 

Yj(t)  =  n°(t)  +  n^t)  (24) 


where  the  time  dependence  of  n.  and  ni  of  Eq.  21  is  indicated.  If 
and  n  are  taken  for  the  amounts  of  each  nuclide  in  the  gas  and 
liquid  phases,  respectively,  the  material  balance  for  each  element  is 


£a  VtJ  =  f A  njA(t)  +^A  V*1  (25) 

If  k?  is  used  for  kj/[”(n(£)/v)RT],  substitution  of  Eq.  21  in  Eq.  25 
(Eq.  21  holds  for  all  nuclides  of  element  j)  gives 


1 6 


(26) 


Ea  yjA(t>  =  (1+kj)  £a  njA^t} 


and  since  (l+k°)  is  the  multiplier  for  each  of  the  n.^(t)  terras,  the 
terms  for  eacn  nuclide  can  be  separated  out.  The  separated  terms 
give  the  amount  of  each  nuclide  condensed,  which  is 


“JA( 


t)  =  yjA(k) 


l+k9 


(27) 


The  fraction  of  the  mass  chain  condensed,  r  (A,t),  is  defined  b; 


(2b) 


Substitution  of  Eq.  27  and  replacing  the  y..(t)/Y.  ratios  by  y.:(A,t) 
in  Eq.  28  gives  J 


r0(A,t)  =y 

^ j  l+k? 

J 


(29) 


in  which  yj(A,t)  is  the  fraction  of  the  chain  yield  of  element  j 
with  mass  number  A.  According  to  Eq.  29,  the  fraction  of  the  chain 
condensed  in  the  liquid  phase  depends  only  on  kj  and  the  fractions 
of  the  chain  yield  of  the  elements  present.  The  use  of  Eq.  29 
therefore  requires  values  of  the  independent  yields  of  each  member 
element  of  the  chain  at  times  from  a  few  seconds  after  fission,  or, 
at  least  for  the  time  of  the  end  of  the  first  period  of  condensation. 

The  experimental  radiochemical  "R"  factor  for  a  given  sample 
or  quantity  of  mixed  fission  products  is  usually  defined  as  the  ratio 
of  the  number  of  atoms  of  a  given  mass  present  to  the  number  of  atoms 
of  mass  number  99  present,  divided  by  the  expected  value  of  their 
ratio  for  thermal  neutron  fission  of  U235.  in  mathematical  notation, 
this  is 

n(A) 

K0(A)  n(99)  (30> 
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Ordinarily*  count-rate  ratios  are  used  along  with  the  appropriate 
decay  corrections  from  the  time  of  analysis  to  zero  time  with  a 
corresponding  value  of  Kq(A)  that  has  "been  previously  determined 
from  an  analysis  of  a  sample  of  U^35  bombarded  with  thermal  neu¬ 
trons.  For  analyses  made  more  than  several  days  after  fission, 
only  the  last  or  last  two  members  of  a  decay  chain  of  most  mass 
numbers  will  be  present  in  appreciable  amount.  Since,  as  mentioned 
in  the  introduction,  the  observed  "R"  factors  give  an  overall  meas¬ 
ure  of  variation  from  U^35  thermal  fission,  another  factor  is  needed 
to  account  only  for  that  part  due  to  the  first  period  of  the  conden¬ 
sation  process.  For  this,  let 


r.(A) 


n(A) 

K(a)  n(99) 


(31) 


in  which  K(A)  is  the  true  yield  ratio,  of  mass  number  A  to 

99,  which  varies  with  the  kind  of  fissioning  nuclide  and  incident 
neutron  spectrum  (the  same  definition  holds  if  some  mass  number 
other  than  99  is  selected  as  the  reference  nuclide).  The  value  of 
the  ratio,  n(A)/n(99).>  should  be  the  same  as  the  ratio  of  the  sum 
over  j  of  the  respective  njA(t)  terms.  Making  the  substitution  in 
Eq.  31  gives 


r(A) 


^~‘jn  ja(  i) _ 

h(A)EjnJ99(t) 


(32) 


or,  with  Eq.  27, 


■(A)  = 


£,  qA(t)/(i+k°) 


and  finally,  replacing  K(A)  with  Y^/Y^,  gives 


E.  yj(A,t)/(l+k9) 

Ej  y-j(99,t)/(l+h°) 


(33) 


(3*0 
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Thus  r(A)  is  equal  to  rQ(A, t)/rQ(99j t)  of  the  tneory  and  equal  to 
Kq(A)r99(a)/k(A)  for  experimentally  determined  "R"  factors. 

2.1.5  Condensation  by  Compound  Formation  With  Carrier  Material 


Compound  formation  of  an  element  with  the  carrier  in  the  gas 
phase,  followed  hy  condensation  of  the  heavier  gas  molecule  into  a 
liquid  solution  with  the  melted  carrier  material,  would  be  described 
by  use  of  Henry's  law  for  the  dilute  solution.  However  when  the  com¬ 
pound  is  formed  with  the  (bulk)  liquid  carrier,  then  the  free  energy 
of  formation  of  the  compound  and  its  vaporization  and  that  of  the 
carrier  is  involved.  The  overall  reaction  for  a  direct  combination 
of  fission  product  element  A  with  carrier  material  B  to  form  the 
compound  AB  may  be  written  as  follows: 


A(g)  4-  B(X) -» AB(£),  dilute  solution  in  B(£.)  (35) 


This  reaction  may  be  written  as  a  sum  of  three  or  four  separ¬ 
ate  reactions  depending  on  whether  the  combination  of  A  and  B  occurs 
as  A  condenses  or  with  B  in  the  vapor  phase  prior  to  condensation. 
There  should  be  no  difference  in  the  two  processes  with  regard  to 
the  total  change  in  free  energy  between,  the  same  initial  and  final 
states.  The  difference,  if  any,  would  he  in  the  kinetics  of  the 
process.  In  the  first  process,  the  separate  reactions  and  the  stand¬ 


ard  free  energy  changes  are: 

1.  A(g)~»A  (l)  (36a) 

A =  KT  /n  pA  (36b) 

2.  A(£)  +  B(£)-*AB(£)  (37a) 

AF2=-RT^nKAB  (37b) 

where  K^g  is  the  equilibrium  constant  for  the  formation  of  AB(£)  from 

A(jL)  and  B(£.)  at  the  temperature,  T. 

3.  AB(£) — ?AB(£),  dilute  solution  in  B (£)  (3&a) 

A*3  “  -RT  l n  (38b) 
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where  a^  is  the  thermodynamic  activity  of  AB(Z)  in  b(£)  and  is 
equal  to  NAB  k^B  if  %£  is  the  mole  fraction  and  is  the  Henry's 
law  constant.  The  sum  of  the  standard  free  energy  changes  for  the 
three  reactions,  or,  the  free  energy  of  the  reaction  given  by  Eq. 
35,  is 


A?°  =  RT/n 


NAB  kAB  KAB 


(39) 


In  the  second  process,  the  separate  reactions  and  the  standard 
free  energy  changes  are: 


1. 

A(g)  +  B(g)  =  AB(g) 

(40a) 

Ap?  »  -RT  in  PAB 

1  papb 

(40b) 

o 

b(£) — *B(g) 

(4la) 

/\f°  =  -RT  In  Pg 

(4lb) 

3- 

AB(g) — ^AB(Z) 

(42a) 

Af°  -  RT  U  Pab 

(42b) 

4. 

AB(£)— )  AB(>t) ,  dilute  solution  in  B(£) 

(43a) 

tV°  -  -®h». m 

(43b) 

The  sum  of  the  standard  free  energy  changes  for  these  reactions  is 


Af°  =  RT  In 


%B  kAB 


m 


The  free  energy  changes  given  by  Eqs.  39  and  44  would  be  equal 
if  Pa/Nab  kAB  of  Eti''  39  is  equal  to  Pa/nab  of  44,  or  if  % B  is 
unity.  If  the  perfect  gas  law  is  used  to  estimate  the  number  of  moles 
of  A,  iu,  in  the  gas  volume  from  pa  in  Eqs.  39  or  44  and  replacing  Nab 
by  xipg/n(i),  the  two  equations  become 
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nA 


=  &A  kAB  KAB  e 


/^F°/RT 


(n(£)/V)  RT 


and 


nA 


^f°/rt 

(n(j ty/v)  RT 


=  nA  kAB  e 


(^5) 


(U6) 


Af°/rt/t,  ,  \  .  Af°/rt 

respectively.  Replacing  either  K^g  e  or  e 

(Eq.  46)  by  k^g  reduces  the  two  equations  to  the  same  form  as  Eq.  21. 
The  same  summation  formulae  following  Eq.  21  for  the  material  balances 
would  therefore  apply  except  that  the  standard  free  energy  functions, 
if  available,  could  be  used  to  calculate  the  free  energy  change  for 
the  reaction  given  by  Eq.  35  or  for  those  of  the  various  separate 
reactions . 

The  reaction  described  by  Eq.  35  did  not  include  any  reaction 
between  the  atmospheric  oxygen  and  the  element  A  in  the  condensation 
process.  Of  course,  even  in  the  presence  of  oxygen  many  of  the  fis¬ 
sion  product  oxides  are  partially  or  completely  dissociated  in  the 
vapor  phase.  For  those  that  are  completely  dissociated  at  the  tem¬ 
peratures  where  the  carrier  material  exists  as  a  liquid,  the  reaction 
described  by  Eq.  35  is  applicable.  Also,  if  the  oxides  are  completely 
associated,  A(g)  may  be  taken  to  represent  the  oxide  molecule,  and 
the  reaction  of  Eq.  35  can  be  applied.  However,  for  the  elements  in 
which  the  oxide  molecule  is  partially  dissociated  in  the  vapor  state 
and  associates  further  with  oxygen  in  the  condensed  state,  the  oxygen 
partial  pressure  will  influence  the  relative  amounts  of  the  element  in 
each  state. 

The  overall  reaction  for  this  condensation  process  is 


A(g)  +  xC^g)  +  B(jt) — *A02x  B(£),  dilute  solution  in  B(^)  (47) 


in  which  x  is  the  number  of  oxygen  molecules  that  combine  with  each 
atom  of  element  A.  Separate  reactions  for  this  overall  reaction  can 
be  set  up  in  the  same  way  as  for  Eq.  35*  For  the  reaction  of  Eq.  47 
to  be  different  from  Eq.  35  (or  competing  with  it),  the  gas  atoms  of 
A  and  02  must  be  in  equilibrium  with  the  oxide  of  element  A  in  the 
vapor  state.  The  oxide  molecules  then  either  react  with  B(/)  as 
they  condense  or  react  with  B(g)  in  the  vapor  and  then  the  larger 
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molecules  condense  to  form  the  dilute  solution.  In  the  first  pro¬ 
cess,  the  separate  reactions  and  the  standard  free  energy  changes 
are: 


1*  A(g)  +  x02(g)— >A02x(g) 

(48a) 

Af°  =  RT  in  Pao/paPq 

(48b) 

2.  A02x(g) - >A02x{L) 

(49a) 

AF°  =  RT  In  pAQ 

(49b) 

3.  AO^Ol)  +  B(Jt) 

(50a) 

Ap°  -  -RT  A  kaob 

(50b) 

where  K^q-d  is  the  equilibrium  constant  for  the  formation  of  AO^  B (X) 
from  the  two  liquid  compounds  at  the  temperature,  T. 


4 .  A02x  — ^A02x  b(£)>  dilute  solution  in  b(,£)  (51a) 

Af°  -  RT  in  Naob  kA0B  (51b) 


The  sum  of  these  reactions  and  standard  free  energy  changes 
give,  for  Eq.  47, 


/\F°  -  RT  in  2^2 - - 

nAOB  kA0B  KAOB 


(52) 


In  the  second  process,  the  separate  reactions  and  standard 
free  energy  changes  are : 


1.  A(g)  +  x02(g)-7>A02x(g) 

Af°  -  -RT  in  PA0/PAP o 


(53a) 

(53b) 
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2. 

A02x(g)  +  B(g) — »A02x  B(g) 

(54a) 

AF°  =  -RT  /n  Paob/Pa^B 

(54b) 

3- 

b(jl)  y  B(g) 

(55a) 

AF°  =  pb 

(55b) 

4. 

A02x  B(g) — >  AOgx  BU) 

(56a) 

Af°  -  HT  U  pA0B 

(56b) 

5- 

A02x  B(JL)  — >  AOg*  B (i),  dilute  solution  in  B(£) 

(57a) 

Af°  -  -HI  £n  Baob  tAaB 

(57b) 

The  sum  of  the  five  standard  free  energy  changes  for  Eq.  47  is 


/\p°  ■  m£n  PAP° _ 

aOB  aob 


(58) 


Equations  52  and  58?  when  solved  for  n.  "by  use  of  the  perfect  gas  lav 
give 


and 


Af°/rt 

nA  =  nAOB  kA0B  KAOB  e _ 

(n U)/V)  RT  pg 

Af°/rt 

nA  a  nAOB  kA0B  e 

(n (l)M  RT  p£ 


(59) 


/ 

\V\>) 


respectively.  The  values  of  n^  in  Eqs.  59  and  60  are  sensitive  to 
the  oxygen  partial  pressure.  If  Pq  is  greater  than  1  atmosphere* 
the  value  of  will  "be  decreased  (most  with  largest  x  value)  and 
the  amount  of  element  A  condensed  is  increased.  If  Pq  is  less  than 
1  atmosphere*  the  value  of  nA  will  be  increased  (most  with  largest 
x  value)  and  the  amount  of  element  A  condensed  is  decreased.  If 
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X  o 

Pq  is  incorporated  into  the  k^oB  in  ‘t^ie  Bame  way  aB  was  done  for  the 
factors  of  Eqs.  45  and  46,  then  Eqs.  59  and  60  can  also  he  reduced 
to  the  same  form  as  Eq.  21  for  material  balance  and  other  general 
summation  formulae. 

Equations  35  to  60  are  perhaps  more  rigorous  in  the  definition 
of  the  condensation  process  than  those  of  Sections  2.2  and  2.3  that 
describe  the  process  only  in  terms  of  Henry's  law.  In  a  real  case*, 
the  carrier  material  is  not  an  inert  substance  but  capable  of  forming 
compounds  with  many  of  the  fission  product  elements.  Also,  the  above 
equations  will  be  applicable  for  elements  that  react  with  the  carrier 
in  the  solid  state  (see  Section  2.6);  for  these,  the  notations  only 
need  be  changed  from  ($.)  to  (3)  to  refer  to  the  solid  rather  than 
the  liquid  state. 


2.2  THE  SECOND  PERIOD  OF  CONDENSATION 


2.2.1  Description  of  the  Formation  Process 


In  the  second  period  of  condensation,  the  carrier  material  is 
in  the  solid  state  and  the  condensation  of  the  remaining  radioelements 
can  only  occur  on  the  surface  of  the  carrier  particles.  The  latter 
can  consist  of  unmelted  particles  and  the  smaller  particles  that  had 
been  melted  but  have  not  fallen  out  of  the  cooling  fireball.  Although 
this  period  of  condensation  may  continue  indefinitely  with  the  exchange 
of  rare  gas  elements  that  are  formed  by  decay  or  that  decay  into  alk¬ 
ali  metal  daughter  elements,  it  should  he  essentially  complete  with 
respect  to  most  of  the  radioactive  elements  between  5  and  7  minutes 
after  detonation.  For  a  given  group  of  particles,  it  should  he  com¬ 
plete  when  they  fall  away  from  the  volume  containing  the  radioactive 
gases.- 


Of  the  elements  that  condense,  or  adsorb,  on  the  surface  of 
the  larger  particles,  some  consist  of  decay  products  of  the  more 
volatile  elements  (present  at  the  end  of  the  first  period  of  conden¬ 
sation)  and  would  tend  to  be  soluble  in  rain  water  and  salt  solutions. 
Others,  consisting  of  the  more  reactive  elements,  may  form  compounds 
with  the  solid  carrier  and  even  diffuse  into  the  interior  of  the  par¬ 
ticles.  In  either  case,  the  elements  condensing  during  this  period 
would  be  more  available  for  solution  than  those  condensed  into  the 
liquid  phase. 

Since  a  surface  condensation  is  indicated  for  this  process, 
the  effective  specific  activity  (i.e.  activity/gram)  carried  by 
various  particle  sizes  should  be  inversely  proportional  to  the 


24 


particle  "diameter".  In  a  sample  of  fallout  containing  particles 
from  both  types  of  condensates,  the  measured  specific  activity  would 
not  correspond  to  any  of  the  curves  of  Fig.  1;  the  flat  portion  of 
the  curve(s)  indicated  for  the  smaller  particles  would  instead  show 
an  increase  in  the  value  of  n^(p)/MjnQ(/,p)  as  the  particle  diameter 
decreased. 

2.2.2  Condensation  Mathematics 


The  mathematical  equations  for  describing  the  processes  that 
can  occur  in  this  period  of  condensation  have  not  been  fully  developed. 
Completion  of  the  current  work  by  C.E.  Adams^-1  should  aid  the  treat¬ 
ment.  However,  tentative  treatment,  of  two  rather  simplified  types 
of  interactions  allow  approximate  formulae  to  be  given.  For  the 
less  volatile  elements  and  a  large  amount  of  available  surface,  the 
simplified-  form  of  Langmuir's  adsorption  isotherm  (low  partial  pres¬ 
sures)  is  given,  in  the  terms  used  previously,  by 

nj  =  aj(e)nj  (6l) 

in  which  n!  is  the  amount  of  element  adsorbed,  a^(s)  is  a  constant 
including  the  surface  area,  and  n^  is  proportional  to  the  partial 
pressure  of’ the  element  j.  The  value  of  the  constant,  aj(s),  would 
presumably  depend  on  the  kind  of  gas  molecule  that  is  adsorbed.  If 
n°  is  to  be  the  same  as  that  given  by  Eq.  21  for  the  amount  of  ele¬ 
ment  j  remaining  at  the  end  of  the  first  period  of  condensation, , 
then  the  adsorption  reaction  described  by  Eq.  6l  must  occur  within 
a  very  short  time.  Otherwise  its  value  must  be  adjusted  to  account 
for  decay.  This  adjustment  would  not  be  required  if  all  the  ele-  , 
ments  in  the  chain  are  adsorbed  and  remain  so.  For  a  short-time 
condensation,  n°  may  be  replaced  by  kjn^  and  the  "R"  factor  for 
this  condensation  process  can  be  written  as . 

r'(A)--  b  °J(s)  kJ  m 

Lj  aj(s)  kj  yj(99>t)/(l+kj) 

If  aj(s)  is  the  same  for  all  elements,  Eq.  6 2  reduces  to 


r’(A)  «  ^j  k3  yj(A,t)/(l+fc3) 
yj(99,t)/(l+k°) 


(63) 
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Although  the  sum  of  the  numerators  and  denominators  in  Eqs.  34 
63  will  both  be  equal  to  1,  the  sum  of  r(A)  and  r‘(A)  will  not  equal 
1  (unless  one  is  aero). 

Another,,  and  probably  more  realistic,  process  is  to  consider 
the  adsorption  condensation  as  one  in  which  the  relative  amount  of 
each  element  condensed  is  related  to  its 'sublimation  pressure  (as  if 
the  solid  carrier  would  act  as  pure  compound  of  each  of  the  condens¬ 
ing  species) .  The  computational  values  from  such  a  process  will 
reflect  the  relative  volatility  of  constituent  molecules  at  all 
temperatures  at  which  this  kind  of  condensation  can  occur.  If  an 
excess  of  solid  surface  area  is  present,  the  amount  condensed  by 
the  process  (assuming  it  to  be  reversible)  at  any  time  after  solid¬ 
ification  of  the  carrier  is  given  by 


00 


(64) 


in  which  n*  is  the  amount  of  element  j  condensed  on  the  surface  of 
the  solid  particles,  n?  is  the  amount  of  the  element  that  had  not 
condensed  in  the  liquid  carrier,  and  n^°  is  the  amount  in  the  vapor 
phase.  All  three  quantities  depend  on  the  time  after  solidification 
because  of  decay.  For  a  perfect  gas 


(65) 


in  which  V  is  the  volume  containing  the  n°,°  moles  of  the  gaseous 
species  and  p®  is  the  sublimation  pressure  and  is  given  by 

p».  eAS0/B  e-AHs/ffl  (66) 

in  which  Asg  is  the  entropy  of  sublimation  sndAHs  is  the  heat  of 
sublimation  at  the  temperature,  T.  If  uniform  mixing  in  the  fireball 
volume  is  assumed  for  the  particles  and  gaseous  species,  then  V  is 
the  fireball  volume.  The  material  balance  for  element  j  (all  time- 
dependent)  is  1 

■  j  ' 

Yj  -  n^  +  nj  +  (67) 

so  that,  Eq.  64,  with  Eqs.  21,  65,  and  67,  becomes 
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(68) 


or 


n*  =^1 


ring 


53 

RT 


njA- 


ia_E 


i  +  k3 


A  yjA 


_  v_  r  s 

“  RT  ^A  PjA 


(69) 


Separating  out  the  equal  terms  for  each  nuclide  gives 

i  _  kjyJA  -  VP,jA 
iA  "  ,  ,  o  RT 

The  "R"  factor  for  this  process  is  now 

_  -  (v/RTiZj  (A,t) 

Li  yj(99,t)k°/(i+k°)  -  (v/rt)£j  Pj(99>t) 


(TO) 


(71) 


in  which  pj(A,t)  is  equal  to  PjaAa  an^  yA  'to'tal  chain  yield. 

For  elements  with  very  low  sublimation  pressures,  Eq.  71  reduces  to 
Eq.  63.  For  chains  with  very  volatile  elements,  r’(A)  could  approach 
zero . 


2.3  SUMMARY  OF  FALLOUT  PARTI CLE  FORMATION  PROCESS 


A  general  description  of  the  over-all  fallout  particle  formation 
process,  from  the  foregoing  treatment,  is  as  follows.  The  larger 
particles  that  fall  out  of  the  fireball  earliest  and  land  near  ground 
zero  should  contain  radioelements  that  were  condensed  during  the  liq¬ 
uid  phase  and  very  little,  if  any,  radioelements  that  were  condensed 
during  the  second  period  of  condensation.  The  smaller  particles  that 
enter  late  and  stay  in  the  rising  fireball  and  cloud  for  longer  times 
should  carry  radioelements  that  were  condensed  on  their  surfaces. 

The  smallest  of  these  particles  would  make  up  the  world-wide  fallout 
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and  the  others  are  deposited  at  large  distances  from  ground  zero.* 
The  intermediate  size  particles  that  deposit  at  intermediate  dis¬ 
tances  from  ground  zero  should  contain  radioelements  that  were 
condensed  during  both  periods  of  condensation.  It  is  reasonable 
to  assume  that  individual  particles  at  these  locations  could  con¬ 
tain  condensates  from  either  or  both  types  of  condensation.  In  a 
radiochemical  determination  of  an  "R"  factor  for  a  given  mass  chain 
from  a  gross  sample  of  the  fallout  in  which  a  given  fraction  of  a 
mass  chain  had  condensed  during  the  first  period  of  condensation 
and  the  remainder  during  the  second  period;  the  expected  "ft”  factor 
would  be  given  by 


R(A)  =  k(a) 


(1-x)  GjA  x 

(1-x)  £jnj99  +  x 


(72) 


in  which  x  is  the  fraction  of  element  j  condensed  during  the  second 
period  of  condensation.  Substitution  of  the  appropriate  equalities 
in  Eq.  -72  gives 

yj(A,t)/(l+k°)  -  x[jC3  yj(A;t)(l-kj)/(l+kj) 

Ej  yj(99>t)/(l+k°)  -  x[Zj  (99,t) (l-k°)/(l+k° ) 

(73) 

+  (V/RT)Xj  Pj(A;t)j 
+  (v/rt)£\  Pj(99;t)] 

The  variation  of  x  from  0  to  1  corresponds  respectively  to  the  close- 
in  and  world-wide  (or  far  distant)  fallout  compositions  discussed 
above. 


In  all  the  condensation  processes,  no  variation  of  fractionation 
with  particle  size  (for  sizes  larger  than  a  few  tenths  of  a  micron) 
was  indicated.  Any  variation  in  the  amounts  condensed  to  give  a  vari¬ 
ation  in  fractionation  with  particle  size  would  be  caused  by  a  com¬ 
bination  of  time  of  particle  fall  from  the  gas  volume  and  the  value 
of  x.  It  is  possible  that  some  amount  of  the  larger  particles  would 
fall  o.ut  of  the  gas  volume  or  move  through  it  like  projectiles  at 


*The  present  world-wide  fallout  from  air;  seawater,  tower,  and  sur¬ 
face  bursts  also  contains  a  large  fraction  of  vapor-condensed  par¬ 
ticles  with  the  activity  more  or  less  uniformly  distributed  through 
their  volumes. 


28 


temperatures  (and  times)  between  the  boiling  and  melting  point  of 
the  carrier.  The  description  of  this  computationally  as  a  continu¬ 
ous  process  via  Eq.  73*  for  example,  would  be  to  initiate  the  cal¬ 
culations  for  the  phase  of  condensation  at  the  boiling  point  of  the 
carrier  (with  x  =  0  until  the  melting  point  is  reached)  and  continue 
the  process  up  to  the  time  of  6  or  7  minutes  after  detonation.  Such 
a  process  should  give,  in  relative  terms,  the  overall  variation  of 
fractionation  with  size  and  downwind  distance  from  ground  zero. 


SECTION  3 


THE  VAPOR  PRESSURE  GF  SOME  OXIDES,  COMPOUNDS,  AND 
ELEMENTS  OF  FISSION  PRODUCTS  AND  OTHER  SUBSTANCES 


If  the  Henry’s  law  constants  for  all  the  radionuclides  and  a 
silicate  soil  material  (as  the  carrier)  were  known,  the  condensation 
and/or  fractionation  "R”  factors  of  Section  2  could  he  calculated, 
also  provided  some  estimate  of  the  term  (n(j£)/v)RT  could  he  made. 

The  necessary  data  has  not  yet  heen  measured.  Estimates  of  the 
Henry’s  law  constants  can  he  made,  however,  if  appropriate  proper¬ 
ties  are  assumed  for  the  carrier  material  and  the  gaseous  fission 
products.  For  such  calculations,  presented  in  Sections  4  and  7, 
the  following  properties  were  stipulated:  (l)  the  gaseous  species' 
of  each  fission  product  element  (with  a  few  exceptions)  was  the 
same  as  that  in  equilibrium  with  its  own  liquid  (or  solid)  oxide 
in  the  presence  of  oxygen;  (2)  the  liquid  mixture  formed  a  perfect 
( ideal)  solution;  (3)  the  carrier  material  was  non-reactive  so  that 
no  compound  formation  was  involved  (alternately,  if  a  compound  was 
formed  its  vapor  pressure  was  the  same  as  the  oxide);  and  (4)  the 
melting  point  of  the  carrier  was  taken  to  he  l400°C . 

For  the  formation  of  one  mole  of  gas  from  either  the  liquid 
(vaporization)  or  solid  ( sublimation) ,  the  'standard  free  energy  is 
given  by 

Af°-=-Ah°  -  tAs°  (74) 


in  whichAH0  is  the  difference  in  heat  content  of  the  gas  and  the 
liquid  or  solid,  and  A 3°  is  the  change  in  entropy.  Also, 


Af°  =  -RtAiK  =  -RTinp 


(75) 
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in  which  p  is  the  vapor  pressure  of  the  gas  in  contact  with  the  liquid 
or  solid  and  K  is  the  equilibrium  constant  for  the  vaporization,  or 
sublimation  process.  Combining  Eqs.  57  and  58  gives 


in  which 


and 


log  P  =  -  |  +  B 


A 


B 


-  Ah° 

”  1^57^ 


-  As° 

“  ^^575 


(76) 

(77) 

(78) 


The  notation  Ay,  By  will  be  used  for  vaporization  and  As,  Bs  for  sub¬ 
limation.  In  general,  both  Ah°  and  As0  vary  with  temperature  so 
that  A  and  B  are  not  actually  constants.  They  were  considered  con¬ 
stant  here  however  and  evaluated  from  the  data  at  the  boiling  point 
(pel  atmos)  and  the  temperature  at  -which  the  total  vapor  pressure 
is  10"3  atmos.  For  the  oxides  that  decompose  on  vaporization  or 
sublimation,  the  constants  were  then  adjusted  to  a  constant  oxygen 
pressure  of  1  atmos.  For  example,  sodium  oxide  decomposes  to  Na(g) 
and  02(g)  as  follows; 

2  Na20?=i4Na(g)  +.  Og(g)  (79) 


For  this  reaction, 

p  =  p(Na)  +  p(02)  (80X 

and 

p(Na)  =  4p(02)  (81) 

so  that 

p  -  5p(o2)  (82) 
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(83) 


Also  , 

K  =  p4(Na)  p(02)  *  pj(l) 


where  pj_  is  the  sodium  vapor  pressure  when  the  oxygen  pressure  is  1 
atmosphere.  Combining  the  above  equations  gives 


log  px  =  log  4  -  (.5/4)  log  5  +  (5A)  log  p  (84) 
log  px  =  -  0.272  -  I.25A */T  +  1.25B'  (85) 


For  oxygen  pressures  other  than  1  atmos.,  p(Na)  can  be  determined 
from  the  ratio  Pl/p1/^(02) •  The  constants,  for  Eqs.  76,  via  equations 
similar  to  Eq>  85,  and  the  various  vaporization  and  sublimation  reac¬ 
tions  and  melting  points  for  the  oxides  or  other  assumed  gaseous 
species  are  summarized  in  Table  1.  Some  values  of  the  entropies  of 
vaporization  and  sublimation  and  even  of  the  melting  points  and  boil¬ 
ing  points  were  obtained  from  correlations  or  were  otherwise  estimated, 
Most  of  the  data,  however,  should  give  values  fairly  close  to  the  true 
pressures  between  10 “3  and  1  atmcs , ,  as  well  as  for  pressure  of  sev¬ 
eral  atmospheres  for  temperatures  higher  than  the  boiling  points. 

For  most  mass  chains  in  the  condensation  process,  this  range  of  vapor 
pressure  is  the  most  critical  one  if  the  melting  point  of  the  carrier 
material  falls  within  the  corresponding  temperature  range.  If  the 
vapor  pressure  of  the  element  or  its  oxide  is  either  very  low  or 
very  high  at  the  carrier  melting  point  temperature,  the  fraction  con¬ 
densed  will  he  essentially  one  or  zero,  respectively;  the  exact  value 
of  the  pressure  is  then  not  required. 

For  a  perfect  solution,  the  relation  between  the  vapor  pressure 
of  the  solute  and  its  mole  fraction  is  given  by  Raoult's  law  which, 
in  terms  of  the  notation  given  here,  is 


Pj  «  N.p0  (86) 

where  p0 , is  the  vapor  pressure  of  the  pure  compound  and,  for  a  per¬ 
fect  solution,  is  the  same  as  kj  for  Henry's  law.  Ko  estimate  can 
be  given  for  the  magnitude  of  the  deviations  that  might  be  expected 
from  perfect  solution  behavior  in  the  dilute  solutions  at  high  tem¬ 
peratures  in  the  silicate  systems  of  all  of  the  fission  product 


33 


TABLE  1 


.Vapor  Pressure  Data  for  the  Ox  idea  of  Fission  Products  and  other  Elements 
Constanta  for  log  pa  »  (a/t)  +  B 


Oxide 

Reaction 

M.P. 

(°K) 

*v 

»v 

pJ 

KagO 

%o 

Kb?0 

Cs^O 

a +  o,(g) 

1193 

mr 

763 

l4, 320 
12,600 
13,190 
11,250 

6.75 

6-93 

7*97 

7.23 

16,060 

14,000 

14,440 

12,360 

8.21 

8.38 

9.42 

8,68 

Pi/p^Coa) 

HgO 

CaO 

SrO 

BaO 

KO-»«o(e) 

3075 

2860 

2730 

2196 

21,830 

24,750 

20,1*50 
20, 980 

6.51 

6.51 

5.84 

6.99 

29,090 

31,520 

26,890 

26,100 

8.87 

8.&I 

8.20 

9*35 

pi 

CdO 

FeO 

2M0->2M(b)  +  02(g) 

dl750b 

1641* 

l6,46o 

31,90) 

8.98 

8.98 

19,  no 
24,500 

10.50 

10.50 

p1/p1/2(o2) 

><30^-1  3H(g)  +  20, 

1870 

84,350 

40.44 

96,35c 

46. 85 

Pl/p2/3(o2) 

y2°3 

l^2^3 

rr&i 

Hd^Oj 

p=eo3 

2»20j~>  liMO(g)  +  02(g) 

2690 

2590 
-  (2420) 
(2380) 
(2300) 

1*0,500 

39,700 

38,900 

30.500 

37.500 

8.61 

(8.61) 

(8.611 

(e.6i) 

(8.61) 

56, (too 
50,400 
40,800 
48,300 
47,000 

15.33 

(15-33) 

15.33) 

15.335 

(15*33) 

Pi/p^tOa) 

Aig03. 

Ina03 

«203_,M20(g)  ♦  02(8) 

2300 

1998 

1560 

58,000 

44,300 

31,600 

14.868 

( 14*87) 

73,900 

41,800 

21.25 

Pl/P(0g) 

A“2°3 

ay>j_>llMO(g)  ,  OjCg), 
t>i4oo°k 

588 

8,700 

8.61 

- 

■ 

Aa2°3 

588 

3,670 

5.01 

5,4oo 

7.96 

Pi 

8bg°3 

2Kjp3^  kMO(e)  ,  02(g), 
T>2000°g 

928 

16, 4oo 

8.61 

_ 

. 

Sb203 

ayj->Hi,o6(e) 

928 

9,170 

5.01 

11,900 

7*96 

Cu 

Pd 

Afl 

«->H(g) 

1356 

1823 

1234 

24,l*CO 

19,600 

13,3*0 

9.45 

5-70 

5.44 

20,670 

6.29 

pi 

2.rOa 

CeOo 

UOg 

M^-gMC^g) 

2960 

(3000) 

2600 

32,360 

29.500 

26.500 

6.73 

St:Si 

36,950 

34,200 

30,700 

8.28 

as 

U300 

2M308  ♦  03(g)_> 
(m-jie) 

(2450) 

11,030  • 

3-55 

19,600 

7.05 

Plp1/6(oi!) 

SlOo 

GeOg 

3n0g 

2K0a->2M0(g)  t  05(g) 

(2040) 

1389 

1098 

39,250 

20,900 

26,750 

12.55 

11.85 

11.91 

39,700 

22,730 

29,230 

12.79 

13.15 

13.22 

p/pi/?(o?) 

RuOo 

RbO~ 

«  +  Ojtg) ->l«2(8) 

2700 

2240 

5,330 

10,300 

0.52 

0.44 

6,670 

11,670 

1.02 

0-95 

PiP<0p) 

Ru 

Rh 

M-vM(g) 

2700 

2240 

?<;?» 

IM 

3;Z88 

?:& 

pi 

SeOg 

a»j-^3IO(g)  .  03(e), 
T>1500°K 

613 

13,060 

12.74 

14,030 

14.29 

p/p1/2(02) 

SeOg 

"Og-lHOjtg) 

613 

3,970 

6.78 

4/920 

8.33 

Pi 

TeOg 

aBj-vawte)  .  Oj(e), 

T  >  10OO°K 

1006 

18,500 

11.35 

19,850 

12.68 

.p/pl/e(02> 

TeOg 

H02-,«02(g) 

1006 

11,300 

7  •  37 

12,000 

8.07 

Pi 

Mb2Q5 

“2°5  ->>¥>5(8) 

1785 

22,300 

6.75 

27, 700 

9.77 

KsO-3 

*°3  -*  *305(8)  ,T  <  J-5T5 

1068 

7,620 

5.64 

17,430 

i4,8e 

h*°3 

BOj-^KOjtj) 

1068 

10,200 

6.66 

13,590 

9.73 

Tc^ 

^-.HgOfts) 

’39? 

3,570' 

6.12  — 

7,205 

15 .40 

MaBr 

Nal 

HX— >HX(g) 

1028 

924 

8,830 

6, 

5.30 
-  5.50 

- 

- 

Kr 

Xe 

x-*x(g) 

485 

524 

4.o4 

4.00 

a.  Values  In  parentheses  are  estimated  values, 
to.  d  for  decomposition  point. 
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elements.  For  many  of  them,  the  deviations  should  certainly  "be 
large.  In  assuming  a  liquid  silicate  as  the  condensing  surface, 
then  those  elements  whose  ionic  character  differs  most  from  bi+^ 
will  probably  deviate  most  from  ideal  behavior.  The  elements  that 
would  be  expected  to  deviate  most  would  therefore  be  those  of  the 
alkali  metals.  However,  similar  oxides  should  give  similar  devi¬ 
ations  so  that  the  relative  values  of  pj  for  similar  oxides  should 
be  in  the  same  numerical  order  as  the  kj  values.  The  fact  that 
the  solutions  are  dilute  does  not  necessarily  indicate  that  the 
solution  is  perfect  or  make  a  good  estimate  of  kj. 

The  melting  point  selected  is  approximately  that  of  a  silicate 
soil  containing  fairly  large  amounts  of  AlgO-^  and  KgO  or  NagO  along 
with  SiOg.  For  a  surface  detonation,  large  amounts  of  this  material 
would  be  vaporized.  The  AlgOcj,  with  a  boiling  point  of  3950°K  (l 
atmos-  Og),  could  start  vapor-condensing  at  some  temperature  leBS 
than  this  (depending  on  the  partial  pressure)  into  liquid  drops 
and  co-condense  some  of  the  less  volatile  radioelements.  The  latter 
cannot  vapor-condense  by  themselves  because  their  actual  vapor  pres¬ 
sures  are  lower  than  those  which  would  be  calculated  from  the  data 
in  Table  1.  These  small  vapor-condensed  particles  of  AlgO^  will 
solidify  when  they  cool  to  2300°K  and/or  combine  with  SiC>2  alu- 
mino-sllicate  liquid  particles  that  can  be  formed  between  about 
3100°K  and  2040°K.  The  influx  of  larger  liquid  particles  at  these 
and  lower  temperatures  can  dissolve  the  very  small  vapor -condensed 
particles  upon  collision  as  well  as  condense  more  of  other  elements 
still  in  the  gaseous  state.  The  establishment  of  equilibrium  vapor 
pressures  for  elements  that  decay  from  a  non-volatile  to  a  volatile 
element  in  the  liquid  state  may  be  slower  than  for  the  gas-to-liquid 
equilibrium  since  diffusion  through  the  liquid  to  the  surface  is 
slower  than  diffusion  through  the  gas  to  the  liquid  surface. 

A  carrier  material  such  as  a  silicate  soil  contains  elements 
and  compounds  that  could  form  other  compounds  with  the  fission  pro¬ 
duct  elements  whose  vapor  pressures  and  solubilities  in  the  liquid 
state  could  be  different  from  that  calculated  for  the  oxide  (aside 
from  use  of  Raoult's  law).  For  example,  Na  and  K  form  compounds 
with  all  reactive  non-metals  and  most  metallic  elements  form  sili¬ 
cate  compounds.  However,  many  also  form  silicate  glasses  as  solu¬ 
tions  and  for  these  the  assumed  process  applies.  The  vapor  pressure 
data  for  the  possible  gaseous  species  of  the  fission-product  elements 
over  silicate  solutions  have  not  been  measured.  The  exact  gaseous 
species  need  not  be  known  if  the  Henry's  law  constants  for  the  liquid 
solutions  are  known  as  a  function  of  temperature. 
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SECTION  4 


EVALUATION  OF  THE  u(j£)/V  RATIO 


4.1  CALCULATION  OF  THE  RATIO  a (£)/V  FROM  RADIOCHEMICAL  DATA 


Radiochemical  data  from  samples  of  fallout  that  give  the  content 
of  out;  or  more  radionuclides  per  gram  of  fallout  (preferably  consist" 
ing  of  the  larger  fused  particles)  can  he  used  to  make  estimates  of 
the  concentration  of  the  carrier  material  in  the  fireball  during  the 
first  period  of  condensation.  In  order  to  do  this,  however,  the 
assumption  of  uniform  mixing  in  the  fireball  must  be  made,  the  vol¬ 
ume  of  the  fireball  at  the  time  when  the  particles  solidified  must 
be  known,  and  the  total  yield  of  the  radionuclide(s)  (i.e.  total 
chain  yield)  is  required.  Also,  the  molecular  weight  of  the  carrier 
material  must  be  known.  The  mole  fraction,  Nj,  can  then  he  calcu¬ 
lated  in  terms  of  the  ratio  of  the  number  of  moles  of  the  mass  num¬ 
ber,  A,  to  the  number  of  moles  of  the  carrier.  With  this  information, 
the  ratio,  n(/)/v,  can  he  estimated  from 


yj(A,t) 

Ti  i 
1  *  X’ajrm 


(87) 


in  which  YA  is  the  total  weapon  yield  of  the  mass  number  A  and  V  is 
the  volume  of  the  fireball  at  the  time  the  carrier  melting  point 
occurs.  The  best  estimates  of  n(£)/v  will  result  from  analyses  for 
mass  chains  that  are  neither  condensed  in  large  amount  nor  missing 
from  the  sample  in  large  amount.  Equation  87  can  be  solved  by  suc¬ 
cessive  approximations  and,  if  data  for  several  mass  numbers  axe 
available,  the  final  value  of  the  ratio  should  be  selected  from 
those  in  which  the  summed  term  has  values  between  0.2  and  0.8. 
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Ho  complete  set  of  unclassified  data  have  been  reported  for  use 
of  Eg..  87.  However,  Kimura^  gives  values  of  the  specific  activity 
of  the  fallout  from  the  1  March  195^  detonation  at  the  Eniwetok  Prov¬ 
ing  Grounds.  His  values  are  0.37  mc/gm  of  fallout  on  23  April  (IH*53) 
and  8  x  10"°  gm  fission  products/gm  of  fallout.  In  terms  of  unfrac¬ 
tionated  fission  products  where  the  activity  at  53  days  is  3-7  x  10"° 
d/s/fission, ^3  the  first  value  can  he  converted  to  7  x  10-^  fission 
product  atoms  per  gram  of  fallout.  The  second  value  converts  to 
4  x  1CP-1*  fission  product  atoms  per  gram  of  fallout.  The  report  also 
gives  data  on  the  relative  yields  of  some  of  the  individual  radio¬ 
nuclides  present.  However,  the  carrier  material,  CaOj  litis  a  much 
higher  melting  point  than  siliceous  soil  minerals  and  the  end  of 
the  first period  of  condensation  would  occur  much  earlier  for  the 
CaO  carrier,  for  the  yield  of  that  device.  The  gross  fission  pro¬ 
duct  concentrations  given  above  indicate  at  least  the  carder  of  mag¬ 
nitude  of  the  specific  activity  of  fallout  from  large  yield  weapons. 
The  information  in  the  remainder  of  this  Section  will  show  that  large 
variations  in  the  specific  activity  in  the  fallout  from  weapons  of 
different  yield  is  not  expected  (i.e.  for  surface  detonations)  if 
the  fraction  of  fission  yield  is  held  constant.  It  should,  however, 
vary  directly  with  the  fraction  of  fission  yield.  For  example,  if 
the  1  March  1954  detonation  were  50  percent  fission,  then  a  fission 
weapon  might  be  expected  to  produce  fallout  with  specific  activities 
between  8  x  1011*'  and  14  x  1CH-^  fisBion  product  atoms  per  gram  of 
fallout . 

According  to  data  from  The  Effects  of  Nuclear  WeaponsfEHW) , 
the  crater  radius  for  a  surface  detonation  is  62W1/3  feet  (W  in  KT) 
and. its  depth  is  25 ft.  For  a  cone-shaped  crater,  which  is  more 
representative  of  the  crater  shape  for  a  large  yield  detonation  than 
is  an  ellipsoid  of  revolution  (as  given  in  and  a  soil  dens¬ 

ity  of  110  lb/cu  ft,  the  mass  of  material  thrown  out  of  the  crater 
is  given  by 


M  -  5-03  x  lO^W0*92  gm 


(88) 


If  all  the  activity  produced  in  a  20-KT  fission  weapon  were  mixed 
uniformly  with  the  ejected  mass,  the  fission -product  concentration 
would  be  about  7-3  x  10^3  atoms/gm;  for  a  1-MT  fission  weapon,  it 
would  be  about  1.0  x  lO1^  atoms /gm;  and  for  a  15-MT  fission  weapon, 
it  would  be  about  1.2  x  10^.  With  the  assumption  of  uniform  mix¬ 
ing  of  the  fission  products  with  all  the  crater  material,  the  con¬ 
centration  changes  rather  slowly  with  yield  (i.e.,  only  as  W®*®”). 
However,  even  for  the  15-MT  yield,  the  "uniform  mixing"  concentra¬ 
tions  are  lower  by  a  factor  of  7  to,  12,  assuming  50  $  fission  for 
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the  1  March  195^  detonation,  or  3  to  6,  assuming  100  $  fission,  than 
those  found  by  Kimura.  Thus  these  percentages  indicate  that  only 
the  order  of  5  to  10  percent  of  the  crater  material  must  come  in 
contact  with  the  radioactive  material  and  is  involved  in  the  for¬ 
mation  of  radioactive  fallout  particles . 


k.2  CALCULATION  OF  THE  RATIO  n(/)/v  FROM  THERMAL  DATA 
AND  EMPIRICAL  SCALING  EQUATIONS 


4.2.1  Description  of  Calculation 

The  scaling  equations  and  other  data  for  an  air  burst  given 
in  BNW  will  be  used  along  with  modifications  and  idealized  scaling 
assumptions  for  a  description  of  the  fireball  volume  and  tempera¬ 
ture,  with  time,  for  a  surface  detonation.  The  air  burst  model 
assumes  detonation  at  sea  level  without  an  air-soil  interface.  For 
the  surface  detonation,  the  air-soil  interface  will  be  added.  The 
reference  time  of  the  second  maximum  in  the  observed  temperature 
used  in  ENW  will  be  retained  as  a  reference  point  in  time. 

In  the  two  model  detonations,  a  certain  fraction  of  the  total 
energy  released  will  be  estimated  as  still  being  in  the  fireball 
at  the  second  maximum  in  the  form  of  gaseous  atoms  of  high  internal 
energy.  Some  of  the  energy  is  still  in  the  blast  wave  (or  trans¬ 
ferred  to  substances  outside  the  fireball  volume  by  the  blast  wave). 
Another  portion  of  the  energy  has  been  radiated  from  the  fireball 
(EKW  calls  this  portion  the  "thermal,  energy") .  Some  of  the  energy 
has  been  used  in  expansion  of  the  hot  gases  against  the  atmosphere. 

At  times  after  the  second  maximum,  the  gases  in  the  fireball  release 
their  internal  energy  and  cool  by  heating  inflowing  air  and/or  soil 
material,  expanding  further  against  the  atmosphere,  and  radiating 
energy  out  into  the  atmosphere.  The  thermal,  balances  for  these  pro¬ 
cesses  will  be  estimated  in  order  to  compute  the  amount  of  soil  that 
can  be  heated  and  melted  when  the  fireball  has  cooled  to  l4C0°C.  In 
addition,  the  energy  balances  will  be  computed  as  a  function  of  yield 
so  that  the  total  amount  of  the  major  constituents  contained  in  the 
volume  of  the  fireball  can  be  estimated. 

In  the  calculations,  the  energy  used  to  heat  the  surrounding 
air  and/or  soil  is  that  required  to  dissociate  these  substances, 
originally  in  their  normal  thermodynamic  state  at  298°K  and  1  atmos, 
to  atomic  gases  and  then  heat  the  latter  (as  ideal  gases)  to  the 
temperature  of  the  second  maximum.  For  cooling  to  ll00°C,  the 
energy  released  comes  from  the  change  in  kinetic  energy  of  the 
gases  and  from  the  association  of  the  gas  atoms  to  gas  molecules 
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or  to  liquid  molecules  at  l400°C.  Except  for  consideration  of  the 
heat  lost  by  radiation*  ionization  energies  will  net  be  taken  into 
account . 

4.2.2  Partition  of  Energy  at  Second  Maximum 

In  EM  (p.  6),  the  distribution  of  energy  for  the  air  hurst 
is  15  percent  for  nuclear  radiation,  35  percent  for  thermal  radia¬ 
tion  (i.e.  radiant  energy),  and  50  percent  for  blast  and  shock. 

The  time  or  time  period  of  the  detonation  at  which  these  values 
apply  is  not  given,  nor  is  energy  allotted  for  the  processes  given 
in  4.2.1.  Of  the  energy  distribution  values  given  in  Elffi,  the  15 
percent  for  nuclear  radiation  will  be  assumed  to  be  unavailable 
or  lost  from  the  fireball  with  respect  to  the  processes  discussed 
above  for  both  the  air  burst  and  surface  burst.  The  remaining  85 
percent  is,  then,  to  he  distributed  among  (l)  the  energy  in  the 
fireball  at  the  second  maximum,  (2)  the  energy  lost  by  radiance, 
and  (3)  the  energy  carried  beyond  the  fireball  perimeter  by  the 
blast  wave. 

The  energy  content  of  the  fireball  at  the  second  maximum  will 
be  taken  to  include  the  dissociation  energy  for  air  (O.8N2  +  O.2O2), 
the  change  in  internal  energy  for  a  temperature  rise  of  8000°K,  and 
the  work  for  expanding  the  gases  to  the  fireball  volume  at  the  second 
maximum  against  a  1.0  atmosphere  external  pressure.  The  energy  con¬ 
tained  in  the  fireball  for  this  process  is 


=  2nx(3  x  8000  +  79,200)  +  PQV2  -  n-jHT  -  (89) 

in  which  is  the  number  of  moles  of  undissociated  air,  the  value 
3  is  used  for  heat  capacity  (Cv)  in  cals/mile  of  an  ideal  gas,  the 
value  79,200  cal/mole  is  the  dissociation  energy  (  Ae)  for  one  mole 
of  dissociated  gas  atoms  at  298°K,  V2  is  the  fireball  volume  at  the 
second  maximum,  H]_RT0  is  equal  to  P0V0  where  VQ  is  the  original 
volume  of  the  heated  air  molecules,  and  T0  is  taken  as  298°K  (stan¬ 
dard  conditions  at  sea  level) .  For  a  20-KT  yield  air  burst  the 
value  of  V2  as  calculated  from  EM  for  a  spherical  fireball  is  2.5 
X  10-43  cm3  so  that  P0Vg  ^»5  X  10^3  cm3-atmos  or  6.1  X  10-^  cals 
(l  cal  =  41.29  cm3-atmos).  With  these  values,  Eq.  89  reduces  to 


^  =  2.06  X  105ni  +  6.1  x  1011 


(90) 
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The  energy 
p.  333 >  u-P  to  the 
the  conversion  of 
a  20 -KT  air  hurst 


lost  from  the  fireball  by  radiance,  from  BMW, 
second  maximum  is  0.2  x  W/3*  or  6.7  $.  Using 
1.11  X  1012  cals/KT,  the  radiant  energy  lost  for 
is  1.5  X  lO^2  cals. 


The  energy  content  of  the  blast  wave  as  it  separates  from 
the  fireball  should  consist  of  the  internal  energy  content  of  the 
air  molecules  expressed  as  a  temperature  rise,  the  compression 
energy  (i.e.,  the  dynamic  pressure  of  the  pulse),  and  a'potential 
energy  or  work  terra  expressed  in  terms  of  outward  velocity  of  the 
wave.  The  measured  fireball  temperature  decrease  prior  to  the 
second  maximum  (see  EHW,  p.  69)  will  be  taken  as  the  decrease  of 
the  air  temperature  within  the  blast  wave  itself.  The  temperature 
curve  in  this  period  of  time,  expressed  as  an  increase  in  the  air 
temperature  from  ambient  can  be  represented  by 


At  =  304t“°**23 


(91) 


for  t  in  seconds  and  T  in  °K.  For  a  breakaway  time  (i.e.  the  time 
the  blast  wave  leaves  the  fireball)  of  0.015  Bees, AT  is  1800°. 

For  this  temperature  rise, A E  for  air  is  10,860  cals/mole. -*-5  The 
peak  overpressures  from  EKW,  p.  109,  can  be  represented,  in  part, 
for  a  20  KT  air  burst  by 


p  =  1.20  X  lO^r"2^1,  r  -  2.4  X  101*  to  5-0  X  10^  (92) 


for  r  in  cm  and  p  in  psi  overpressure .  The  variation  of  the  fire¬ 
ball  radius  with  time  up  to  the  breakaway  from  EKW,  p.  65,  can  be 
represented  by 


r20  "  loW,  t  =  10-*  to  0.015  secs  (93) 

for  r  in  cm  and  t  in  secs.  At  0.015  sec,  r  is  1.04  X  10^  cm,  and, 
extrapolating  Eq.  92  to  1.04x  10*  cm,  p  is  620  psi  (42  atmos). 

The  corresponding  value  of  the  peak  dynamic  pressures  is  1330  psi 
which  is  equal  to  2.19  cal/cm3  (l  psi  »  1.65  X  10”3  cal/cm3) .  For 
a  temperature  of  2100°K  (l800  +  300),  the  number  of  moles  of  gas 
atoms  per'  cm3  at  a  pressure  of  43  atmospheres,  using  the  perfect 
gas  law,  is  2.4  X  10"*.  The  energy  absorbed  by  the  gas  due  to  the 
compression  is  then  9*1  X  10 3  cals/mole.  The  outward  velocity  of 
the  blast  wave  is  obtained  by  differentiating  Eq.  93  to  give 


(94) 


1.85  x 


icAt"0*^ 


Using  the  gross  kinetic  energy,  l/2  rav^,  as  measure  of  the  potential 
work  of  the  wave  in  moving  the  air  molecules  outward  and  converting 
to  number  of  moles  of  air  and  to  caloric  units,  the  energy  content 
at  0.015  secs  is 


Qj.  =  2.28  X.  lO^ng 


(95) 


where  ng  is  number  of  moles  of  air  in  the  blast  wave  at  the  breakaway. 
The  total  energy  in  the  blast  wave  is 


02  =  1.09  X  lO^tig  +  0.91  X  10^  +  2.28  X;  lO^ng 
=  4.28  X  10^n2 


(96) 


The  sum  of  and  Qg  should  be  equal  to  0-78  X  2.22  X  10^3,  or  1.73 
X  10^-3,  cals  for  the  20  KT  air  burst.  Thus, 

1.73  X  1043  b  2.06  X  10^  +  6.1  X  1011  +  4.28  X  10^  (97) 


In  order  to  solve  Eq.  97  for  m  and  n g,  a  relation  between  and  ng 
was  required.  Since  none  was  directly  available,  it  was  assumed 
that  all  the  gas  molecules  in  the  fireball  at  the  second  maximum 
and  in  the  blast  wave  at  breakaway  originally  occupied  the  volume 
enveloped  by  the  fireball  at  the  breakaway.  That  is  to  say,  the  air 
molecules  outside  the  blast  wave  do  not  absorb  any  of  the  Qj_  or  Qg 
energies  until  the  blast  wave  hits  them.  With  this  assumption,  the 
perfect  gas  law  estimate  is  that  the  volume  originally  contained 
1.94  X  10$  moles  air  and 

+  ng  b  1.94  X  108  (98) 

With  Eq.  9$,  the  value  of  n-,  is  5*15  X  10?  moles  and  ng  is  1.42  x  108 
moles.  The  corresponding  value  of  Qj.  is  1.12  X  10^-3  cals  or  0.51  of 
the  total  energy;  the  value  of  Qg  is  6.08  x  1012  cals  which  is  0.27  of 
the  total  energy.  The  latter  value  for  the  fraction  of  the  energy, 
carried  by  the  blast  (and  shock)  wave  is  about  half  of  that  given  in 
EHW. 
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For  a  surface  detonation,  about  the  same  energy  partition 
should  occur  except  that  the  fraction  of  the  energy  lost  from  the 
fireball  up  to  the  second  temperature  maximum  by  radiation  should 
be  less.  If,  for  a  given  yield,  the  variation  of  the  temperature 
and  fireball  radius  with  time  before  the  second  maximum  for  the 
surface  burst  is  similar  to  that  for  the  air  burst,  then  the  ratio 
of  radiant  energy  lost  for  the  two  types  of  bursts  should  be  equal 
to  the  ratio  of  the  products  of  the  fourth  power  of  the  tempera¬ 
ture^)  at  the  second  maximum  and  the  surface  area  of  the  fireball 
at  that  time.  Taking  a  spherical  fireball  for  the  air  burst  and 
a  hemispherical  shape  for  the  surface  burst,  the  ratio  of  the  two 
energies  radiated  into  the  air  space  around  the  two  is 

t|(b)  B|(s) 

q(s)/<i(a)  =  — r - ^ -  (99) 

2T^(a)  Bg(a) 

where  To  and  Rg  are  the  temperature  and  fireball  radius  at  the  sec¬ 
ond  maximum,  respectively,  and  (s)  and  (a)  refer  respectively  to  the 
surface  or  air  burst.  If  it  is  further  assumed  that  the  volume  of 
the  fireball  at  the  second  maximum  is  the  same  for  both  types  of 
detonations  (i.e.  and  contains  approximately  the  same  fraction  of 
the  total  energy) ,  then 


(100) 

For  this  ratio  of  radii,  Eq.  99  is 

r 

<l(s)/q(a)  =  0.79  |T2(s)/T2(a)  J  (101) 

The  probable  maximum  value  of  Tg(s)/Tg(a)  is  one;  then  for  a  q(a) 
value  of  0.06TW,  q(s)  is  0.053W.  A  probable  minimum  value  of 
Tg(s)/Tg(a)  might  be  0.5;  the  corresponding  value  of  q(s)  is  O.OO33W. 
The  fraction  of  the  total  energy  in  the  fireball  at  the  second  maxi¬ 
mum  for  the  surface  burst  should  therefore  be  between  0.53W  and 
0.5&7.  These  calculations  assume  that  the  radiant  energy  that 
strikes  the  ground  surface  has  been  utilized  to  vaporize  soil 
whose  vapors  are  contained  in  the  fireball  volume.  Some  of  this 
energy  might  be  lost  in  heating  soil  that  may  enter  at  later  times.  - 
But  on  the  other  hand,  the  shock  wave  would  contribute  some  energy 
in  heating  soil  materials  whose  gaseous  (or  liquid)  products  cannot 


escape  the  fireball  volume  by  outward  movement  as  in  the  air  burst. 
A  mean  value ,  0.55 W,  was  therefore  do  IrUe  go  UjjSa ug  v/± 

energy  contained  in  the  surface  burst  fireball  at  the  second  tem¬ 
perature  maximum.  In  the  calculations  to  follow,  these  energy 
partitions  were  assumed  to  be  invariant  with  yield;  this  assump¬ 
tion,  as  is  shown  in  the  following  paragraphs,  leads  to  a  parti¬ 
cular  set  of  solutions  to  the  derived  scaling  equations  for  the 
fireball  temperature  and  vapor  concentration. 

4.2.3  Fireball  Temperature,  Volume,  and  Energy  Content 
Sealing  Functions  for  Model  Air  Burst 


The  scaling  equations  for  the  air  burst  are  presented  here 
as  a  basis  for  deriving  comparable  functions  for  the  surface  burst 
in  the  following  subsection.  Geometry,  energy  content,  and  assump¬ 
tions  about  the  differences  and/or  similarities  in  the  two  types  of 
bursts  have  been  used  to  alter  or  retain  the  necessary  scaling 
functions.  Also,  some  of  the  scaling  functions  in  the  EMW  have 
been  altered  to  obtain  consistent  sets  of  functions  for  the  model 
bursts . 


The  variation  of  the  fireball  radius  after  the  second  temper¬ 
ature  maximum  was  derived  from  Fig.  2.85,  p.  65  of  ENWj  for  a  20-KT  yield, 
the  values  of  620  ft  at  0.15  sec  (t2 ,  time  Of  the  _ second  maximum)  and 
76o  ft  at  0.75  sec  (t_,  time  of  maximum  expansion  of  fireball)  were 
used.  These  values  give  the  variation  of  the  20  KT  yield  fireball 
radius  with  time  as 


R(20  KT)  =  1.89  x  10^(t/t2)0‘128 


(102) 


for  R  in  cm  and  t  in  secs .  For  the  yield  scaling  of  the  radius , 
the  geometrical  scaling  function,  wV3,  was  used  instead  of  Ww  5 
as  given  in  ENW1  The  higher  power  of  W  given  in  BRW  may  have  been 
used  to  account  for  the  decrease  in  external  air  pressure  for  the 
required  air  burst  height.  For  the  model  air  burst,  however,  the 
external  pressure  is  assumed  to  be  constant  at  1  atmosphere.  Using 
the  20-KT  fireball  radius  at  the  second  maximum  to  solve  for  the 
coefficient  of  W-*-/3  (i.e.,  as  being  representative  of  the  model 
burst  conditions),  Eq.  102  becomes 


R  =  6.96  x  lO3^*333  (t/t2)0clS8,  t2^t^tm  (103) 
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for  W  in  EP.  If  the  time  for  full  vxpansion  of  the  fireball  is  also 
taken  to  scale  according  to  WV3,  then  (from  t-,  =  0-75  see  and  W  => 

20  KT) , 


=  0.276W0*333 


(104) 


From  ENW,  p.  331;  the  scaling  for  tg  (note  change  of  notation  from 
ENW) ,  adjusted  to  give  0.15  sec  for  20  KT,  is 

t2  =  0.0336w°‘^°°  (105) 

The  ratio  t  /to  is 
m'  d- 

t  /t0  ^  8.21W“°‘16T  (106) 

m  t- 

By  eq.  106,  tg  and  tffl  approach  each  other  with  increasing  yield  and, 
if  the  relation  held  for  all  yields,  the  time  of  maximum  expansion 
and  of  second  temperature  maximum  would  coincide  for  a  yield  of  about 
3  X  105  KT  (300  MT)  .  Substitution  of  Eq.  106  in  Eq'.”  103  gives,  for 
the  fireball  radius  at  maximum  expansion. 


R  «  0.12  X  lO3#0’312 
m 


(107) 


for  Rm  in  cm. 

Scaling  relationships  for  the  fireball  temperature  can  be  der¬ 
ived  from  the  information  given  by  ENW,  pp.  300-333*  The  thermal 
power,  P,  is  given  by 


P  -  1.73  X  10-11T^R2  cals/  sec 


(108) 


for  T  in  °K  and  R  in  cm.  The  radiant  energy  lost  up  to  tg  is  given  by 


(109) 
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TABLE  2 


Summary  of  Graphical  Integration  of  Scaled  Thermal 
Power  Function  up  to  the  Time  of  the  Second  Maximum 


y 

rj 

Li 

o.iAz 

J  Z  dy 
*'0 

0.1 

0.05 

0.0022 

0.002 

0.2 

0.12 

0.0085 

0.011 

0.3 

0.24 

0.0180 

0.029 

0 .  4 

0 . 43 

0.0322 

0.061 

0.5 

O.67 

0.0555 

0.116 

0.5 

0.82 

0.0755 

0.192 

0.7 

0,91 

0.0866 

0.278 

0.8 

O.96 

0.0935 

0.372 

0.9 

0.99 

0.0374 

0.469 

1.0 

1.00 

0.0995 

O.569 

Graphical  integration  of  the  first  tern:  of  Eq,  116  gives  O.356;  the 
complete  integral  up  to  y  values  larger  than  1.4  is  then  given  hy 


=  P2t2  jo.356  +  2.17  [(1.4)“°' 


59  _  y-0.59| 


(117a) 


=  Pp-to  2.135  -  2.l7y“°’59| 


According  to  Eq.  108,  the  thermal  power  at  the  second  maximum  is 


F-o  =  1.73  x  10_11Tit 


(118) 


Substitution  of  Eqc.  103^  10 5,.  and  118  into  Eq.  117  gives 

q2(a)  =  2.82  X  10 " Vf 1,167  (2.135  -  2.17y“0,59)  (119) 

At  temperatures  below  about  1000oK,  the  energy  loss  by  radiation 
should  drop  very  rapidly  and  be  negligible  for  T/Ty  less  than  0.1 
In  order  to  estimate  a  value  of  y  for  which  T/T2  is  0 . 1,  Eqs.  10 3, 

107  y  lOo,  110  and  114  ca.11  be  combined  to  give 


T1^  »  T^R|  X  1.28(t/t 2)_1*590;  t>tn 


(120) 


It  may  "be  noted  here  that,  for  the  20  KT  yields  occurs  at  t/^ 
equal  to  5-0  whereas  the  right  hand  side  of  Eq.  120  applies  for  t/tg 
values  as  low  as  1.4  which  would  imply  that,  for  the  curve  given,) 

R  occurred  at  a  tm/t2  value  of  1.4  independent  of  yield  rather  than 
at  the  value  given  by  Eq.  106  for  20  KT.  It  is  likely  that  the  P/P2 
curve)  at  least  up  to  tjjj/tg,  is  yield-sensitive  and  that  the  curve 
in  ENW,  if  an  observed  curve,  is  for  a  yield  much  greater  than  20  KT 
At  any  rate,  for  times  near  t-m  and  longer,  the  curve  indicates  T  is 
proportional  to  (t/tg)  .  Complete  substitution  of  the  indicated 

equations  into  Eq.  120  and  solving  for  T  gives 

T  =  0.929W0,0108  T2  (t/t2r°'398>  t~tm  (121) 

For  a  large  range  of  yields  the  value  of  0.929^*^^^  is  near  one; 
hence  the  value  of  t/tp  or  y  for  T/Tg  equal  to  0.1  is  about 
325(2. 17y"0*59  =  0.072).  Using  this  value  in  Eq.  119  gives 


q2(a)  =  5.82  X  3jO“5Tg  W1*167 


(122) 


The  total  radiant  energy  is  then 

qT(a)  =  7.43  X  10” ^Tp  W1,167  (123) 

The  fraction  of  this  amount  up  to  tg  is  0.22.  Nov,  if  q^a)  is  taken 
as  1/3  W  in  KT  (or  3-70  X  lO11  W  in  Cals),  then 

T2  =  8,390w"°*°42  (124) 

Equation  124  shows  that  when  the  energy  partition  is  postulated  to 
be  independent  of  yield,  the  temperature  at  the  second  maximum  de¬ 
creases  with  yield. 

For  a  yield  of  20  KT,  Eq.  124  gives  7400°K  as  the  value  for 
Tpj  this  value  is  almost  1000°K  lower  than  shown  in  ENW,  p.  68  and 
69.  A  lower  temperature  results  even  when  the  fireball  radius  equa¬ 
tions,  p.  66,  arc  used  for  R  in  Eq.  120.  If  qT(a)  is  taken  as  O.58 
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W(3CT),  as  has  been  found  in  the  previous  section;  then  Eq.  123  would 
give  a  second  maximum  temperature  of 


9,64a/ 


-0.042 


(125) 


For  a  yield  of  20  KT,  Eq.  125  gives  8500°K  as  the  value  for  T2  which 
is  in  better  agreement  with  the  values  given  in  EHW,  p.  68  and  69. 

To  obtain  the  temperature  given  by  Eq.  125,  all  of  the  energy  in  the 
fireball  is  used  as  radiant  energy;  none  of  the  energy  has  been  allo¬ 
cated  to  further  expansion  of  the  fireball  and  atomic  cloud  or  for 
the  work  against  gravity  in  raising  the  gas  and  particle  mass  to  high 
altitudes.  However,  there  are  two  additional  sources  of  energy  that 
have  not  been  mentioned  that  could  easily  provide  this  "extra"  energy; 
these  are:  (l)  the  buoyant  force  of  the  atmosphere  to  lift  the  warm 
gas  mass,  and  (2)  the  kinetic  energy  and  the  energy  of  cooling  and 
expansion  of  the  air,  originally  at  1  atmos  and  298°K,  to  a  lower 
final  pressure  and  temperature  at  higher  altitudes;  this  is  a  net 
gain  from  the  air  molecules,  themselves,  that  were  enveloped  in  the 
fireball.  Also,  some  of  the  blast  energy  will  be  reclaimed  as  pre¬ 
heated  air  in  the  rising  cloud.  Thus,  Eq.  125  may  well  not  be  the 
maximum  value  for  Tg  with  respect  to  the  fraction  of  the  total  bomb 
energy  lost  by  radiation. 

If  qp(a)  is  taken  as  equal  to  1.11  x  1012a  W  calories  (a  being 
the  fraction  of  the  total  energy  released  which  is  lost  by  radiance), 
then  T2  from  Eq.  112  is 

t2  =  l.io  x  10V-/4  w-°-042  (126) 


Substitution  of  Eq.  126  into  Eq.  121  gives 

T  =  1.02  X  lO^GT1/11  w'°'°31  (t/t2)-0'398,  t>tm  (127) 


And  further  substitution  of  Eq.  105  for  tg  gives 

T  =  266oarL//1+  W0'168  t‘0,39y  t2t 

rn 


(128) 


In  view  of  the  possibility  of  variations  in  the  scaled  power  curve 
of  EM-/  due  to  yield,  it  may  be  assumed  that  the  lower  time  limit 
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for  application  of  Eqs.  127  and  128  is  t£  rather  than  tffl.  With  this 
assumption,  the  substitution  of  Eqs.  103  (or  10?)  and  12?  in  Eq.  108 
gives 

p  =  9.16  x  io12aw°-^3  (t/tg)"1*33^,  tg&t^  (129) 

and 

P  =  1.57  X  lO^a  W° ' 500  (t/t2)"1,590,  t^trn  (130) 


For  these  approximations  of  the  thermal  power.,  the  energy  lost  by 
radiance  after  tg  is 


qg(a)  =  9.16  X  10 


12a«°-^t2L-33h  t-1<334  at 

J%r 


J2 


+  1.57  x  lo^aw0’500  tg'590^ 


t~lA9°  at 


(131) 


m 


If  the  final  limit  of  integration,  tp,  is  again'  defined  as  the  time 
when  t/t?  is  0.1,  then  by  use  of  Eq.  126  and  Eq.  127  or  128,  tp  is 
given  by 


tf  =  9.09W°'52T 


(132) 


Integration  of  Eq.  131  up  to  this  limit  yields 
q2(a)  =  9.22  x  10ua  W1,0^3  (1  -  0.235U°'056  -  0 .0357W“° -°59)  (133) 


Thus,  even  though  the  scaling  equations  were  derived  on  the  basis  of 
q2(a)  being  proportional  to  T.7 ,  the  combination  of  the  T  and  R  func¬ 
tions  from  EKW  (with  some  alterations  for  R)  gives  qp(a)  as  not  being 
directly  proportional  to  W.  However,  the  p over  of  W  is  not  far  from 
one  so  that  q2(a)  from  Eq.  133  conforms  approximately  to  the  original 
postulation.  The  main  difference  is  that  most  of  the  yield  dependence 
in  p/?2  occurs  in  the  period  from  tp  to  tfU  when  the  fireball  is  still 
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expanding;  the  extrapolation  of  the  variation  of  T  with  t/t2  hack  to 
this  time  period  may  also  account  for  some  of  the  difference  in  the 
yield  dependence. 

4.2.4  Estimate  of  Difference  in  Fireha.il  Temperature  at 

Second  Maximum  Between  Model  Air  and  Surface  Bursts 


The  main  assumptions  that  were  made  to  obtain  an  estimate  of 
the  difference  in  the  fireball  temperature  at  the  second  maximum 
between  the  model  air  and  surface  hursts  are:  (l)  in  the  surface 
burst,  l/2  of  the  energy  in  the  fireball  is  utilized  to  heat,  dis¬ 
sociate,  and  expand  air  molecules  and  l/2  of  the  energy  tu  vapunze, 
heat,  dissociate,  and  expand  molecules  of  soil;  (2)  the  air  and  sur¬ 
face  bursts  are  equally  efficient  in  the  heating  process  with  respect 
to  the  number  of  gas  atoms  heated  to  the  temperature  of  the  second 
maximum  per  unit  of  yield,  and  (3)  the  fraction  of  the  total  yield 
utilized  at  the  second  temperature  maximum  is  0.51  for  the  air  burst 
and  0.55  for  the  surface  burst. 


The  first  law  of  thermodynamics  for- the  utilization  of  0.51 
of  the  total  energy  at  the  second  maximum  for  the  model  air  burst  is 


(13^) 


in  which  n«  is  the  number  of  air  atoms  (w  and  0  atoms),  Cv  is  the 
heat  capacity  of  (an  ideal)  monatomic  gas  (3  cal/ mole -deg) ,  Tp  is 
the  temperature  at  the  second  maximum,  T0  .is  the  original  tempera¬ 
ture  of  the  gas  atoms  (298°k),Ae°  is  the  heat  of  dissociation  of 
air  molecules  at  298°K  (79,200  cat/moleof  gas  atoms),  PQ  is  the 
external  pressure  (1  atmosphere),  Vg  .is  the  volume  of  the  fireball 
at  the  second  maximum,  VQ  is  the  original  volume  of  the  air  mole¬ 
cules  at  PQ,  V  ,  and  T0,  and  is  O.5IW  in  KT  or  5-66  x  lO^-W  in 
cal.  For  a  spherical  fireball,  the  volume  from  Eq.  103  is 


V  *  1.41  X  1012W  (t/t2)°’385  cm3  (135) 


so  that 

V2  =  1.4l  X  1012W  (136) 


Using  the  perfect  gas  law,  the  value,  of  P0VQ  for  the  original  air 
molecules  is 
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PoVo  - 


(137) 


Substitution  of  the  appropriate  terms  in  Eq.  134  (l  cal  =  4l.29  cm3 
atoms)  gives 

5.66  X  10nW  a  n23(T2-298)  +  79*200  +  3-42  X  1010W  -  300n£  (138) 


so  that 


/„  5.32  X  loll 

VW  “  f'3  To^OOOj 


(139) 


In  order  to  make  a  similar  computation  for  the  model  surface 
burst,  the  thermal  properties  of  the  "ideal"  soil  must  be  specified. 
The  nonreactive  "ideal"  soil  will  be  assigned  the  thermal  properties 
of  anorthoclase  (NapO'A^Oj’SSiOg)  to  be  consistent  with  the  selected 
melting  point  (l400^C).  One  mole  of  the  soil  upon  complete  dissocia¬ 
tion  gives  26  moles  of  gas  atoms;  the  molecular  weight  is  524.  In 
estimating  the  energy  computations,  the  bonding  energy  between  the 
three  oxides  was  assumed  to  be  small  compared  to  the  oxygen  bonding 
to  the  metal  atoms .  The  changes  in  energy  for  any  process  involving 
the  whole  molecule  therefore  raay  be  equated  to  the  sum  of  the  energy 
changes  for  the  individual  oxide  molecules  indicated  in  the  formula. 
The  entropy  of  fusion  was  taken  to  be  29  cal/mole/deg;  the  heat  of 
fusion  is  then  48,000  cal/mole. 

The  change  in  internal  energy  of  the  atoms  and  molecules  was 
estimated  from 


As  -  Ah  -A(PV) 


(140) 


and  when  gas  atoms  are  involved,  the  perfect  gas  law  was  used  and 
the  quantity AnRT  substituted  for  Z\(PV)  for  a  constant  temperature 
process.  The  dissociation  energy  of  the  soil  at  298°K  was  calculated 
as  follows  from  the  data  of  Shell  and  Sinkel5  (s  =  solid;  g  =  gas; 

ZiH  and  An  values  in  cal  for  indicated  reaction  in  number  of  moles): 


24, 500 


(1) 

2  Na(s)  +  1/2  02(g) — 

— ?  MagO  ( s  ) , 

Ah° 

(2) 

2  Na(s) — 

2Na(g), 

&H° 
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(3)  1/2  02(s)  — >o(g),  Ah0  =  59 >100 

Sum:  (3)+(2)-(l)  Na20(s) - }2Na(g)  +  0(g),  Ah°  =  144,700 

Ae°  =  144,700  -  3  X  2  X  298  *  142,900 

(1)  2  Al(s)  +  3/2  02(g) - ^Al203(s),  Ah°  =  -490,200 

(2)  2  Al(s) - >2Al(g),  Ah°  =  154,300 

(3)  3/2  02(g) - *30  (g),  Ah°  =  177 > 300 

Sum:  (3)+(2)-(l)  Al203(s)  - - »?2Al(g)  +  30(g) ,Ah  =  731,800 

Ae°  =  731,800  -  5  X  2X  298  =  728,800 

(1)  Si(s)  +  02(g)  - )  Si02(s),  Ah°  =  -210,200 

(2)  Si(c)  - >Si(g),  Ah°  «  90,000 

(3)  02(g)  - >20  (g),  Ah°=  118,200 

Sum:  (3)+(2)-(l)  Si02(s)  - - >  Si(g)  +  20(g) ,  Ah°  =  4l8,400 

Ae°  =  731,800  -  3  X  2  X  298  =  416,600 


For  the  process , 

Ka20*Al202’6Si02(s) — »2Na(g)  +  2  Al(g)  +  6Si(g)  +  l6  0(g)  (l4la) 

Ae°  *  142,900  +  731,800  +  (6  x  4i8, 4oo) 

=  3,370,000  cal/mole  of  soil  (l4lb) 

or, 

Ae°  =  129,600  cal/mole  of  gas  atoms  formed  (l4lc) 


The  first  law  of  thermodynamics  for  the  utilization  of  0.55 
of  the  total  energy  at  the  second  maximum  for  the  model  surface  hurst 
is 
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(142) 


n2  [cv(T^“T0)  +  72,200]  -  (n|/2)RTQ 

+  4  [cv(t£-tq)  +  129,600]  +  3.42  X  1010W 

t  ” 

in  which  n2  is  the  number  of  moles  of  air  atoms  in  the  fireball, 
is  the  number  of  soil  atoms  in  the  fireball,  Tl  is  the  temperature 
of  the  gas  mixture  at  the  second  maximum,  and  is  O.99W  in  KT  or 
6.10  X  lQ-^W  in  cals.  Making  the  substitutions  indicated  and  let¬ 
ting  the  sum  of  n^  plus  be  represented  by  n  gives,  for  Eq.  142, 


5.76  X  101XW  =  n(3  T^+78,000)  +  91 , 600ng  (143) 

The  model  surface  burst  was  defined  as  one  in  which  half  the 
energy  in  the  fireball  at  the  second  maximum  was  utilized  in  dis¬ 
sociating  and  heating  air  molecules  and  half  for  the  soil  molecules. 

If  each  of  the  two  component  gases  are  expanded  from  their  original 
volumes  to  Vg  (same  as  for  the  air  burst),  then 

n2  [gv(T2-T0)  +  79,200]  +  [pqV2  -  (n»/2)RTo] 

=  n2  [Cv('^-Tc)  +  129,600]  +  PCV2 

Substitution  of  (n  -  nip  for  rig  and  solving  for  n^  gives 

(30^78,000) 
n2  “  n  (6T^+207,600) 

With  Eq.  14-5,  Eq.  14-3  becomes 

'  .  91,600 

X  6t^+207,600) 

The  values  of  iv^/W  for  Eq.  139  and  n/w  for  Eq.  146  are  summar¬ 
ized  in  Table  3  and  the  values  of  n/w  is  plotted  as  a  function  of  Tg 
in  Fig.  2.  The  values  of  Tg  and  Tl  for  equal  values  of  rvj/w  and  n/w 
are  given  in  Table  4.  The  assumption  of  equal  values  of  the  "vapori¬ 
zation  efficiency"  (n/w)  for  the  two  types  of  burst  gives,  approximately, 
a  constant  difference  in  the  temperatures  at  the  second  maximum  of  about 
3800°K. 


9.76  X  lO^W  =  n  (3T£+78,000) 


(144) 


(14-9) 


(146) 
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TABLE  3 


Values  of  n2/w  and.  n/w  for  Several  Values 
of  T2  and  Tg,  Respectively 


ro3  I 
"o°  1 

rii 

.  :  n2  h 

(l0omoles/RT) 

:  g/W 

(lpP  ffioles/KT 

.4,000 

5-91 

.  5.24 

.  e  nAn 

J 

c.  rro 
j  *  i  *- 

5-09 

6,000 

5.54 

4.95 

7,000 

5-37 

4.82 

8,000 

5.21 

4.70 

9,000  _ 

5-Q6 

4.58 

.  10,000 

4.93 

4.47 

11,000' 

4.79 

- 

12,000 

.  4.67 

“ 

TABLE  4 


Values  of  T2  and  Tg  for  Equal  Values  of  rig/w  and  n/W 


n/w 

or  n?/w 
(moles/KT) 

t2 

(°K) 

T2 

(°k) 

t2-t'  ; 

(°K) 

T'/t  0 

2  2 

4.7 

11,750 

8,000 

3,750 

0.681 

4.8 

10,910 

7,150 

3,760 

0.658 

4.9 

10, 150 

6,370 

3,780 

0.628 

5-0 

9,420 

5,640 

3,780 

0.599 

5-1 

0,730 

!4930 

3,800 

O.565 

5.2 

8,050 

4,250 

3,800 

0.528 

Av. : 

3,780 
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SECOND  MAXIMUM  TEMPERATURE  (°K) 


Fig-  2  Variation  of  the  Moles  Gas  per  KT  in  the  Fireball  With 
the  Second  Maximum  Temperature 
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If  the  value  of  T2  for  the  model  air  burst  for  a  yield  of  20 
KT  is  8300°K  as  given  in  EI~!W ,  p.  68  and  69,  then  T g  is  4500°K.  This 
value  of  Tg  seems  to  he  too  low  to  he  consistent  with  the  energy 
utilization  and  energy  partition.  The  temperature  of  4500°K  is  not 
high  enough  to  completely  dissociate  ail  the  air  molecules  and  many 
of  the  more  stable  oxide ’molecules .  If  the  dissociation  energies 
in  Eq.  142  were  reduced,  a  larger  value  of  n/w  at  a  given  value  of 
would  result.  This,  in  turn,  would  give  a  smaller  difference 
between  Tg  and  Tg;  and,  if  the  original  value  of  Tg  is  retained 
a  higher  value  of  TA  would  be  indicated.  A  higher  value  of  Tg  would 
require  use  of  the.  dissociation  energies  (i.e.  Eq.  142)  and  a  return 
to  the  above  temperature  values.  This  inconsistency  can  be  removed 
if  Tg  is  increased. 

The  data  from  ENW  used  in  the  previous  sections  to  obtain  the 
various  scaling  functions  was  associated  with  the  model  air  hurst. 
However,  as  was  noted,  EHW  does  not  state  the  type  of  burst;  except 
for  the  blast  curves,  p.  109,  no  mention  of  the  burst  type  is  given. 
Many  of  the  inconsistencies  found  in  the  ENW  information  with  respect 
to  the  fireball  temperature,  energy  jpCLx,"bl  "b  ion,  and  thermal  power  would 
he  resolved  if  it  were  assumed  that  the  data  actually  applied  to  a 
tower-mounted  explosion.  If  this  assumption  is  applied  and  if  it  is 
also  assumed  that  the  temperature  at  the  second  maximum  for  a  tower 
shot  is  midway  between  that  for  a  surface  burst  and  an  air  burst, 
then  the  latter  two  temperatures  can  be  estimated  from 


and 


Tg  +  T£  =  16,600  (147) 

Tg  -  Tg  =  3,600-  (148) 


The  resulting  values  of  T2  and  TA  are  10,200  °K  and  6400  °K,  respec¬ 
tively. 


Tiie  assumed  scaling  equations  for  R  and  tg  and  the  assumption 
that  the  radiant  energy  ms  proportional  to  the  yield  led  to  the 
result  that. Tg  was  proportional  to  w “0*042  Sxnce  qm(a)  was  propor¬ 
tional  to  Tg  w^*167  when  the  yield  dependence  of  the  term  R^tg  were 
substituted.  If,  however,  instead  of  assigning  the  excess  over 
1.000  in  the  exponent  of  W  to  T2  alone,  the  excess,  0.l6Y,  is 
divided  equally  and  assigned  to  the  scale  functions  for  Tg  and  tg, 
and,  if  qg(a)  is  adjusted  to  0.5IW  (in  KT)  for  the  20  KT  yield 
model  air  burst,  the  following  scaling  functions  result: 
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R  =  5*74  X  lo3w°*333  (t/t2)^*'''28  cm,  tg™t— tm 

(149a) 

R  3  8.58  X  lo3w°*280  t0,128  cm,  t2-s.fs.trfl 

(149b) 

t2  =  0.0431*  W0,1*16  sec 

(150) 

tm=  0.286***  w0-321**  sec 

(151) 

Rm  =  7.32  x  103w°-321**  cm 

(152) 

T  =  1.02  X  10  V0*021  Ct./-fc2)“0*^00  °K,  t2tg 

(153a) 

T  =  2.90  X  lO^H0,1^  t-0.4O0  0^  fgztg 

(153b) 

T2  »  1.02  X  10  V0*021  °K 

(154) 

tf  =  13.63W0,lj;l6  sec,  for  T  =  0.1  t2 

(155) 

P  =  9-06  X  1010W1,li+3  t'1’3^4  cals/sec,  t2£t5tm 

(156) 

P  =  6.59  X  lO10!*/1*22'’  t"1,600  cals/sec,  t^tS-^' 

(157) 

qq(a)  =  7*77  X  1010W  cals  (7  #  of  total) 

(158) 

02(a)  =  7.76  X  10i;LW(l  -  0.222W0’°31  -  0.0295W"0,025) 

(159) 

For  Eq.  159;  q2(a)A>r  in  same  units)  is  0.52. for  1  HI; 

0.50  for  100  Id,  0.49  for  1  MT,  and  0.48  for  10  MT.  This  relative 
decrease  with  yield  in  the  loss  of  fireball  energy  by  radiance,  due 
to  the  second  tern  in  the  pare;  nt  he  pis  of  Sq,/159;  is  a  result  of 
the  scaling  function  for  tm.  If  the  t,,,  scaling  function  exponent 
on  W  were  the  same  as  that  "for  tg,  this  term  in  Eq.  159  would  not 
be  yield-dependent.  The  indicated  drift  of  t  towards  tg  with  in¬ 
creasing  yield  also  results  in  a  decrease  in  the  difference  between 
Rg  and  Rm  with  yield  and  hence  less  relative  increase  in  the  surface 
area  of  the  fireball  from  which  the  energy  is  radiated. 

Besides  the  difference  in  temperature  at  the  second  maximum, 
the  main  difference  in  the  above  scaling  functions  from  those  given 
in  the  proceeding  paragraph  is  that  due  to  the  18  $  decrease  in  the 
fireball  radius.  This  decrease  resulted  from  the  adjustment  of  dg(a) 
to  0.5IW  (in  KT)  at  20  KT.  This  "used  up"  all  the  energy  so  that  the 


^Adjusted  to  0.15  secs  for  20  KT. 

**Taken  to  be  equal;  determined  from  Eq.  121b. 
***Ad justed  to  0.75  secs  for  20  Id’. 
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fireball  radius  thus  determined  is  a  maximum  value.  Lapple1^  indi¬ 
cates  that  the  fireball  radius  may  actually  be  only  about  0-5  of 
the  values  given  by  EEW.  The  decrease  in  the  radius  does  not  result 
in  a  significant  change  in  difference  in  temperature  at  the  second 
maximum  between  the  model  air  burst  and  model  surface  burst;  the 
difference  is  decreased  by  only  about  100°K.  The  reestiuate  of  Tp 
is  therefore  6500°K. 

1.2.5  Scaling  Functions  for  Temperature-,)  •  Fireball  Radius , 
and  Utilization  of  Fireball  Energy  for  Model  Surface 
Burst 


The  scaling  functions  for  the  model  surface  burst  were  based 
on  those  for  the  model  air  burst  given  in  the  previous  subsection 
and  altered  to  account  for  assumed  differences  between  the  two  types 
of  bursts . 

Assuming  equal  fireball  volumes  for  the  two  and  a  hemispherical 
shape  for  the  surface  burst  fireball  up  to  maximum  expansion  results 
in  a  radius  for  the  latter  that  is  2V 3  times  that  for  the  air  burst; 
or 


R  =  7.23  X  10^.7° '333  (t/t2) 0,128  cm,  t2-tS:tia  (160) 

The  dependence  of  the  time  of  the  second  maximum  on  yield  for  the 
model  surface  hurst  will  be  assumed  to  be  intermediate  between  that 
given  by  BIIW  and  that  given  by  Eq.  150  for  the  model  air  burst;  at 
1  KT  the  times  at  which  the  second  maximum  occur  will  be  taken  to 
be  equal  for  the  two  types  of  burst.  In  other  words,  tla  change(s) 
in  temperature  with  time  for  the  surface  burst  is  assumed  to  be 
just  a  little  slower  for  the  surface  burst  than  for  the  air  burst. 
The  value  of  t2,  according  to  the  above  assumptions,  is  then 

tg  =  0.0431W0,460  secs  (l6l) 

Substitution  of  Xq.  lol  into  Eq.  l60  gives 

R  =  1.0b X  loH/0*2^  t0,128  cm,  t2  —  t  — t ..  (162) 

By  again  letting  the  thl  and  R.a  dependence  yield  have  the  same  expon¬ 
ent  on  V  and,  if  the  value  of  ti:,  for  1  KT  is  taken  to  be  the  same 
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as  for  the  model  air  burst  (F,q.  15l) ,  the  time  at  which  maximum 
expansion  of  the  fireball  occurs  is  given  by 


tm  =  0.286w^‘^^  sec 


(163) 


Substitution  of  Eq.  163  in  Eq.  162  gives 


p  _  a .  9?  y  i  o3w°  •  3l4 
in  "  —  " 


cm 


(164) 


The  temperature  at  the  second  maximum  for  the  surface  burst 
will  be  assumed  to  be  independent  of  yield  so  that,  from  the  pre¬ 
vious  subsection 


t2  =  6,500  °IC 


(165) 


For  times  just  beyond  t2,  the  decrease  in  temperature  will  be  assumed 
to  be  proportional  to  (-fc/t^) ^  rather  than  (t/t2)“0-  100  •  This 
initial  slower  cooling  was  taken  because  of  thermal  reflection  and 
possible  entry  of  preheated  materials  from  the  earth’s  surface  into 
the  fireball  after  t2-  The  slower  cooling  curve  was  assumed  to  apply 
until  the  fireball  cooled  to  2500  °K  at  which  temperature  the  cooling 
rate  is  increased  by  use  of  Newton's  law  of  cooling  to  approximate 
the  shape  of  the  fireball  cooling  curve  at  later  times.  The  fire¬ 
ball  temperature  is  then 


*)  j  ') 

T  =  6,500(t/t2)‘‘  °IC,  t2-  l* 


(1 66) 


or,  with  Eq.  l6l 


T  =  2,23a/0,153  t”1/3  °K,  t^tSt. 


(167) 


where  t..  is  the  time  at  which  T  is  25OO  °K.  Substitution  of  25OO 
for  T  in  Eq.  167  gives 


K  =  0.759V/0 ‘h6°  sec 


(163) 
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For  times  longer  than  tx,  where  Newton's  law  of  cooling  is  used, 

T  =  Tq  +  (2500-To)e“K^t_tx)  (l 69) 

in  which  T0  is  the  external  temperature  and. K  is  the  cooling  constant. 
The  value  of  T0  will  he  taken  to  he  273  °K“*(0°C)  and  the  value  of  K 
can  he  determined  from  the  differential  of  Eq.  167  evaluated  at  tx 
and  25OO  °K  so  that  Eqs.  167  and  169  join  smoothly  at  that  tempera¬ 
ture  .  Namely 

K  =  sec-1  at  t  =  t.  (170) 

2,227 

When  the  appropriate  substitutions  from  Eqs.  167  and  168  are  made 
in  Eq.  169,  Eq.  170  is 

:  K  =  0.50ftr0*4^°  sec"1  '  (171) 


Use  of  Eqs.  168,  169,  and  171  to  determine  the  time  at  which  the  fire¬ 
ball  temperature  is  l400°C  (1673  °k)  gives 

t  (l4oo°C)  =  1.69W0,46°sec  (  (172) 

For  t  (l400°c)  equal  to  60  sec,  W  is  2.3  PIT .  If  Eq.  167  were  assumed 

to  apply  up  to  60  secs,  the  yield  would  he  about  1.0  MT*  Of  the  two 

yields,  the  larger  value  associated  with  the  more  rapid  decrease  in 
temperature  below  2500  °K  is  probably  the  better  of  the  two  estimates. 

The  thermal  power  and  energy  lost  from  the  fireball  for  the 
model  surface  burst,  for  times  after  the  second  thermal  maximum,  is 
calculated  for  energy  transmission  from  a  hemispherical  fireball  up 
to  the  time  of  maximum  expansion.  At  about  this  time  or  very  shortly 

thereafter,  the  fireball  should  be  breaking  clear  of  the  ground  sur¬ 

face  in  its  upward  rise  (see  ENW,  p.  23,  where  for  1  MT  the  top  of 
the  cloud  is  at  10,000  ft  in  about  16.5  sec,  at  20,000  ft  in  40  sec, 
and  at  25,000  ft  in  60  sec,  etc.).  For  calculating  the  distribution 
of  the  thermal  energy  during  the  fireball  rise  (until  it  cools  to 
l400°C),  it  is  more  convenient  to  "remove"  the  ground  surface  and 
allow  the  fireball  to  remain  at  constant  volume  with  a  constant  ex¬ 
ternal  pressure.  In  this  procedure,  the  shape  of  the  fireball 
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changes  from  hemispherical  to  spherical  with  a  decrease  in  radius. 
The  fireball  volume,  however,  is  assumed  to  remain  constant  until 
the  temperature  falls  "below  500  °K  or  at  least  until  the  fireball 
separates  from  the  surface.  The  spherical  fireball  radius  at  sep¬ 
aration  time,  ts,  with  the  constant  volume  assumption,  is 

Rs  =  7.32  X  lo3w°*3lijh  cm,  titB  (173) 


The  height  of  tine  top  of  the  fireball  at  this  time  would  then  be 
j.4. 0*+  a  10-JW  cm  ana  since  in  was  y.d.c 1  a  xu-'W’' '  cm  an  t^, 

its  height  is  increased  by  5-42  X  103^P  • 314  cm  between  tm  and  ts. 


Using  the  value  of  Rm  as  the  height  of  the  cloud  at  tm  and 
the  data  of  ENW,  Fig.  2.12,  for  the  ,1-MT  burst,  the  height  of  the 
cloud  top  can  be  represented,  with  a  fair  degree  of  precision,  by 


and 


\  =.  4.24X  101*  t°‘TOl+  cm,  i^sSt  ^60  sec  (174a) 

h^.  =  6.46  X  10^  ^.600  cm,  go  secr£.t^390  sec  (174b) 


Since  hj.  is  equal  to  2Rg  at  ts,  the  use  of  Eq.  173  can  be  made  to 
find  ts  for  the  1-MT  yield.  It  is  4.82  sec.  Assuming  the  equation 
form  atn  for  the  dependence  of  R  on  t  between  and  ts  and  solving 
by  use  of  F,qs.  164  and  173  for  a  1-MT  yield,  gives  a  value  of  n  of 
-0.354.  Use  of  the  general  equations  for  tm  and  Rm  with  this  power 
of  t  gives,  for  the  fireball  radius, 

R  =  5.92  X  io3#-li'25  t-0.354  cm^  tmrs.t^ts  (175) 


The  time  that  the  fireball  separates  from  ground  surface  is  then 


ts  »  0.550W0-311*  sec  (176) 

The  thermal  power  equations  for  the  time  limitations  on  the  temper¬ 
ature  and  radius  functions  are  four  in  number.  For  the  hemispherical 
shape,  the  energy  radiated  between  t2  and  tm  from  the  portion  of  the 
fireball  in  contact  with  the  ground  is  lost  from  the  fireball.  The 
four  functions  are  given  by 
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P  e  4.11  X  lO1^1*161  t-1,0''7  cals/sec,  t2^t^tm 

(177a) 

P  =  2.l8  X  io^°W-L-3^  t"^'585  cals/sec,  tmi£t~ts 

(177b) 

P  =  2.51  X  1010W1,241  if’1  * 333  cals/sec,  ts^.t^-tx 

(177c) 

and 

P  =  5.16  X  io^w0,628  [1  +  32.6e“I^t“i^  +  399e"2K^t"tx) 

+  2170e"3K^t_tx^  +  4430e“4K(t“tx)]  cals/sec,  t^tx 

(177d) 

Eq.  17Tb  was  determined,  by  joining  Eq.  177a  at  to  Eq.  177c 
at  ts  with  a  straight  line  in  a  logarithmic  coordinate  system,  as  an 
approximation  of  the  change  in  radiant  energy  emission  occurring  with 
the  change  both  in  temperature  and  surface  area  while  the  shape  of  the 
fireball  changes  from  hemispherical  to  spherical  without  change  in 
total  volume  as  it  rises  from  the  ground. 

The  radiant  energy  lost  from  the  fireball,  from  integration 
of  Pdt  for  Eq.  176,  is 

q2(a)  =  6.81  x  l01;w1,126(l  -  0.693W0,010  -  0.121W"0,038) 

+  5.16  X  10  V' 628  jt  +  W°-46°(k)60  -  64.2e"K(t"tx 
“  393e-2K(t~tx}  -  l420e"3K(t_t^  -  2l80e“^t~tx^]  cals  (178) 

For  T  equal  to  1673  °K  (l400°c),  substitution  of  the  values  for  K 
and  (t~tx)  into  Eq.  178  gives 

q2(a)  =  6.81  X  1013W1,126(1  -  0.693W0'010  -  0.0971W0'038)  (179) 


For  a  yield  of  2.3  MT,  Eq.  179  gives  a  value  of  qg(a)  that  is  0.29 
of  the  total  weapon  yield.  It  may  be  noted  that  Eqs.  177  through 
179  do  not  support  the  postulate  that  requires  qgCa)  to  be  directly 
proportional  to  W.  Since  generalized  yield  limits  were  used  to  solve 
for  (32(a),  Eq*  179  is  valid  for  a  range  of  yields  for  the  radiant 
energy  lost  between  X,£  and  the  time  at  which  the  fireball  tempera¬ 
ture  is  l400°C . 
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4.2.6  Estimation  of  n(/,)/V  for  Model  Surface  Burst 


The  estimate  of  n(£)/v  for  the  model  surface  hurst  was  made 
for  the  soil  melting  point  of  l400°C .  The  time  after  detonation  of 
60  sec  was  used  to  minimize ,  to  a  degree ,  errors  in  the  relative 
abundance  of  the  short-lived  fission  product  nuclide  composition 
(many  of  the  half-lives  are  estimated)  and  so  that  computations 
would  apply  to  a  surface  detonation  in  the  megaton  yield  range. 

The  discussion  in  the  previous  subsections  was  concerned  with 
establishing  the  energy  content  and  gas  composition  in  the  fireball 
at  the  time  of  the  second  temperature  maximum.  In  order  to  estimate 
the  molar  concentration  of  the  soil  in  the  fireball  at  the  time  and 
temperature  stated  above,  estimates  of  the  energy  utilization  between 
the  time  of  the  second  maximum  and  60  secs  are  required.  In  this 
period,  the  following  processes  were  considered:  (1)  cooling  of  the 
gases  from  T£  to  l400°C  and  liquifaction  of  the  vaporized  soil  with 
a  resultant  release  of  energy,  (2)  association  of  the  gas  atoms  with 
a  resultant  release  of  energy,  ( 3)  loss  of  energy  due  to  radiation 
from  the  fireball,  (4)  loss  of  energy  due  to  expansion  of  the  fire¬ 
ball  gases  against  the  atmosphere,  (5)  loss  of  energy  due  to  heating 
and  expansion  of  additional  air,  and  (6)  loss  of  energy  due  to  heating 
and  melting  of  a  given  amount  of  additional  soil.  These  losses  and 
gains  in  the  energy  among  the  fireball  constituents  and  external  sub¬ 
stances  have  to  be  estimated  and  balanced  in  order  to  estimate  the 
amount  of  soil  present  in  the  fireball  that  can  be  melted  at  1400°C 
in  process  (6) . 

For  a  temperature  of  6500  °K  at 
gives,  for  the  surface  burst 

n  »  4.88  X  10  moles 

The  value  of  the  ratio  of  np(gas  atoms  from  the  soil)  to  n( total  gas 
atoms)  from  Eq.  l4p  is  0-395  so  that 

=  2.95  X  10 moles  of  air  gas  atoms  (l8l) 

and 

np  =  1-93  X  10 moles  of  soil  gas  atoms  (182) 

Upon  complete  condensation  of  the  soil  vaporized  at  tg,  the  number 
of  moles  of  soil  formed  would  be : 26  atoms  per  molecule 


the  second  maximum,  Fig.  2 
of  gas  atoms  (l80) 


n"(l)  =  7.42  X  io\  moles  soil 


I 


(183) 


The  computational  procedure  for  estimating  the  energy  releases 
and  losses  for  the  first  four  processes  is  given  by  the  following 
steps:  (l)  cool  the  gases  from  65OO  °K  to  1673  °K  at  constant  volume, 
(2)  recombine  the  air  gas  atoms  to  gas  molecules  and  the  soil  gas 
atoms  to  liquid  soil,  {3)  subtract  the  radiant  energy  emitted,  and 
(4)  expand  the  gas  volume  from  V2  to  Vm  against  the  external  pressure 
of  1  atm.  The  gases  will  be  assumed  to  behave  like  an  ideal  mona¬ 
tomic  gas  so  that  the  energy  released  in  cooling  to  1673°  is,  using 
3  cal  mole“!  deg"l  for  the  heat  capacity, 

-Q^  =  7*07  X  1010W  cals  (184) 

For  the  soil  composition  of  N&20 'Al^Cu ‘6Si0g  (as  given  previ¬ 
ously)  ,  the  single  oxide  energies  were  calculated  from  Reference  15 
data  as  follows: 


(1)  Na20(s), 

(2)  Na20(s), 

(3)  0(g), 

(4)  2Na(g), 

(5)  Na20(s), 


298°K— >  2Na(g),  298°K  +  0(g),  298°K,  Ah  =  144,700 


2y8°K->  NagO^),  1673% 

Ah  =  24,750 

298°K — >0(g) ,  167 3°K  ■ 

Ah  =  6,880 

298°K — »  2Ka(g),  1673°K, 

Ah  «  13,750 

1673  -^HagOtf),  167 3°K, 

Ah  =  7,200 

Sum:  (l)-(2)+(3)+(4)r(5);  Na2o(£),  l673°K-^2Na(g),  l673°K 

+  o(g),  167  3°K,  Ah  =  133,380 


cal/mole 
cal/mole 
cal/mole 
cal/2  moles 
cal/mole 

cal/mole 


Ae  =  133,320  -  3  X  2  X  1673  =  123,300  cal/mole 


Similar  computations  fer  Al203  and  SiOg  give,  for 


A1203(£),  l673°K^2A/(g),  l673°K  +  30(g),  1673°K,Ah  =  696,700  cal/mole 

Ab  =  680,000  cal/mole 


65 


and;  for 


SiOg(£),  l673°K-s>Si(g),  167 3°K  +  20(g),  167 3°K,  Ah  =  412,200  cal/mole 

Ae  =  402,200  cal/mole 

For  the  dissociation  of  the  molecule,  Na20*Alg03'6Si02(£.)  as  pure 
single  component  oxides  to  gas  atoms  at  l673°K,  Ae  =  123i700  cal 
per  mole  of  gas  atoms  formed.  The  energy  of  dissociation  of  the  air 
molecules  0.8  N2(g)  and  0.2  02(g)  at  l673°K  is  79,400  cal/mole.  Thus, 
for  the  second  process, 

-  Qg  *  2-95  x  10%  x  7.94  x  10 ^  +  1.93  X  10%  X  1.24  X  105 

=  4.73  X  101;LW  cals  (185) 

The  radiant  energy  lost  for  the  2.3  MT  yield  surface  burst  is 

=  0.290  X  1.11  X  1012W  (186) 

=  3.22  X  101:LW  cals 

Eq.  186  applies  only  for  W  equal  to  2.3  MT;  for  other  yields  Eq.  178 
(divided  by  1.11  X  101%)  can  used  to  estimate  Q^. 

The  fireball  volume  at  the  second  maximum,  from  Eq.  160,  is 

V2  53  T-92  X  101:LW  cm3  (187) 

and  the  volume  at  maximum  expansion,  from  Eq.  164,  is 

Vm  «  1.64  X  1012W0’9^2  cm3  (188) 

The  change  in  volume  is 

Av  -  1.64  x  io12w  [w"0,058-o.482]  (189) 
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or;  for  the  2.3  MT  yield  (w‘°*°58  a  0.638), 

■  Vm  =  1.05  X  1012W  cm3  (19O) 


and 

Av  =  2.56  X  101;LW  cm3  (191) 


The  wurk  2.5 sociaheu.  with  "the  expuue ion  of  "the  ij.ireiJti.Lx 

gases  against  1  atmosphere  is  (l  cal  =  41.29  cm3 -atmospheres) 

%  =  po  Av  =  0.062  X  lO1^  cals  (192) 

When  the  fireball  reaches  its  maximum  volume  in  expanding 
against  the  atmosphere,  the  total  internal  pressure  in  the  fireball 
should  be  very  close  to  the  external  pressure,  which,  for  the  model 
surface  burst,  has  been  taken  to  be  1  atmosphere.  The  total  number 
of  air  molecules  (plus  a  small  amount  of  gaseous  detonation  products) 
in  the  fireball  at  l400°C,  by  use  of  the  perfect  gas  law,  is  then 

lUj,  «  1.20  X  loV^2  moles  (193) 

or,  for  2.3  MT  (W-0”^'’^  _  0.638), 

n^  «  7*66  X  10^W  moles  (194) 

With  the  exception  of  the  NagO  component,  the  partial  pressure  of  the 
"soil"  at  its  melting  point  should  be  negligible  in  comparison  to  a 
total  pressure  of  1  atmos. 

The  number  of  air  molecules  that  have  entered  the  fireball  up 
to  this  time,  neglecting  the  NagO  soil  component,  is 

n^  =  n^p  -  n£/2  (195*0 

a  1.20  X  107W  [r°-°58  -  0.123]  moles  (195b) 
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or,  for  2,3  MT  (w-0*°58  -  0.638) 


=  6.18  X  10^W  moles 


(196) 


This  number  of  moles  of  air  is  about  8l  $  sf  the  gas  molecules 
present  in  the  fireball. 

The  change  in  internal  energy  of  the  air  molecules  in  being 
heated  from  298  °K  to  1673  °K  is  10,800  cals/mole.  The  total  change 
in  internal  energy  for  this  is 

Ae5  =  1.30  X  10i;Ltf  [w-°-°58  -  O.123]  cals  (197) 


For  the  ng  moles  of  air  initially  at  298  °K  and  1  atmos,  the  energy, 
or  work,  required  for  their  expansion  to  v’m  (against  1  atmos)  is 


q5  =  3.27  X  1010W  [w-°-°58  +  0.027]  cals 


(198) 


.  Some  energy  is  also  required  for  the  mechanical  mixing  of  the 
incoming  gases  with  the  gases  already  present  in  the  fireball.  The 
postulated  assumption  that  the  fireball  volume  remains  constant  after 
tm  until  the  gases  are  near  ambient  (i.e.  500  °K  or  lower)  would  re¬ 
quire  either  that  external  gases  enter  before  tm  or  that  the  internal 
pressure  he  larger  than  1  atmos.  The  yield- independent  perfect  gas 
law  estimate  of  the  pressure  at  tp  is  3*3  atmos.  If  no  additional 
gases  enter  before  the  pressure  at  would  be  0.84  y0*10o  a^rao s , 
or  1.9  atmos  for  the  2.3  MT  yield.  The  temperature  at  this  time  is 
3.46  X  103  w0.048  °Kj  or  about  5000  °K  for  2,3  MT  so  that  recombina¬ 
tion  of  the  gas  atoms  would  be  occurring  at  this  time.  By  the  time 
the  temperature  decreased  another  1000  °K,  most  of  the  0  and  N  atoms 
would  have  recombined  (note  that  the  formation  of  nitrogen  oxides 
has  not  been  considered  in  this  treatment)  to  form  Op  and  Np  atoms; 
also,  many  of  the  metallic  atoms  would  have  formed  oxide  molecules. 
The  recombination  should  result  in  a  drop  in  the  pressure  (between 
4000  and  5000  °K)  of  about  30'  $  not  considering  the  drop  in  tempera¬ 
ture  .  At  about  2500  °K,  the  soil  constituents  will  start  condensing 
into  liquid  droplets.  This  process  should  also  tend  to  reduce  the 
pressure  further. 

Assuming  complete  recombination  when  the  fireball  has  cooled 
to  4000  °K,  the  number  of  gas  molecules  of  the  material  originally 


68 


present  at  tg  would  be  3.4l  X  l(Af  moles  so  that  the  number  of  addi¬ 
tional  air  molecules  would  be,  from  use  of  the  perfect  gas  law. 


n3  (4000°K)  -  5.00  X  lO^W  (PW-0-°58  -  0.682)  moles  (199) 

where  P  is  the  total  pressure  of  the  gases  in  atmospheres.  For  no  ' 
equal  to  zero,  the  pressure  would  be  about  1.1  atmos  for  the  2.3  MT 
yield.  If  the  additional  gas  were  as  much  as  25  %  of  the  total 
present,  the  pressure  would  be  about  1.4  atmos  (2.3  MT) . 

For  a  reliable  estimate  of  the  mixing  energy,  the  variation  of 
the  gas  concentration  of  content  in  the  fireball  (e.g.,  the  total 
pressure)  with  time  is  required.  Since  that  cannot  be  deduced  from 
the  scaling  equations,  the  end-point  composition  will  be  used  to  esti¬ 
mate  the  entropy  and  neat  of  mixing  for  a  process  in  which  the  gases 
are  at  the  same  initial  temperature  and  pressure  when  mixing  occurs. 
For  the  computation,  the  mixing  process  will  be  assumed  to  occur  at 
3300  °K  (mid-temperature  of  5000  °K  at  tm  and  1673  °K)  the  heat  of 
mixing  is  then  given  by  ' 


q*  =  3300  rvp  Asffl  (200) 

where  As^  is  the  entropy  of  mixing  at  constant  temperature  and  pres¬ 
sure.  This  computation  can  be  interpreted  either  as  a  recognition 
that  the  gases  are  nonideal  and  some  heat  of  mixing  is  involved,  or, 
that  the  gases  are  nearly  ideal  and  the  entropy  change  due  to  mixing 
at  the  selected  temperature  is  essentially  a  measure  of  the  reversible 
work  required  to  do  the  mixing.  The  overall  value  of  Asm  (constant 
T  and  P)  for  the  estimation  is  defined  by 

Asm=  -S  [(n^/n^Zn^/r^)  +  (n^n^n^/n^)]  (201) 

For  the  final  mixture,  Asm  is  -O.98  cals/mole  and 

q,L  =  3.88  X  1010W0’^2  cals  (202) 

The  total  energy  for  the  heating,  expansion,  and  mixing  of  the  air, 
is  the  sum  of  Eqs.  197,  198,  and  202  and  is  given  by 
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=  2.02  X  101;lW  [w“°*058-O.OT5] 

or,  for  2.3  MT  (w=0'°58  =  O.638) 

Q*  =  l.l4  X  10i:LW 
2 


(203) 


(204) 


The  summation  of  Eqs.  184,  185,  186,  192  and  204  gives  a  residual  of 
1.02  X  lO-1--1-  W  cals  for  heating  and  melting  soil.  If  the  temperature 
were  uniform  throughout  the  fireball  and  if  the  soil  particles  melted 
instantaneously,  only  liquid  soil  would  exist.  This  condition  un¬ 
doubtedly  does  not  occur  because  of  the  low  thermal  conductivity  of 
soils .  Many  of  the  particles  that  enter  the  fireball  at  the  time  the 
highest  temperature  for  the  existence  of  a  liquid  phase  is  reached 
(about  2500  °K)  would  be  completely  melted  (it  takes  about  33  secs 
for  the  temperature  to  drop  from  25OO  °K  to  1673  °K  for  the  2.3  MT 
yield)  while  those  that  were  entering  just  as  the  melting  temperature 
occurred  would  not  be  melted  at  all.  The  latter  would,  however,  be 
heated  to  some  degree  and  therefore  would  have  absorbed  some  of  the 
available  energy  —  as  would  the  solid  phase  in  the  interior  of  other 
particles  that  entered  at  intermediate  times  and  were  partially  melted. 
Considering  that  more  soil  would  be  vaporized  between  tg  (1-5  sec)  and 
the  time  at  which  the  temperature  has  decreased  to  2500  °K  at  t„  (27 
secs)  together  with  the  melting  of  additional  soil  from  2500  °K^to  . 
1673  °K,  the  fraction  of  available  energy  for  melting  should  lie 
between  0.25  and  0.75*  Thus,  the  final  proportionment  of  the  energy, 
like  others  previously,  is  an  assigned  amount ]  the  fraction,  0.5  of 
the  . residual  energy,  or  O.5I  X  10^  W  cals'for  the  2.3  MT  yield,  was 
taken  to  make  the  estimate  of  the  number  of  moles  of  liquid  soil  pre¬ 
sent  in  the  fireball  at  the  selected';:soii  melting  point. 

‘i 

The  change  in  heat  content  of  soil  minerals  such  as  the  ideal 
soil  defined  previously,  for  a  rise;  in  temperature  from  298  °K  to 
1673  °K,  would  be  about  200,000  cal /mole .  The  heat  of  fusion  would 
be  about  50>000  cal/mole.  The  nuniber  of  moles  of  melted  soil  is  then 


n(£)  =  2.04  X  10^W  +  0.74  X  10^W  moles 
-  2.78  X  lo^w  moles 


(205a) 

(205b) 


And,  for  the  2.3  MT  yield,  the  average  concentration  of  the  liquid 
phase  in  the  fireball  is 
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n(£)/v  =  2.65  X  IQ"?  moles/cm3  (206) 

The  value  of  (n(/.)/v)RT  (T  =  l673°K)  is  then  3*64  X  10“^  atmos¬ 
pheres;  it  is  used  in  Section  7  to  compute  the  fraction  of  each  fission 
product  element  that  has  condensed  into  the  liquid  phase  of  the  ideal¬ 
ized  carrier  material. 

Not  considering  the  amount  of  soil  vaporized  initially  (up  to 
which  by  Eq.  205  was  about  0.27  of  n(j£),  only  4.6  $  of  the  total 
energy  was  used  to  melt  soil  for  the  2.3  MT  yield  surface  burst.  For 
a  low  yield  tower  shot,  Adams0  estimated  that  about  3  of  the  total 
energy  was  used  in  heating  soil  and  tower  materials.  The  total  energy 
content  of  the  liquid  soil  at  its  melting  point  therefore  represents 
about  6.3  <f„  of  the  total.  According  to  Eq.  205 >  the  total  mass  of 
liquid  soil  is  3-4  X  10-^  gm  (about  380,000  tons);  this  mass,  from 
Eq.  88  represents  about  5*4  °jo  of  the  mass  thrown  out  of  the  crater 
by  the  2.3  MT  yield  detonation. 

No  argument  can  be  given  regarding  the  accuracy  of  the  fire¬ 
ball  scaling  functions  and  the  thermal  data  used  to  arrive  at  the 
final  value  of  n(/)/v.  Certainly  many  improvements  can  be  made  with 
more  detailed  treatment  of  the  many  parameters  when  more  of  the  needed 
data  become  unclassified.  The  treatment  of  the  data  in  this  section, 
although  preliminary  in  scope,  has  identified  some  types  of  information 
and  parameters  required,  and  has  attempted  to  place  them  in  a  concep¬ 
tual  framework  with  respect  to  their  role  in  the  fallout  formation 
process. 
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SECTION  5 


•nT^n-TA»T  nr*p‘nTT'n  -rv^Ti  nTOOTAM  r\TT»  tt235 

r  xoDluiM  xT.hitUj  run  r  iuuxuw  vr  u  ; 


U‘ 


■238 


AND  Pi 


,239 


One  major  factor  that  determines  the  amount  of.  a  fission  product 
that  condenses  during  the  first  period  of  condensation  is  the  amount 
of  that  element  present  in  each  of  the  mass  chains  when  the  carrier 
solidifies.  According  to  the  equations  given  in  Section  2,  the  frac¬ 
tional  chain  yield  of  each  element  is  required  for  estimating  the 
"R"  factor  for  the  chain  relative  to  the  number  of  fissions,  or  to 
the  yield  of  Mo99  or  some  other  reference  nuclide. 

The  independent  yields  of  all  the  radionuclides  have  been  cal¬ 
culated  by  Bolles  and  Ballou^  for  thermal  neutron  fission- of  U235 
according  to  the  theories  of  independent  fission  yields  of  Glendenin1' 
and  of  Present. The  fractional  chain  yields  from  these  calculations 
are  given  from  0  to  189  sec  in  Table  5  for  mass  89  and  in  Table  6  for 
mass . 140;  each  chain  contains  a  rare  gas  (Kr  and  Xe)  element.  The 
two  mass  chains  selected  lie  on  the  outer  edge  of  the  two  yield  peaks . 

In  each  case,  the  independent  yield  distribution  according  to  the 
Glendenin  theory  is  the  broader,  with  the  short-lived,  lowest-Z 
elements  more  heavily  weighted.  This  difference  diminishes  as  the 
mass  number  approaches  the  value  118. 

The  fractional  yields  themselves  indicate  what  might  be  expected 
,  as  to  the  fractionation  of  the  chain  members  during  condensation. 

For  example,  if  the  rare  gas  member  only  is  considered,  the  maximum 
loss  ( minimum  amount  condensed)  should  occur  for  mass  89  between  19 
and  28  secs  (Glendenin) ;  and  for  mass  lUO,  the  fraction  not  condensed 
(for  both  theories  of  yield)  should  decrease  as  the  time  period  of 
the  condensation  process  increases. 

Similar  calculations  have  not  been  made  fox-  the  independent  yields, 
for  fission  of  U235  and  Pu239  with  fission- spectrum  neutrons  nor  for 
U238  fis  sion  with  8  Mev  broad-band  spectrum  neutrons  which  would  be 
more  applicable  to  nuclear  detonations.  Moreover,  the  chain  yields 
themselves  for  Pu239  and  U238  fission  are  not  very  well  known.  Most 
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TABLE  5 

Fractional  Chain  Yields  for  Mass  89  for  U^35  Fission 
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of  the  available  data  and  estimates  of  unmeasured  yields  are  summar¬ 
ized  in  Table  7  for  the  fission  of  11^35,  1J238,  and  Pu239 .  The  bulk 
of  the  data  were  taken  from  the  summary  by  Katcoff .^-9  The  yield 
curve  for  thermal  neutron  fission  of  U235  gives  an  average  value 
of  2. 5  neutrons  per  fission.  For  fission-spectrum  neutron  fission 
of  U238,  an  average  value  of  3-1  neutrons  per  fission  is  obtained. 

These  values,  together  with  the  referenced  data  in  Table  1,  were 
used  to  obtain  a  set  of  stylized  yield  curves  for  fission-neutron 
fission  of  U235  and  Pu239  and  the  8  Mev-neutron  fission  of  U238. 

In  each  case,  the  fission  by  fission-neutrons  was  assumed  to  yield 
about  3  neutrons  per  fission  and  the  8  Mev-neutron  fission  to  yield 
4  neutrons  per  fission.  For  this  change  in  neutron  yield  per  fission 
with  neutron  energy,  l4-Mev  neutron  fission  would  give  about  4-5  neu¬ 
trons/fission.  The  increase  in  the  neutron  energy  from  thermal  to 
fission-spectrum  energies  shifted  the  right  side  of  the  heavy-element 
peak  to  heavier  masses  by  about  l/4  mass  unit  (for  the  same  yield). 

Where  there  was  no  data  to  indicate  possible  changes  in  the  "fine" 
structure  shape  of  the  yield  curve  at  the  peak,  the  general  shape 
of  the  thermal-neutron  yield  curves  was  retained  but  adjusted  in 
height  so  that  the  peak  sum  was  reasonably  near  100.0.  The  large 
discrepancy  in  the  rare  earth  yields  for  thermal  neutron  fission 
of  Pu239,  between  Katcoff 's  values  and  those  given  by  Bunney,20 
may  be  due  to  persistent  errors  in  the  counter  calibration.  The 
yield  values  given  by  Bunney  are  lower  by  a  factor  of  1-5  at  mass 
number  144  and  approach  Katcoff' s  values  as  the  mass  increases;  at 
mass  number  156,  Bunney* s  value  is  only  about  20  percent  lower.  The 
values  given  by  Katcoff  give  a  peak  sum  nearer  100  and  were  therefore 
used  in  Table  7 • 

Comparison  of  the  cumulative  chain  yields  at  the  two  peaks  (A  = 

90  to  100  and  A  =  131  to  144)  for  the  various  types  of  fission  in¬ 
dicates  that  no  very  large  differences  in  the  decay  rates  should  be 
expected.  Of  the  three  mentioned  yield  curves,  the  yields  for 
fission-neutron  fission  of  Pu239  appear  to  differ  most  from  those  of  ther 
mal  neutron  fission  of  U235.  por  mass  numbers  such  as  140  and  95  whose 
radioisotopes  may  contribute  more  than  80  percent  of  the  total  gamma 
radiation  at  specific  times  after  fission,  the  largest  percentage 
difference  is  5*0/6.44,  or  28  percent  (lower)  for  mass  number  iMj 
and  5*6/6.27,  or  11  percent  (lover)  for  mass  number  95*  The  yield 
of  mass  number  90,  however,  is  significantly  lower  for  U238  (8  Mev 
neutrons)  and  Pu239  (fission  neutrons);  the  percentage  differences 
are  36  percent  (lower)  and  48  percent  (lower),  respectively .  The 
yields  of  mass  number  137  are  all  more  nearly  alike.  The  largest 
differences,  of  course,  occur  for  the  mass  numbers  in  the  valley 
between  the  peaks  and  those  at  the  highest  and  lowest  mass  numbers . 

But  even  for  U238  fission  (8  Mev  neutrons),  v/here  the  yields  of  the 
mass  numbers  near  118  are  more  than  a  factor  of  ten  larger  than  for 
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U^35  fission  (thermal  neutrons) ,  the  contribution  of  these  elements 
will  be  a  small  fraction  of  the  total  activity. 

The  ,!R,!  factors  due  to  the  different  yields  will  be  more  sen¬ 
sitive  to  small  changes  in  abundances  than  the  gross  decay  curve. 
Using  the  U^35  fission  with  thermal  neutrons  as  a  standard,  the 
yield  "R"  factors  for  mass  l4h,  for  example,  are  O.89  (U?35,  fis¬ 
sion  neutrons),  0.7^  (u238,  Q  Mev  neutrons),  and  O.87  (Pu239, 
fission  neutrons);  for  mass  number  115  the  respective  yield  "R" 
factors  are  1.6,  14,  and  6.8. 


For  the  computations  presented  in  the  next  two  sections,  the 
independent  yields  as  calculated  by  Bolles  and  Ballou! 3  for  thermal 
neutron  fission  of  U^35  were  converted  to  fractional  chain  yields. 
The  same  fractional  chain  yields  were  then  used  for  all  fission 
yield  curves . 


Although  it  appears  that  Glen&enin's  symmetrical  charge -distri¬ 
bution  curve  is  generally  applicable  for  all  fissile  nuclides  in 
low-energy  neutron  fission,  it  has  been  shown  that  the  most  probable 
charge,  Zp,  for  a  given  mass  split  shifts  toward  stability  with  in¬ 
creasing  neutron  energy;  i.e.  the  higher  fractional  yields  would 
appear  farther  to  the  right  (higher  z)  in  any  decay  chain. 21  Pappas22 
used  a  discontinuous  function  for  Zp  and  considered  the  primary  frag¬ 
ments  before  neutron  boil-off';  Wahl  2 3  has  shown  empirically  however 
that  the  Zp  function  in  thermal  fission  of  U235  is  continuous,  as 
originally  postulated  by  Glendenln  et  al . 24  Herrington25  proposes 


two  charge-distribution  curves,  one  for  even-neutron  nuclides  and 
another  showing  lower  yield  for  odd-neutron  nuclides. 


It  is  clear  that  there  is  no  unequivocal  choice  in  methods  for 
estimating  the  independent  yields  of  the  chain  members,  even  for 
thermal  neutron  fission  of  lj235.  For  higher-energy  fission  the 
experimental  data  on  independent  yield  are  rarer  yet.  It  was  there¬ 
fore  assumed  that  the  fractional  independent  yields  calculated  by 
Bolles  and  Ballou  on  Glendenin’s  postulate  for  thermal  fission  of 
tj23  5  would  not  be  too  inappropriate  for  any  kind  of  low-energy 
fission.  When  more  data  on  the  total  chain  yields  and  the  inde¬ 
pendent  yields  are  available  for  each  fissile  nuclide,  the  indicated 
corrections  can  he  applied  to  the  computations. 
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TABLE  7 


Cumulative  Mass-Chain  Yields  of  Fission  Products 
(Values  are  in  percent  of  fissions) 


Mass 

u235 

U238 

Pu239 

Number 

Thermal 

Neutrons* 

Fission 

Neutrons 

Fission  8-Mev 

Neutrons  Neutrons 

Thermal 

Neutrons 

-  Fission 
Neutrons 

72 

1.6x10" 5 
1.1x10  , 

/  ^  ~  ~  _iL\  PL 

4. 6x10 

5. 0x10" 6 

1.2x10-4* 

73 

0.0012 

3-7x10-5 

- 

2.2x10-4 

- 

74 

0 . 0034 

1 .  Ixl0":!' 

0.001 

i.  -i  *.  x>_Ll 

/  • 

A  AAn  1 

U  .  U'Jll 

75 

(8.8x10-4) 

0.0062 

8.3xlO"4 

o.oo4o 

7 • 5xl0"4 

0.0023 

76 

(0.0029) 

0.012 

0.0012 

0.0078 

o.ooi4 

0.0051 

77 

O.OO83 

0.023 

O.OO38* 

o.oi4 

0.0026 

0.011 

78 

0.021 

0.048 

O.OO95 

0.026 

0.0049 

0.025 

79 

(0.041) 

0.096 

0.019 

0.053 

0.0090 

0.043 

80 

(0.077) 

0.19 

0.045 

0.096 

0.016 

0.075 

81 

o.i4 

0.21 

0.088 

0.18 

0.030 

0.l4 

82 

(0.29)  ■ 

0.50 

0-.20 

0.35 

0.056 

0.23 

83 

0.544 

0.80 

0.4o* 

0.66 

0.10 

0.37 

84 

1.00 

1-3 

O.85* 

1.02 

0.17 

0.60 

85 

1.30 

I.85 

0.80 

1.45 

0.28 

0.92 

86 

2.02 

2.5 

1.38* 

1-9 

0.45 

1.15 

87 

(2-94) 

3-3 

1.90  . 

2.25 

0.73 

1*5 

88 

(3*92) 

4.2 

2.45 

2.7 

1.2 

1-9 

89  4-.  79 

Continued 

5-1 

2.9* 

3-17 

1.9* 

2.4 

*Seymour  Kstcoff,  Fission-Product  Yields  From  U,  Th  and  Pu, 
Nucleonics,  Vol.  16,  No.  4,  p.  78-85  ( 1958)4 


**L.R.  Bunney,  E .M.  Scadden,  J.O.  Abriara  and  N-E.  Ballou,  Radio- 
chemical  Studies  of  the  Fast  Neutron  Fission  of  U^35  and  u2j8^ 
Second  UN  International  Conference  on  the  Peaceful  Uses  of 
Atomic  Energy,  A/Conf.  15/P/643,  USA,  June  1958. 


**■*0 .p .  Ford,  J.S.  Gilmore,  et  al,  Fission  Yields,  LADC-3083* 
****L.R.  Bunney,  E.M-  Scadden,  J.O*  Abriam,’  N-E-  Ballou.  Fission 
Yields  in  Neutron  Fission  of  Pu239,  USNRDL-TR-268,  19 58,  Unci. 


a.  Parentheses  indicate  estimated  values  or  where  Katcoff ’ s  value 
was  altered  in  order  to  adjust  the  yields  to  a  gross  sum  of  100 
in  each  peak. 

b.  Line  indicates  division  of  two  peaks  that  was  used  for  individual 
peak  sums. 
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TAULti  7  (Cont'd) 


Cumulative  Mass-Chain  Yields  of  Fission  Products 
(Values  are  in  percent  of  fissions) 


Mass 

Number 


U235 

Thermal  Fission 
Neutrons*  Neutrons 


Neutrons  Neutrons 


Pug39 _ 

Thermal  Fission 
Neutrons  Neutrons 


90 

5-77 

5.8 

3-2* 

3-7 

2.4 

3-0 

91 

5.84 

5.85 

3.6 

4.3 

3.0 

3.7 

92 

6.03 

6.0 

4.1 

4.8 

3-7 

4.4 

93 

6.45 

6.4 

4.85 

5 ‘2 

4.6 

5-0 

94 

6.40 

6.4 

5-3 

5.45 

5-5 

5.4 

95 

6.27 

6-3 

5.7* 

5.6 

5.9* 

5.6 

96 

6-33 

6-3 

5-8 

5.7 

5-7 

5.3 

97 

6.09 

6.1 

5-7 

5-64 

5-6* 

5-2* 

98 

5-78 

5.8 

5-7 

5.6 

5.4 

5.4 

99 

6.06 

6.1** 

6.3* 

6.2** 

5.9* 

5.9* 

100 

6.30 

6.7 

6.1 

6.4 

6.0 

6.4 

101 

'5.0 

5-3 

5-5 

C.  c 

6.0 

5-9 

102 

4.1 

2-9 

5.6 

5-9 

5-9 

5-3 

103 

3-0  . 

1-7 

6.6 

5.0 

5.8* 

4.6 

104 

1.8 

0.95 

5-4 

3-2 

5.0 

3-5 

105 

0.90 

0.54 

3-9 

2.2 

3-9* 

3-2 

106 

0.38 

0.30 

2.7* 

1.5 

5.0* 

3.6 

107 

0.19 

0.17 

1-35 

1.0 

4.0 

3-1 

108  1 

(0.085) 

0.095 

0.67 

0.70 

3-0 

2.6 

109  1 

(0.039) 

0.053*** 

0.32* 

0.48 

1.5* 

1.9*- 

110  1 

(0.020) 

0.030 

0.15 

0-33 

0.65 

0.81 

ill  1 

(0.015) 

0.022*** 

0.073* 

0.23*** 

0.27* 

0.34 

112  1 

(0.013) 

0.020*** 

0.046* 

0.19 

0.10* 

o.i4* 

113  1 

(0.012) 

0.018 

0.043 

0.17 

0.055 

0.090 

114  1 

(0.011) 

0.017 

0.04l 

0.16 

0.046 

0.075 

115 

0.0104 

0.017*** 

0.040* 

0.15*** 

0.041* 

0.069* 

116  1 

f 0.010)6 

0.017b 

0.039 

0.14 

0.039 

0.065 

117  1 

,0.010) 

0-017 

0.039^ 

0.146 

0.038 

o.o64 

118  I 

'0.010') 

0.017 

o.o4o° 

0Tl4 

O.O386 

0.0646 

119  < 

(o.oii) 

0.017 

o.o4i 

0.14 

0.039 

0.064 

Continued 
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TABLE  7  (Cont'd) 


Cumulative  Mass-Chain  Yields  of  Fission  Products 
(Values  are  in  percent  of  fissions) 


U235 _  y233  _ Pu239 


Mass 

Thermal 

Fission 

Fission 

8  Mev 

Thermal 

Fission 

Number 

Neutrons* 

Neutrons 

Neutrons 

Neutrons 

Neutrons 

Neutrons 

120 

(0.011) 

0.018 

0 .042 

0.15 

0.041 

0.065 

121 

(0.012) 

0.020 

0.044 

O.lb 

o.o44* 

0.066 

122 

(0.013) 

0.022 

0.046 

0.17 

0.047 

O.O69 

123 

(0.015) 

0.030 

0.050 

0.19 

0.052 

0.076 

124 

(0.017) 

0.053 

0.055 

0.23 

0.058 

0.082 

125- 

0.021 

0.095 

0.072 

0.33 

0.072* 

0.l4 

126 

(0.058) 

0.17 

0.175 

0.48 

0.175 

0-35 

127 

(0.145) 

0.30 

0-39 

0.70 

0-39* 

0.80 

128 

0-37 

0.54 

0.77 

1.0 

0.77 

1*9 

129 

0.90 

0.95 

1-45 

1*5 

1.45 

2.5 

130 

2.0 

1=7 

2*5 

2.2 

2-5 

3*2 

131 

(2.88) 

2-9 

3.2* 

3*2 

3*8* 

3*8 

132 

(4-31) 

4-3 

4.7* 

4.4 

5.0 

4.6 

133 

(6.48) 

6.1 

5.5* 

5.4 

5.27* 

4.9 

134 

(7-80) 

7-3 

6.6* 

6.5 

5.69* 

5*2 

135 

(6.40) 

6-3 

6.0* 

5*9 

5*53* 

5*1 

136 

(6.36) 

6.4 

5.9* 

5.8 

5.06* 

5*3 

137 

(6.05) 

6.0 

6.2 

5.85 

5-24* 

6.4* 

138 

5.74 

5=7 

6.4 

5*9 

5*5 

5.4 

139 

(6.34) 

6.4 

6.5 

6.0 

5*7* 

5*2 

140 

6.44 

6.4 

5.7* 

5.6 

5.68* 

5.0* 

l4l 

(6.30) 

6.3 

5*7 

5*5 

5.2* 

4.7 

142 

(5.85) 

5-9 

5*7 

5*4 

6.69* 

4.9 

143 

(5-87) 

c*  Q 

J  .KJ 

c  e 

S  '  s 

4.97 

5.4* 

.n 
✓  ' 

144 

5-67 

5.1** 

4.9* 

4.3** 

5.29* 

4.8 

145 

3-95 

4.2 

3*7 

3*7 

4.24* 

4.4 

146 

3-07 

3.3 

3*1 

3*17 

3*53* 

3*7 

147 

2.18 

2.5** 

2.6** 

2.7** 

2.92* 

3-0 

148 

1.70 

1.85 

2.0 

2.27 

2.28* 

2.36 

149 

1.13 

1.3** 

1.45 

1-9** 

1*75 

1.86 

Continued 


TABLE  7  (Coni' d) 


Cumulative  Mass-Chain  Yields  of  Fission  Products 
(Values  are  in  percent  of  fissions) 


Mass 

u235 

U238 

Pu239 

Number 

Thermal 

Fission 

Fission 

8  Mev 

Thermal 

Fission 

Neutrons* 

Neutrons 

Neutrons 

Neutrons 

Neutrons 

Neutrons 

150 

0.67 

0.80 

1.05 

1.45 

1.38* 

1.48 

151 

0.45 

0.50 

0.74 

1.02 

1.08 

1.16 

152 

O.285 

0.31 

0.50 

0.66 

O.83* 

0.92 

153 

0.15 

0.19** 

0.32 

0.4i** 

0-52 

0.60 

154 

0.077 

O.O96 

0.19 

0.25 

0.32* 

0.37 

155 

O.O33 

0.048 

0.11 

0.15 

0.20 

0.23 

156 

0.014 

0.023** 

0.066* 

0.092** 

0.12* 

o.i4 

157 

0.0078 

0.012 

0.034 

0.057 

0.064 

0.075 

158 

0.002 

0.0062 

0.016 

0.032 

0.034 

0.043 

159 

0.00107 

0.0034** 

0.0090** 

0.017** 

0.020**** 

0.025 

i  /iV\ 

JLU^J 

3-5x10-4 

0.0012 

O.OO36 

0.0085 

0.0092 

0.011 

l6l 

7-6x10-5 

4 . 6x10 -***9 • kxlO-4 

o.oo44** 

O.OO38****  0.0051 

*Seymour  Katcoff,  Fission-Product  Yields  From  u,  Th  and  Pu, 
Nucleonics ,  Vol.  l6,  No.  V,  p  .  78-85  ’(1958)  *  °~ 

**L-R-  Bunney,  E«M  Scadden,  J.O.  Abrieun  and  N.E.  Ballou,  Radio¬ 
chemical  Studies  of  the  Fast  Neutron  Fission  of  U235  and  U238, 
Second  UN  International  Conference  on  the  Peaceful  Uses  of 
Atomic  Energy,  A/Conf.  15/p/ 643,  USA,  June  1958. 

***G.P.  Ford,  J.S.  Gilmore,  et  al,  Fission  Yields,  LADC-3083- 
****L-R-  Bunney,  E-M.  Scadden,  J.O.  Abriam,  N.E.  Ballou.  Fission 
Yields  in  Neutron  Fission  of  Pu239,  USNRDL-TR-268,  1958,  Unci. 


a.  Parentheses  indicate  estimated  values  or  where  Katcoff ?s  value 
was  altered  in  order  to  adjust  the  yields  to  a  gross  sum  of  100 
in  each  peak. 

b.  Line  indicates  division  of  two  peaks  that  was  used  for  individual 
peak  sums . 
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SECTION  6 


CALCULATION  OF  DISINTEGRATION  RATES,  PHOTON  EMISSION  RATES, 
PHOTON -ENERGY  EMISSION  RATES,  AND  AIR -IONIZATION  RATES 
FROM  NORMAL  FISSION  PRODUCTS  OF  U235,  U23°,  AND  Pu239 


6.1  DESCRIPTION  OF  COMPUTATIONS 

In  making  the  computations ,  the  fractional  chain  yields  were 
taken  to  be  the  same  for  all  fissile  nuclides  as  were  discussed  in 
Section  5 •  For  this  stipulation,  the  disintegration  rate  per  fis¬ 
sion  of  a  given  fission  product  nuclide  from  each  fissile  nuclide 
differs  only  by  a  constant.  It  was  therefore  convenient  to  use 
the  disintegration  rates  calculated  by  Bolles  and  Bailout  for 
thermal  fission  of  U235  and  multiply  them  by  the  ratio  of  the 
yields  from  the  fissile  nuclide  of  interest  to  that  from  the  ther¬ 
mal  fission  of  U235j  in  order  to  obtain  the  disintegration  rates 
desired.  These  ratioB  are  the  same  as  the  yield  "R”  values  given 
in  Section  3*  The  summation  over  all  masses  time-wise  gives  the 
total  fission  product  disintegration  rate. 

Where  decay  schemes  of  individual  fission  products  are  known, 
gamma  ray  characteristics  may  also  be  computed.  If  the  particular 
characteristics  are  evaluated  per  disintegration  of  each  radionuclide, 
each  disintegration  rate  value  can  be  multiplied  by  the  appropriate 
factor  time -wise  and  the  products  for  all  nuclides  summed  to  give 
ule  bOuSU.  euax-ue  wi  is  ujuc  vx  me  mixture - 

The  disintegration  rates,  photon  emission  rates,  photon  energy 
emission  rates,  and  air  ionization  rates  were  computed  for  times 
extending  from  45-8  minutes  to  25-7  years  after  fission  for  the  3 
fissile  nuclides.  The  first  three  rate  curves  were  calculated  for 
ICh  fissions;  the  last  was  for  the  air  ionization  3  ft  above  a  smooth 
infinite  plane  uniformly  contaminated  with  10^  fissions/sq.  ft.  The 
required  basic  decay  scheme  information  was  taken  from  a  previous 
report26  and  a  current  revision^  including  nuclear  data  through 
July  1959 •  Details  of  the  method  of  calculation  of  the  air 
ionization  rates  have  been  published.  28 
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6.2  RESULTS  OF  COMPUTATIONS 


The  results  of  the  computations  are  summarized  in  Table  8  as 
decay  data  for  thermal-neutron  fission  of  U235,  fission -neutron 
fission  of  U^35  8  Mev-neutron  fission  of  U238,  thermal -neutron 

fission  of  Pu239,  and  fission-neutron  fission  of  Pu239.  The  use 
of  Katcoff 's  yields!9  (adjusted)  for  thermal  fission  give  disin¬ 
tegration  rates  that  are  almost  identical  with  those  of  Bolles 
and  Ballou.  The  air  ionizations  (u^35,  thermal)  are  also  very 
close  to  those  of  Reference  19)  at  2.4  hr,  they  are  about  8  fj  higher 
and  at  2.6  y  they  are  about  19  $  lower  (maximum  fluctuation). 

Dolan^9>30  has  calculated  the  disintegration  rates  and  photon 
emission  rates  for  lit-  Mev  neutron  fission  of  U238;  ratios  of  his 
values  to  those  in  Table  8  for  the  8  Mev  neutron  fission  of  U238 
are  given  in  Fig.  3*  It  may  be  noted  that  the  disintegration 
rates  are  within  5  $  of  each  other  from  1  to  about  350  hours;  the 
agreement  in  the  photon  emission  rates  is  not  quite  as  good;,  with 
Dolan's  values  being  more  than  10  °jo  lower  after  40  hrs.  The  maxi¬ 
mum  spread  is  +5  $  (75  h)  to  -12  $  (2500  h)  for  the  d/s  computations 
and  +3  io  (7-5  h)  to  -17  $  (1200  h)  for  the  photons/sec  computations. 
A  few  more  photons  were  counted  in  the  method  by  which  the  data  in 
Table  8  were  obtained  than  by  the  method  used  by  Dolan  (chiefly 
in  the  energy  range  0  to  20  Kev) . 


The  air  ionization-rate  curves  from  each  type  of  fission  are 
of  chief  interest;  these  are  compared  in  Figs.  4  and  5  in  terms 
of  an  air  ionization  "RM  factor.  The  factor ,  rfp,  is  the  ratio 
of  the  air  ionization -rate  from  one  type  of  fission  to  that  from 
thermal  neutron  fission  of  U235.  The  fluctuation  in  the  curves  of 
Figs.  4  and  5  reflect  the  relative  prominence  of  the  important 
gamma,  emitters  in  each  mixture .  The  deviation  in  .  from  the 
value  1  is  a  measure  of  the  difference  in  the  ionization  rate  from 
that  of  the  U^35  thermal  fission  reference  curve.  The  order  in  the 
rfp  deviations,  from  least  to  most,  is:  (l)  fission  neutron  fission 
of  U235.  (2)  8  Mev  neutron  fission  of  u23°,  (3)  fission  neutron 
fission  of  Pu239,  and  (4)  thermal  neutron  fission  of  Pu239.  The 
maximum  relative  deviation  for  the  first  three  (combined)  is  from 
-l6  j>  (2.5  h)  to  +5-5  70  (llO  h)  between  1  and  7000  hours  after 
fission.  However,  between  2  and  3  years  after  fission  the  U23.8 
(8  Mev)  r^  value  is  almost  1.6  and  the  Pu239  (fission)  r-fp  value 
is  almost  2.2,  reflecting  the  higher  yields  for  the  rare  earth 
elements  (heavy  mass  peak) . 


The  two  main  factors  that  determine  the  gross  decay  of  the 
normal  product  mixtures  (besides  the  half-lives  and  the  individual 
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TABLE  8 


Decay  of  Normal  Fission  Products  From  U235,  U238  an^  pu239 
1.  dis/sec  for  10^  fission  (Glendenin) 


Age 

U235 

U238 
(8  Mev) 

Pu239 

Years 

Days 

Hours 

Thermal 

Fission 

Thermal 

Fission 

O.763 

1.618 

I.615 

1.598 

1.558 

1.536 

i.12 

1.072 

1.076 

1.043 

1.008 

1.001 

1.64 

0.6860 

O.6933 

O.6578 

0.6253 

O.6269 

2.40 

O.4361 

6.4454 

0.4152 

O.3867 

0.3936 

3-52 

0.2818 

0.2908 

0.2691 

0.2473 

0.2568 

5.16 

0-1847 

0.1916 

0.1780 

0.1640 

0.1734 

7.56 

0.1228 

0.1271 

0.1195 

0.1117 

0.1186 

ll.l 

(1)  8171 

(1)  8393 

(1)  8000 

(l)  7661 

(1)  8036 

16.2 

(1)  5280 

(1)  5373 

(1)  5201 

(l)  5107 

(1)  5245 

23.8 

(1)  3311 

(1)  3341 

(1)  3291 

(l)  3304 

(1)  3318 

1.45 

34.8 

(1)  2037 

(1)  2042 

(1)  2049 

(1)  2092 

(1)  2062 

2.13 

51-1 

(1)  1223 

(l)  1219 

(1)  1250 

(1)  1297 

(1)  1260 

3.12 

74.9 

(2)  7417 

(2)  7381 

(2)  7698 

(2)  8062 

(2)  7751 

4.57 

(2)  4787 

(2)  4778 

(2)  4978 

(2)  5205 

(2}  4982 

6-70 

(2)  3239 

(2)  3236 

(2)  3323 

(2)  3448 

(2)  3290 

9.82 

(2)  2227 

(2)  2222 

(2)  2247 

(2)  2318 

(2)  2198 

14.4 

(2)  1523 

(2)  1502 

(2)  1501 

(2)  1548 

(2)  1453 

21.1 

(2)  1025 

(2)  1001 

(2)  1002 

(2)  1036 

(3)  9635 

30.9 

(2)  6823 

(3)  6603 

(3)  6662 

(3)  6910 

(3)  6385 

45-3 

(3)  4456 

(3)  4277 

(3)  4365 

(3)  4544 

(3)  4191  , 

66.4 

(3)  2872 

(3)  2743 

(3)  2814 

(3)  2961 

(3)  2733 

97-3 

(3)  1838 

(3)  1764 

(3)  1790 

(3)  1925 

(3)  1780 

143 

(3)  1117 

(3)  1088 

(3)  1082 

(3)  1170 

(3)  1080 

208 

(4)  6162 

(4)  6078 

(4)  5938 

(4)  7133 

(4)  6507 

301 

(4)  3192 

(4)  3174 

(4)  3140 

(4)  4175 

(4)  3722 

1.2 

438 

(4)  1676 

(4)  1666 

(4)  1710 

(4)  2506 

(4)  2185 

1.78 

650 

(5)  9854 

(5)  9766 

(4)  1025 

(4)  1538 

(4)  1341 

2.60 

(5)  6010 

(5)  5978 

(5)  6236 

(5)  8855 

(5)  7926 

3.80 

(5)  3715 

(5)  3732 

(5)  3740 

(5  4617 

(5)  4406 

5.58 

(5)  2479 

(5)  2536 

(5)  2359 

(5)  2332 

(5)  2479 

8.18 

(5)  1915 

(5)  1987 

(5)  1731 

(5)  1449 

(5)  1683 

12.0 

(5)  1509 

(5)  1589 

(5)  1323 

(5)  1036 

(5)  1253 

17.6 

(5)  1189 

(5)  1261 

(5)  1027 

(6)  7922 

(6)  9719 

25-7 

(6)  9079 

(6)  9624 

(6)  7814 

(6)  6061 

(6)  7465 

G  ontinued 
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TABLE  8  (Cont’d) 

Decay  of  Normal  Fission  Products  From  U235;  and  Pu239 

2.  betas/sec  for  1q4  fissions  (Glendenin) 


Age 

U235 

U238 
(8  Mev) 

Pu'239 

Years 

Days 

Hours 

Thermal 

Fission 

Thermal 

Fission 

0.763 

1.544 

1-543 

1-527 

1.482 

1.466 

1.12 

1.009 

1. 

,015 

0.9816 

0.9432 

0.9398 

1.64 

O.6358 

0.6444 

0.6076 

0 . 5724 

0.5767 

2-40 

0.3983 

0.4081 

0.3763 

0.3457 

0.3542 

3-52 

0.2547 

0.2637 

0.2402 

0.2170 

0.2272 

5.16 

0.1655 

0.1722 

0.1571 

0.1420 

0.1515 

7.56 

0.1088 

0.1130 

0.1042 

(D 

9556 

0.1025 

ll.l 

(1)  7159 

(1) 

7377 

(1)  6932 

(D 

6518 

(1)  6892 

16.2 

(1)  4581 

(1) 

4680 

(1)  4496 

(D 

4346 

(1)  4484 

23-8 

(l)  2866 

CD 

2905 

(1)  2860 

CD 

2840 

(1)  2853 

1.45 

34.8 

(1)  1781 

(1) 

1697 

(l)  180 4 

CD 

1831 

(l)  1802 

2-13 

51-1 

(1)  1093 

CD 

1095 

(l)  1125 

(D 

1165 

(1)  1130 

3-12 

74.9 

(2)  6837 

(2) 

7117 

(2)  6847 

(2) 

7450 

(2)  7166 

4-57 

(2)  4532 

(2) 

4558 

(2)  4690 

(2) 

4911 

(2)  4709 

6-70 

(2)  3102 

(2) 

3130 

2  3155 

(2) 

3268 

(2)  3128 

9.82 

(2)  2143 

(2) 

2158 

(2)  2121 

(2) 

2180 

(2)  2076 

14.4 

(2)  1455 

(2) 

l46o 

(2)  l4o4 

(2) 

l44o 

(2)  1360 

21.1 

(3)  9707 

(3) 

9701 

(3)  9216 

(3) 

9454 

(3)  9856 

30.9 

(3)  6373 

(3) 

6350 

(3)  5999 

(3) 

6158 

(3)  5740 

45-3 

(3)  4108 

(3) 

4083 

(3)  3852 

(3) 

3963 

(3)  3694 

66.4 

(3)  2632 

(3) 

2608 

(3)  2459 

(3) 

2559 

(3)  2390 

97-3 

(3)  1697 

(3) 

1684 

(3)  1581 

(3) 

1690 

(3)  1579 

143 

(3)  1025 

(3) 

1025 

(4)  9610 

(3) 

1075 

(4)  9900 

208 

(4)  5919 

(4) 

5920 

(4)  5587 

(4) 

6756 

(4)  6167 

301 

(4)  3113 

(4) 

3110 

(4)  3033 

(4) 

4070 

(4)  3616 

1.2 

438 

(4)  1639 

(4) 

1629 

( 4)  1 666 

(4) 

2469 

(4)  2140 

1.78  650 

(5)  9542 

(5) 

9455 

(5)  9922 

(4) 

1511 

(4)  1307 

2.60 

(5)  5708 

(5) 

5680 

(5)  5921 

(5) 

8595 

(5)  7605 

3-80 

(^)  ^Atq 

(5) 

344? 

(5)  3444 

(5) 

4364 

(5)  7605 

5  •  58 

(5)  2194 

(5) 

2257 

(5)  2080 

(5) 

2089 

(5)  2181 

8.18 

(5)  1843 

(5) 

1723 

(5)  1469 

(5) 

1218 

(5)  1401 

12.0 

(5)  1255 

(5) 

1343 

(5)  1082 

(6) 

8220 

(6)  9919 

17.6 

(6)  9621 

(5) 

I04l 

(6)  8124 

(6) 

6002 

(6)  7378 

25.7 

(6)  7129 

(6) 

7737 

(6)  5974 

(6) 

4412 

(6)  5456 

Continued 
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TABLE  8  (Cont'd) 

Decay  of  Normal  Fission  Products  from  U235,  u238,  and  Pu239 
3.  photons/sec  for  lCr  fissions  (Glendenin) 


Age 

u235 

U238 
(8  Mev) 

Pu 

239 

Years 

Days 

Hours 

Thermal 

Fission 

Thermal 

Fission 

0.763 

1.937 

1.934 

1.928 

1.965 

1.921 

1  IO 

1.280 

1  070 

“i  otrr» 

*1  0  Cr\ 

i  otrr* 

1.64 

0.8049 

0.8068 

0.7821 

0.7688 

0.7658 

2.40 

0.4908 

0.4977 

0.4751 

0.4495 

0.4581 

3.52 

0.2971 

0.3060 

0.2904 

0.2654 

0.2798 

.5.16 

0.1815 

0.1893 

0.1810 

0.1628 

0.1776 

7.56 

0.1149 

0.1204 

0.1165 

0.1052 

0.1166 

11.1 

(1)  7584 

(1)  7901 

(1)  7741 

(1)  7156 

(1)  7808 

16.2 

(1)  4989 

(1)  5146 

(1)  5090 

(1)  4866 

(1)  5133 

23.8 

(1)  3239 

(l)  3310 

(1)  3307 

(1)  3268 

(1)  3326 

1.45 

34.8 

(1)  2106 

(l)  2140 

(1)  2151 

(1)  2180 

(1)  2157 

2.13 

51.1 

(1)  1355 

(1)  1372 

(1)  1355 

(l)  1444 

(1)  1398 

3-12 

74.9 

(2)  8799 

(2)  8907 

(2)  9118 

(2)  9597 

(2)  9145 

4-57 

(2)  5886 

(2)  5968 

(2)  6115 

(2)  6461 

(2)  6097 

6.70 

(2)  3927 

(2)  3974 

(2)  4050 

(2)  4267 

(2)  3994 

9.82 

(2)  2582 

(2)  2590 

(2)  2626 

(2)  2753 

. (2)  2551 

l4.4 

(2)  1644 

(2)  1629 

(2)  1648 

(2)  1720 

(2)  1577 

21.1 

(2)  1020 

(3)  9949 

(2)  1019 

(2)  1064 

(3)  9640 

30.9 

(3)  6166 

(3)  5895 

(3)  6212 

(3)  6526 

(3)  5849 

45.3 

(3)  3573 

(3)  3334 

(3)  3665 

(3)  3895 

(3)  3467 

66.4 

(3)  2014 

(3)  2100 

(3)  1832 

(3)  2268 

(3)  2015 

97-3 

(3)  1165 

(3)  1053 

(3)  1213 

(3)  1330 

(3)  1188 

143 

(4)  6400 

(4)  5830 

(4)  6540 

(4)  7310 

(4)  6540 

208 

(4)  3342 

(4)  3179 

(4)  3335 

(4)  38IO 

(4)  3460 

301  1 

(4)  1421 

(4)  1393 

(4)  1431 

'(4)  1728 

(4)  1568 

1.2 

438 

( 5)  5158 

f 5)  5176 

(5)  5631 

(5)  7512 

(5)  6714 

1.78  650 

(5)  2117 

(5)  2158 

(5)  2647 

(5)  3752 

(5)  3334 

2.60 

(5)  1109 

(5)  1145 

(5)  1534 

(5)  2058 

(5)  1884 

0  0s\ 
J  *w 

(6)  6058 

(6)  6348 

(6)  9011 

(5)  1040 

( 5)*  1024 

5.58 

(6)  3762 

(6)  3948 

(6)  5505 

(6)  5079 

(6)  5644 

8.18 

(6)  3014 

(6)  3080 

(6)  3893 

(6)  3199 

(6)  3861 

12.0 

(6)  2709 

(6)  2690 

(6)  3088 

(6)  2583 

(6)  3153 

17.6 

(6)  2413 

(6)  2357 

(6)  2559 

(6)  2238 

(6)  2709 

27.7 

(6)  2074 

(6)  2016 

(6)  2141 

(6)  1919 

(6)  2310 

Continued 
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TABU!  8  (Cont'd) 

Decay  of  Normal  Fission  Products  From  U235,  \jQ3&,  and  Pu239 
4.  photon-Mev/sec  for  K>4  fissions 


Age  U235  U238  Pu£39 


Years 

Days 

Hours 

Thermal 

Fission 

(8  Mev) 

Thermal 

Fission 

0.763 

1 

.856 

1 

•853 

1 

.720 

1 

•633 

1 

605 

1.12 

1 

.243 

1 

.240 

1 

.147 

1 

.084 

1 

065 

1.64 

0 

•7757 

0 

.7770 

0 

.712? 

0 

.6656 

0 

6550 

2.40 

0 

.4550 

0 

.4614 

0 

.4170 

0 

.3813 

0 

3795 

3  >52 

0 

.2584 

0 

.2664 

0 

•2379 

0 

.2118 

0 

2170 

5.16 

0 

.1463 

0 

.1529 

0 

.1368 

0 

.1190 

0 

1271 

7.56 

(1) 

8549 

(1) 

8983 

(1) 

8172 

(1) 

7044 

(1) 

7786 

11.1 

(1) 

5212 

(1) 

5454 

(1) 

5096 

(1) 

4463 

(1) 

4940 

16.2 

(1) 

3121 

(1) 

3239 

(1) 

3U3 

(1) 

2811 

(1) 

3053 

23.8 

(1) 

1817 

(1) 

1868 

(1) 

1844 

(1) 

1731 

(1) 

1826 

1.45 

34.8 

(1) 

1068 

(1) 

1090 

(1) 

1093 

(1) 

1063 

(1) 

1086 

2.13 

51-1 

(2) 

6328 

(2) 

6431 

(2) 

6521 

(2) 

6593 

(2) 

6495 

3.I2 

74.9 

(2) 

3897 

(2) 

3974 

(2) 

4064 

(2) 

(2) 

4234 

(2) 

4o44 

4.57 

(2) 

2570 

(2) 

2643 

(2) 

2695 

2845 

(2) 

2669 

6.70 

(2) 

1751 

(2) 

1810 

(2) 

1817 

(2) 

1918 

(2) 

1778 

9.82 

(2) 

1199 

(2) 

1236 

(2) 

1216 

(2) 

1275 

(2) 

1170 

14.4 

(3) 

8006 

(3) 

8199 

(3) 

7926 

(3) 

8247 

(3) 

7484 

21.1 

(3) 

5202 

(3) 

5273 

(3) 

5068 

(3) 

5253 

(3) 

4726 

30.9 

(3) 

3226 

(3) 

3229 

(3) 

3125 

(3) 

3252 

(3) 

2911 

45-3 

(3) 

1862 

(3) 

1831 

(3) 

1809 

(3) 

1906 

(3) 

1706 

66.4 

(3) 

1038 

(4) 

9996 

(3) 

1011 

(3; 

1086 

(4) 

9784 

97-3 

(4) 

6179 

(4) 

5891 

(4) 

5976 

(4) 

6537 

(4) 

5946 

143 

(4) 

3540 

(4)' 

3410 

(4) 

3390 

(4) 

3790 

(4) 

3500 

208 

(4) 

1995 

(4)' 

1945 

(4) 

1894 

(4) 

2150 

(4) 

1972 

301 

(5) 

8161 

(5) 

8024 

(5) 

7951 

(5) 

9567 

(5) 

8671 

1.2- 

438 

(5) 

2357 

(5) 

2316 

(5) 

2593 

(5) 

3601 

(5) 

3150 

1.78 

650 

(6) 

6623 

(6) 

6425 

(6) 

9611 

(5) 

1542 

(5) 

1312 

2.60 

(6) 

3466 

(6) 

3383 

(6) 

5728 

(6) 

8451 

(6) 

7463 

3-80 

(6) 

2445 

(6) 

2435 

(6) 

3959 

(6) 

4593 

(6) 

4469 

5.58 

(6) 

1942 

(8) 

1957 

(6) 

2Q94 

(6) 

2471 

(6) 

2790 

8.18 

(6) 

1724 

(6) 

1726 

(6) 

2331 

(6) 

1710 

(6) 

2123 

12.0 

(6) 

1579 

(6) 

1562 

(6) 

i960 

(6) 

1444 

(6) 

1821 

17-6 

(6) 

1405 

(6) 

1377 

(6) 

1620 

(6) 

1253 

(6) 

1567 

25.7 

(6) 

1198 

(6) 

1168 

(6) 

1289 

(6) 

1049 

(6) 

1299 

Continued 
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TABLE  8  (Cont'd) 


Decay  of  Normal  Fission  Products  From  U235,  u23®  and  Pu239 
5.  r/hr  at  3  ft  above  an  infinite  plane  for  lCr  fissions  per  sq  ft. 


U235  U238  Pu239 


- 

Years  Days 

Hours 

Thermal 

Fission 

(8  Mev) 

Thermal 

Fission 

0.763 

(8)9977 

(8)9970 

(8)9329 

(8)8907 

(8)8750 

1.12 

(8)6648 

(8)6632 

(8)6172 

(a; 5866 

(8)5Ybl 

1.64 

(8)4149 

(8)4153 

(8)3827 

(8)3592 

(8)3537 

2.40 

(8)2453 

(8)2484 

(8)2256 

(8)2071 

(8)2065 

3-52 

( 8) 1410 

(8)1450 

(8)1303 

(8)1166 

(8)1196 

5.16 

(9)8079 

(9)8418 

(9)7582 

(9)6642 

(9)7098 

7.56 

(9)4786 

(9)5014 

(9)4587 

(9)3986 

(9)4398 

ll.l 

(9)2964 

(9)3094 

(9)2897 

(9)2555 

(9)2821 

16.2 

(9)1804 

(9)1869 

(9)1792 

(9)1626 

(9)1761 

23.8 

(10)9716 

(9)1094 

(9)1073 

(9)1010 

(9)1063 

1.45 

34.8 

(10)6305 

(10)6428 

(10)6393 

(10)6235 

(10)6360 

2.13 

51.1 

(10)3730 

(10)3786 

(10)3817 

(10)3869 

(10)3811 

3.12 

74.9 

(10)2276 

(10)2319 

(10)2365 

(10)2470 

(10)2362 

4.57 

(10)1483 

(10)1524 

(10)1556 

(10)1645 

(10)1546 

6.70 

(11)9986 

(10)1031 

(10)1039 

(10)1099 

(10)1021 

9-82 

(11)6774 

(11)  ^>72 

(11)6899 

(11)7244 

(11)6655 

l4.4 

(11)4490 

(11)4587 

(11)4462 

(11)4650. 

(11)4226 

21.1 

(11)2910 

(11)2940 

(11)2845 

(11)2953 

(11)2660 

30-9 

(11)1813 

(11)1807 

(11)1762 

(11)1837 

(11)1645 

45.3 

(11)1061 

(11)1039 

(11)1034 

(11)1092 

(12)9777 

66.4 

(12)6055 

(12)5807 

(12)5910 

(12)6360 

(12)5728 

97-3 

(12)3676 

(12)3497 

(12)3559 

(12)3896 

(12)3543 

143 

(12)2170 

(12)2090 

(12)2079 

(12)2320 

(12)2200 

■208 

r»J< 

^  J .c.  j  J.X7*t 

/  1  1  1  £)i 

\  J.C./ 

( 1 o\ 1 1 3  3 
\  -*■«-/ 

( 12 ) 12 of 

12 ) 1100 

301 

(13)4874 

(13)4790 

(13)4733 

(13)5707 

(13)5170 

1.2  438 

(13)1399 

(13)1373 

(13)1525 

(13)2135 

(13)1864 

1.78  650 

(14)3884 

(14)3750 

(14)5517 

(14)9083 

(14)7690 

9.60 

(14)2031 

(14)1975 

(14)3160 

( 14)4964 

(14)4352 

3.80 

(l4)l444 

(14)1432 

( 14)2213 

(14)2692 

(14)2594 

5.58 

(14)1154 

(14)1158 

(14)1603 

(].4)l442 

( 14) l6ll 

8.18 

(14)1026 

(14)1021 

(14)1291 

(15)9971 

(14)1225 

12.0 

(15)9432 

(15)9293 

( 14)1094 

(15)8452 

(14)1057 

17.6 

(15)8310 

(15)8211 

(15)9164 

(15)7377 

(15)9160 

25-7 

(15)7183 

(15)6987 

(15)7431 

(15)6219 

(15)7668 
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ACTIVITY  (DOLAN) 
ACTIVITY  (THIS  REPORT) 


TIME  AFTER  FISSION  (HOURS) 


Fig.  3  Comparison  of  14.0  Mev  Neutron  (Dolan29,30) 
to  8  Mev  Neutron  Fission  of  U23o 
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« 


10  100  1000  10,000 
TIME  AFTER  FISSION  (HOURS) 


Fig.  U  Gross  Air  Ionization  Rate  "R"  Factors  for  Various  Types  of 

Fission  (hours  after  fission) 
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1  10  100 

TIME  AFTER  FISSION  (YEARS) 


Fig.  5  Gross  Air  Ionization  Rate  "R"  Factors  for  Various  Types  of 

Fission  (years  after  fission) 
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nuclide  decay  schematics)  are  the  mass  chain  yields  and  the  indep¬ 
endent  yields  of  the  isotopes  in  the  chain.  For  times  after  fis¬ 
sion  of  about  1  hour  and  greater*  the  Bolles -Ballou  calculations-*- 3 
show  that  the  difference  in  the  total  disintegration  rates  based 
on  Present’s  yield  theory  from  those  based  on  Glendenin's  postulate 
is  insignificant.  This  is  due  to  the  fact  that  at  these  times  after 
fission  most  of  the  chains  have  decayed  to  the  last  one  or  two  active 
members  from  short-lived  early  members.  The  displacement  in  the 
curves  of  Figs.  3;  ^ >  and  5  from  1.00  are  therefore  due  to  differ¬ 
ences  in  the  chain  yields .  The  curves  show  that  for  times  between 
about  1  hour  and  1  year*  the  maximum  error  in  the  ionization  rate 
by  use  of  the  data  for  thermal  neutron  fission  of  U235  would  be 
about  15  The  error  in  the  integrated  dose  for  any  time  period 
would  be  less . 

The  H+l  ionization  rates  at  3  feet  above  an  infinite  smooth 
contaminated  plane  for  a  unit  yield  distribution  of  fission  pro¬ 
ducts  per  unit  area  are  summarized  in  Table  9-  The  highest  value 

TABLE  9 

Summary  of  H+l  Ionization-Rates  Per  Unit  Yield  Per  Unit 
Area*  for  3  Feet  Above  an  Infinite  Smooth  Contaminated  Plane 


Type  of  Fission 

H+l  Ionization  Rate 
(Unit  Yield/Unit  Area) 

(r/hr  at  1  hr}/ 
(fiBs/sq.  ft) 

(r/hr  at  1  hr)/ 
(ET/sq.  mi) 

U235  (thermal) 

7.60  x  10"13 

3950 

U235  (fission) 

7.58  x  10_13 

3940 

U23u  (8  Mev) 

6.94  x  10_13 

3610 

Pu239  (thermal) 

6.70  x  10_13 

3480 

Pu239  (fission) 

6.54  x  lO"^ 

3400 

is  for  the  thermal  neutron  fission  of  U235j  the  lowest  is  for 
fission- neutron  fission  of  Pu239.  The  same  value*  1.45  X  10^3 
fissions/KT*  was  used  to  convert  all  the  rates  from  fissions  to 
kilotons-  The  corresponding  ionization  rate  factor  derived  from 
EKW  is  (1240  r/hr  at  1  hr)/(KT/sq.  mi)  or  about  a  factor  of  3  lower 
than  the  values  of  Table  9*  Other  authors  have  made  similar  cal¬ 
culations  and  comparisons  of  these  factors  and  the  decay  curves 
for  the  thermal  neutron  fission  of  u235.31>32*33>3^ 
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SECTION  7 


CALCULATION  OP  THE  DISINTEGRATION  RATE,  AND  AIR  IONIZATION 
RATE  ERuM  FRACTIONATED  EISSTuN  PRODUCTS  OF  U^35 , 


U‘ 


r±iuM  j<j.bs.i.ujw  r-Kuuuuxs  OF  u*--*-,  awn 

CONDENSED  WITHIN  PARTICLES  OF  AN  IDEALIZED  CARRIER 

MATERIAL 


The  properties  of  the  idealized  carrier  material  given  in 
Section  2  were  that  (l)  it  did  not  form  compounds  with  the  fission 
product  elements  or  oxides,  (2)  it  dissolved  them  as  stable  oxides, 

(3)  and  that  it  had  a  melting  point  of  l400°C*  The  last  two  pro¬ 
perties  are  somewhat  similar  to  those  of  some  common  soil  minerals. 
The  selection  of  the  first  tvro  properties  was  made  so  that  a  sample 
calculation  of  the  condensation  process  could  be  made  by  use  of  the 
theory  utilizing  Henry's  law  of  dilute  solutions.  The  process  was 
further  idealized,  by  necessity  due  to  the  fact  that  no  values  of 
Henry's  law  constants  were  directly  available,  by  the  use  of  Raoult's- 
law  of  perfect  solutions;  the  latter  idealization  involved  setting 
Henry's  law  constant  equal  to  the  equilibrium  partial  pressure  of 
the  gaseous  species  of  each  element  over  its  own  oxide. 

The  stipulation  of  the  melting  temperature  is  sufficient  for 
computation  of  the  Raoult's  law  constants  from  the  data  of  Table  1; 
these  are  given  in  Table  10.  The  oxygen  pressure  was  taken  to  be 
1  atmosphere  as  a  substitution  for  the  total  pressure.  The  fact 
that  some  of  the  oxides  show  larger  pressure  than  the  rare  gases 
is  partly  due  to  the  neglect  of  the  difference  in  heat  capacities 
in  the  use  of  the  simplified  vapor  pressure  equation  at  temperatures 
beyond  its  range  of  validity.  It  was  pointed  out  earlier,  however, 
that  errors  due  to  this  over-extrapolation  would  not  influence  the 
fractionation  computation  since  the  amount  condensed  of  the  elements 
with  high  vapor  pressures  would  be  essentially  zero  and  not  sensitive 
to  the  exact  value  of  their  ,va,por  pressures. 

If  no  carrier  material  would  be  present  in  the  fireball  along 
with  the  fission  product  elements,  or  if  they  were  all  in  the  vapor 
state  when  the  temperature  is  l400°C,  their  partial  pressures  can  be 
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TABLE  10 


Summary  of  Raoult’s  Law  Constants  for  the  Oxides  of  Fission 
Product  and  Other  Elements  at  l673°K 


Element  kj  (atmos) 

Element 

kj  (atmos) 

Cu 

7.4  X  io"6 

Pd 

8.7  x  10-7 

Zn 

5.5  X  10-5 

Ag 

2.9  x  10-3 

Go* 

8.5  x  10-4 

Cd 

0.12 

Ge 

0.229 

In 

9.6  x  10"5 

As 

2.6  x  io3 

Sn 

J.2  x  10-5 

Se 

8.7  x  10 ^ 

Sb 

0.339 

Br 

1.05 

Te 

4.17 

Kr 

5.6  x  103 

I 

2.09 

Kb 

i.20 

Xe 

4.9  x  10 3 

Sr 

1-3  x  10-° 

Cs 

3-24 

Y 

6  x  lO"1^ 

Ba 

5.0  X  10-7 

Zr 

2  x  10"14 

La 

2  x  10-15 

Nb 

1.6  x  10-7 

Ce 

8  x  10-13 

Mo 

3-63 

Pr 

1  x  10~14 

Tc 

9.6  x  103 

Wd 

3  x  10-14 

Ru 

1.1  x  10-3 

Pm 

2  x  10-13 

Rh 

9-3  x  10-7 

U 

2.2  x  10-5 

9  6 


estimated  from  the  perfect  gas  law  using  Vm  as  the  volume  and  using 
the  yield  factor  of  0*5  moles  of  fission  products  produced  per  KT- 
The  partial  pressures  are  summarized  in  Table  11  for  some  of  the 
more  abundant  fission  product  elements  as  given  by  Bolles  and  Ballou. ^3 


TABLE  11 


Calculated  Partial  Pressures  and  Equilibrium 
Partial  Pressures  of  Some  of  the  More  Abundant  Fission 
Products  Dispersed  Uniformly  in  Vm  at  l673°K  for  a 
Fission  Yield  of  2*3  MT 


Element 

Percent  No. 
of  Atoms  at 

1  min 

Partial 
Pressure 
in  Volume 
(10"9  atmos) 

Partial  Pressure 
Over  Oxide 
(atmos) 

Rb 

,  5 

4.2 

1.2  (Da 

Cs 

c* 

4.2 

3-2  (1) 

Sr 

8 

6-7 

1-3  x  10-o(s) 

Ba 

7 

5-8 

5.0  x  10“7 (s) 

La 

4  ■ 

3-3 

2  x  10“15(s) 

Ce 

6 

5.0 

8  x  10-13^) 

Sb 

4 

0.84 

0.34  (1) 

Zr 

8 

6.7 

2  x  10-1£|-(s) 

Nb 

6 

2-5 

1.6  x  10_7(s) 

Mo 

8 

6.7 

3-6  (1) 

Te 

6 

5-0 

4.2  (1) 

a.  s  for  solid,  1  for  liquid. 


Of  all  the  elements  listed,  only  La,  Ce,  and  Zr  would  have  condensed 
to  a  solid  state'  at  this  temperature.  The  remainder  would  do  so  only 
at  lower  temperatures .  Use  of  the  equations  of  the  variation  of  the 
gas  volume,  temperature,  and  percent  abundance  of  each  element  with 
time  could  be  made  to  estimate  the  times  that  condensation  would 
begin  for  eac  h.  This  computation,  if  other  vaporized  materials 
(bomb  casing,  soil,  etc.)  were  included,  could  be  made  to  obtain 
more  information  about  the  formation  of  vapor-condensed  particles 
in  the  early  stages  of  the  whole  fallout  formation  process. 

In  computing  the  amount  of  each  mass  chain  that  condensed  in 
the  liquid  soil,  use  was  made  of  Eq.  29  to  calculate  the  condensation 
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"R"  factor  with  reference  to  the  number  of  fissions.  The  computation 
of  r0(A)  by  use  of  the  Bolles -Bailout  values  of  the  independent  yields 
for  both  Present's  and  Glendenin fs  yield  theories  is  illustrated  for 
masses  89 ,  90,  and  l4o  in  Table  12;  the  value  of  k<?  was  determined 
by  dividing  the  values  of  kj  in  Table  10  by  0.0364 .  For  nuclides 
other  than  the  end  numbers  of  the  chain,  the  accumulated  sum  of  the 
y^(A)/(l+k9)  terms  up  to  the  indicated  chain  number,  divided  by  the 
accumulated  percent  yield  of  the  chain,  gives  the  r0(A)  value  for 
that  nuclide.  This  converts  the  partial  sums  to  fractional  multi¬ 
pliers  for  each  nuclide  in  the  chain. 

The  most  exact  method  of  computing  the  amount  of  each  radio¬ 
nuclide  present  (and  its  activity)  for  a  given  number  of  fissions 
would  he  to  use  the  value  of  yj(A)/(l+k9)  for  each  and  the  approp¬ 
riate  decay  formula  for  the  production  of  the  daughter  products. 

But  to  simplify  the  computation  process,  the  calculations  were  made 
by  direct  multiplication  of  the  r„(A)  values  and  the  single  nuclide 
d/s  values  per  10^  fissions  and  r/hr  values  per  10^  fissions/sq  ft. 
Although  this  procedure  gives  higher  values  for  the  daughter  products 
in  each  chain,  the  error  decreases  with  time.  Compared  to  other 
possible  errors  involved  in  the  computation,  the  error  due  to  this 
approximation  is  rather  insignificant.  For  example,  the  initial 
fraction  for  Sr$9  ie  0.010  which  increases,  after  Rb®9  decays  out, 
to  O.OI83;  in  the  exact  method  of  computation,  the  latter  value  is 
reached  by  2.40  hr.  The  ro(90)  value  for  Sr90  i3  valid  from  31.2 
minutes  and  the  ro(l40)  value  for  Ba-^0  is  essentially  valid  at 
60  seconds. 

The  rQ(A)  values  for  all  the  radionuclides  used  in  the  computa¬ 
tions  are  summarized  in  Table  13  for  Loth  the  Present  and  Glendenin 
independent  yield  theories  as  calculated  by  Bolles  and  Ballou. 13 
For  more  than  half  of  the  nuclides,  the  r  (A)  values  from  the  Glendenin 
distributions  are  larger.  The  computed  r(A)  values  relative  to  mass 
99  for  some  nuclides  of  interest  are:  Sr&9,  0.0232  and  0.0200;  Sr90, 
O.195  and  0.193;  I131,  0.914-  and  0.108;  Csl37,  0.00408  and  0.00455; 
and  Ba-1-1^,  0.517  and  0-574.  The  rQ(A)  values  indicate  that  about 
85  percent  of  the  total  Sr90  formed  (for  the  2-3  MT  yield  surface 
burst)  would  be  condensed  on  (small)  solid  particles  during  the 
second  period  of  condensation  along  with  about  99  percent  of  the 
Csl37  and  its  daughter,  Bal37m,  some  fraction  of  these  amounts 
would  presumably  be  soluble  and  biologically  available  after  the 
particle  reached  the  earth's  surface. 

The  results  of  the  decay  computations  in  d/s  per  10^  fissions 
and  r/hr  per  104  fissiens/sq  ft  are  given  in  Table  14.  The  ioniza¬ 
tion  rates  for  the  Glendenin  yields  are  plotted  in  Fig-  6,  with  the 
curve  for  normal  fission  products  from  thermal  fission  of  U235.  The 
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TABLE  13 


Suwnsry  of  rQ(A)  Values  *t  60  Seconds  After  Fission  for  Fission  Product  Nuclides  That  contribute 
to  the  Gross  Activity  at  45.8  Minutes  After  Fission 


Nuclide 


ZnT® 

Zd34 

c.32 


^5 

*73 

„78 


i*36 

..79 


an01 

s<=a87 

66^3 

Br63 


Krp03 

Kra85 
Kr 87 
*•88 

Hb88 

Bb® 

Bb^7 

Sr89 

Sr9° 

Sr91 

Sr32 

3,33 

YS» 

Yl9i 

x29? 

y92 

y93 

y94 


Zr93 

■*97 


o<*> 


Present  Glendenin 


0.996 

0.998 

0-997 

0.998 

0-997 

0.998 

AnOc 

0-986 

0.990 

0.993 

0.600 

0.61*9 

0.183 

O.203 

O.lfcO 

0.142 

0.137 

0.159 

0.126 ' 

0. 13$ 

0.0226 

0.0288 

0.000315 

0.000800 

0.00Oll»3 

0. 000419 

0.0 

0.0 

0.00926 

0.00808 

0.0135 

0.0135 

0.00926 

0.00806 

0.0279 

0.0285 

0.0266 

0.0275 

0.0129 

0.0152 

0.00122 

0.00171 

0.00213 

0.00213 

0.00840 

0-00735 

0.0291 

0.0290 

0.0183 

0.0153 

0.154 

0-148 

O.236 

0.243 

0-921 

0.922 

0.999 

0.998 

0.154 

0.148 

O.236 

0.2U3 

0.236 

0.243 

0-921 

0.923 

0,-999 

0.998 

1:00 

1-00 

1.00 

1-00 

1.00 

1-00 

Nuclide  r=(A'  Nuclide  r0^' 


Present 

Glendenin 

Present 

Glendenin 

Nbi95 

1-00 

1.00 

5b123 

0.998 

0.998 

Nb^S 

1.00 

1.00 

Sb126 

0.967 

0-975 

«b37 

1.00 

1.00 

Sb723 

O.902 

0.947 

Hb38 

1.00 

1.00 

Sb128 

0.109 

o.m 

Mo99 

0.789 

0.765 

Sb7® 

0-0973 

0.0974 

no101 

0.0396 

0.0410 

sb13* 

0.0970 

0.0970 

Mo7“ 

0.0141 

0.0121 

t.^5 

0.998 

0.998 

Tcp39 

0-709 

0.765 

Te^l27 

0.902 

0.947 

TclOl 

0.0383 

0.0293  .  - 

7*2^ 

0.902 

0.947 

Tc102 

0.0135 

0-0112 

Te1129 

0.0958 

0.0968 

Ru703 

Ru7°5 

Ru706 

0.266 

0.869 

0.278 

0.889 

Te2129 

TCj731 

Td2131 

0.0941 

0 .0833 

0-0724 

0.0960 

0.0902 

0.0829 

0.955 

0.960 

t.132 

0.01)26 

0.0563 

Shi103 

0.266 

0-278 

Tep33 

0 .0448 

0.0654 

ShjW 

0.869 

0.890 
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2.  r/hr  at  3  ft  above  an  infinite  plane  for  101*  fissions  per  sq.  ft 
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Fig.  6  Air  Ionization  Rate  at  3  ft  Above  a  Smooth  Infinite 
Plane  Uniformly  Contaminated  With  Fractional  Fission  Pro¬ 
ducts  From  10^  Fissions/sq  ft,  in  Melted  Fallout  Particles 
From  a  2.3  MT  Yield  Surface  Detonation.  Rates  for  normal 
fission  products  from  thermal  neutron  fission  of  U235  arc 
given  for  comparison, 
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AIR  IONIZATION  RATE  (R/HR) 
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gross  fission  product  "R"  factor,  r^p,  which  is  defined  as  the  ratio 
of  the  ionization  rate  for  the  fractionated  products  to  that  for 
the  unfractionated  products  from  thermal  fission  of  IJ235  is  given 
in  Fig.  7-  Up  to  about  5000  hrs  (200  days),  the  differences  in  the 
rfp  values  for  the  different  types  of  fission  are  not  large.  At 
1  hr  after  fission,  the  rfp  values  for  U23°  (8  Mev),  U235  (fission) 
range  from  Q. 37  to  0.40.  The  minimum  at  about  150  hrs  is  where 
I131  and  Ba^O-LalW  are  in  high  abundance  and  the  maximum  at 
about  3500  hrs  is  where  Zr95-Nb95  are  in  high  abundance.  The 
large  r^p  values  after  10,000  hrs  for  Pu239  (fission)  is  due  to 
the  larger  yields  of  some  of  the  rare  earth  nuclides. 

The  decay  curves  were  calculated  only  for  radionuclides  dissol¬ 
ved  in  liquid  soil  particles  and  therefore  apply  only  to  a  group  of 
particles  in  the  fallout  of  a  given  size  range  (say,  those  that 
separated  from  the  fireball  between  about  50  and  70  seconds  after 
burst) .  The  rfp  values  for  smaller  particles  would  generally  be 
larger  and  even  approach  values  that  are  the  inverse  of  those 
shown  in  Fig.  7»  The  rfp  values  for  the  world-wide  fallout  might 
approach  1  minus  the  rfp  values  of  Fig.  7 • 

The  differences  in  the  d/s  and  r/hr  rates  between  those  based 
on  Present's  theory  of  yield  and  Glendenin's  are  insignificant. 

The  calculations  for  both  sets  were  made  because,  when  the  conden¬ 
sation'1  process  is  considered,  the  later  decay  rates  are  very  much 
dependent  on  the  independent  nuclide  yields  present  at  the  time 
of  cessation  of  the  first  phase  of  the  condensation  process.  Thus 
it  appears  that,  for  yields  in  the  MT  yield  range  and  for  soil 
melting  at  temperatures  below  2000°K,  the  gross  decay  of  the  frac¬ 
tionated  fissions  is  not  sensitive  to , the  differences  in  the  inde¬ 
pendent  nuclide  yields  from  the  two  theories.  It  would  be  expected, 
however,  that  the  differences  in  the  gross  decay  curves  would  in¬ 
crease  somewhat  as  the  yield  decreased  and  the  soil  melting  point 
increased. 

The  H+l  air  ionization  rate  in  units  of  KT's  per  square  mile 
for  the  various  type3  of  fission  for  the  s2-3  MT  surface  burst  are 
listed  in  Table  15 •  These  values  cannot  ;be  directly  applied  to  the 
fraction  of  a  weapon  per  unit  area  that  produces  a  given  radiation 
intensity  as  would  be  measured  with  a  survey  meter  three  feet  above 
an  extended  flat  contaminated  area.  First,  some  induced  activities 
will  be  produced;  if  only  the  production  of  u239=jjp239  is  considered, 35 
a  yield  of  1  atom  U^39  pgr  fission  would  increase  the  above  values 
by  Q.3U  X  10"13  (r/hr)/fission/sq  ft)  or  180  (r/hr)/(KT/sq  mi) . 

Second,  a  Co^Q  calibrated  survey  instrument,  with  operator,  will 
give  a  meter  reading  at  3  ft  above  the  surface  about  25  /  lower 
than  the  true  air  ionization  rate. 28  Third,  a  real  surface  is  not 
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1  10  100  1000  10,000 

TIME  AFTER  FISSION  (HOURS) 


Fig.  7  Variation  of  the  Gross  Ionization  Rate  "R"  Factor  for  the 
Fractionated  Fission  Products  With  Him  After  Fission 
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TABLE  15 


Air  Ionization  3  ft  Above  the  Surface  at  H+l  for  Unit 
Yield  Fallout  Distributions  on  an  Ideal  Plane,  for 
Fractionated  Fission  Products  From  a  2*3  MT 
Surface  Detonation 


Type  of  Fission 


Air  Ionization 
(Unit  Y ield'/tTnit  Area) 


kr/hx-  at  1  hr)/ 
(fiss/sq  ft)a 


kr/hr  at  1  hr)/ 
(KT/sq  mi)a 


U235,  thermal 
neutron 
U235,  fission 
neutrons 
U23&,  8  Mev 
neutrons 
Pu239,  fission 
neutrons 


3.00  x  10“13 
2.96  x  lO"^ 
2.85  x  10“13 
2.92  x  10-13 


1,560 

1,540 

1,480 


1,520 


a .  Glendenin.'  s  Postulate . 


smooth  so  that  actual  air  ionization  at  3  ft  above  the  surface  will 
again  be  less  than  that  given  in  Table  15 •  If  the  terrain  rough¬ 
ness  .  vejx-  ?Uch  as  to  reduce  the  ionization  rate  by  20  56,  the  com¬ 
bination  of  the  three  factors  would  give  about  1000  (r/hr)/\  KT/sq  mi 
as  the  "observed"  value  of  the  ratio  for  8  Mev  neutron  fission  of 
U238.  unfortunately,  this  conversion  factor  between  r/hr  and  KT/ 
sq  mi  will  not  be  constant  over  the  whole  region  of  heavy  fallout. 

At  further  downwind  distances  where  the  fractionation  would  be  less, 
the  observed  conversion  factor  may  approach  or  even  exceed  the 
value,  (36IO  +  180)  X  0.6,  or  2270  (r/hr) /(KT/sq  mi)  for  the  U23& 

(8  Mev)  fission.  Even  with  fractionation,  the  values  in  Table  15 
are  higher  than  the  factor  derived  from  EEWl4  (1240  r/hr/(KT/sq  mi)) 
However  no  specifications  as  to  yield  and  type  of  detonation  are 
given  for  the  information  in  ENW;  hence  the  value  from  that  refer¬ 
ence  is  more  appropriately  associated  with  the  data  of  Section  6 
even  though  its  numerical  value  is  much  nearer  to  those  given  here- 

The  factors  given  in  Table  15  do  not  have  the  precise  relation 
to  fissions  or  KT's  the  way  the  factors  given  in  Section  6  for  nor¬ 
mal  fission  products  do.  'The  original  mixture  of  fission  products 
which  is  directly  related  to  the  yield  has  been  altered  in  many 
ways.  It  is  possible  to  compute  the  values  in  Table  15  because 
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the  fraction  of  each  nuclide  produced  and  condensed  has  been  calcu¬ 
lated  on  the  basis  of  the  original  fission  yield.  In  an  experimental 
evaluation  of  the  factor,  the  yield  representation  of  a  square  foot 
of  fallout  contamination  is  unknown .  A  representation  can  be  made 
to  some  degree  6r  accuracy  by  a  radiochemical  analysis  of  the  fall¬ 
out  material  and  a  measurement  of  the  air  ionization  rate  (corrected 
to  H+l  or  other  reference  time) .  The  selection  of  the  radionuclides 
for  analysis  should  depend  on  some  knowledge  of  the  type  of  particles 
present.  For  the  types  of  particles  represented  by  the  calculations 
given  here,  an  analysis  that  gave  the  abundance  of  the  radionuclides 
with  r0(A)  values  of  1.00  (see  Table  13)  should  give  the  same  value 
of  the  factor  as  discussed  for  those  in  Table  1$.  It  may  be  noted 
that  Mo99,  for  the  conditions  of  these  calculations,  would  give 
factors  that  are  27  to  31  $>  high.  In  the  real  fallout  where  a 
mixture  of  particle  sizes  and  types  are  present,  it  would  appear 
that  the  selection  of  the  appropriate  radionuclides  and  their 
yield  representation  should  depend  on  other  factors  than  the 
originalms.es  chain  yield  and  its  chain  decay  characteristics. 


Approved  by: 


EUGENE  ?/.  COOPER  / 
Associate  Scientific  Director 
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SEP  3  0  I960 

From:  Commanding  Officer  and  Director 

To:  Dis  tribution  for  Technical  Report  USNBDL-TR-425 

Subj :  U.  3.  Naval  Radiological  Defense  Laboratory  Report  USNRDL-TR-425 

en  titled  *1  Theory  of  Formation  of  Fallout  from  Land-Surface  Nuclear 

Detonations  and  Decay  of  the  Fission  Products"  by  C.  F.  Miller; 

_ _ 

CA-  tAJUtt  AUA 

1.  It  ic  requested  that  the  following  corrections  be  made  in  subject  report: 

p.  iv  -  last  paragraph,  first  sentence,  add  "are  presented*  after 
■fission  of  Pu^39.» 

p.  11  -  Equation  11  is  to  read  as  follows: 

p.  15  -  Equation  21  is  to  read  as  follows: 

p.  27  -  Equation  70  is  to  read  as  follows: 


KFinp/p  =  42 25 
_o 

J  (n  (i)/V  )ET 


n 1  = 


jl  1  +  kj 


-^4 

RT 


p.  32 

P.  45 
p.  46 

p.  49 

P.  50 


line  3  reads:  "Combining  Eqs.  57  and  58  gives" 

Change  to  read  "Combining  Eqs.  2*  and  21  gives* 

Una  1  -  correct  typographical  error  in  spelling  of  word 
"expansion* 

line  8  reads  "Substituting  Eqs.  103,  104,  and" 

Change  to  read  "Substituting  Eqs.  103,  105.  and" 
line  23  reads  "then  T2  from  Eq.  112  is" 

Change  to  read  "then  T£  from  Eq.  121  is* 

Equation  131  closes  as  follows:  /"tf 

/  -1.490 

t  dt 

J  ** 


1/  rT- 


OCT  5 
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Subj :  U.  S.  Navsl  Radiological  Defense  Laboratory  Report  TJSN3.V.  L-TR-425 

entitled  "JL  Theory  of  Formation  of  Fallout  from  Land-Su"  face  Racism: 
Detonations  and  Decay  of  the  Fission  Products*  by  C.  F,  Miller; 
errata  for 


Change  this  portion  of  the  equation  to  read  as  follows t 


-1.5.90 


p.  58  -  footnote  **  reads  •  •  .  .  j  determined  from  Eq„  121b* 

M 

Change  to  read  *....$  determined  from  Eq.  149b* 
p.  63  -  line  14  reads  *of  Fdt  for  Eq.  176,  is* 

Change  to  read  *of  Pdt  for  Eq.  177.  is* 
p.  64  -  last  line,  add  *so  that*  after  *per  molecule* 
p.  66  -  line  14  reads  *  .  .  .  for  other  yields  Eq.  178* 

Change  to  read  *....;  for  other  yields  Eq.  179" 
p.  98  -  third  paragraph,  seventh  line  reads  *  .  .  . 

0.914  and  0.108* 

Change  to  read  "...  ;  j?»  0.0914  and  O.lOf!;*. 


T,  J.  MATHERS 
By  direction 


